Earth's geochronology map: a work in progress...

Tips for creating and manipulating planet textures for Celestia.
Avatar
Topic author
Fenerit M
Posts: 1880
Joined: 26.03.2007
Age: 17
With us: 17 years 11 months
Location: Thyrrenian sea

Earth's geochronology map: a work in progress...

Post #1by Fenerit » 03.08.2010, 23:35

Having verified its feasibility, I open this thread. The topic will shows progresses, problems and various issues. Premise that I do not have ARC(something) GIS programs, all that you can see here is made through free softwares, graphics and GIS, with a bit of personal Perl scripts. The map show the geochronology of the Earth's ROCKS.
The map is based upon free ESRI shapefile vectors hosted here:

http://nationalatlas.gov/atlasftp.html#geology
for U.S. conterminous and the rest;

http://geogratis.gc.ca/geogratis/en/option/select.do;jsessionid=2B48BF4C258172E4AED10E6566E85417?id=E9627B13-9ED3-4100-EEA1-16817F1F36FC
for Canada (despite its title, the geo ages are within the DBF)

http://certmapper.cr.usgs.gov/rooms/we/index.jsp
for the rest of the world but the western Mexico. In the latter case I don't know where to find one, so I hold for good the one from here (after vectorized):
http://www.lib.utexas.edu/maps/mexico.html

This is all of the Earth's geochrons maps I got for free.

The open source GIS software with which I deal with shapes are
http://www.mapwindow.org/ and Quantum GIS http://www.qgis.org/, but mostly the former.

Shapes have several projections, thus all are first converted in geographic projection (WGS 1984, ellipsoid) with Mapwindow, and then overlaid on its long/lat reference grid. Since the maximum graphic output's resolution is 8k, such grid have been split in four part in order to get, at least for VT, 16k; each part has been processed onto GIMP with the antialias filter and re-joined. The 16k map successively has been tiled with the F-TexTools:
http://www.celestialmatters.org/?q=node/77

Now, like I said, the success of the operations above are the roots to assert that the geochronology map is feasible. The problems are others:
1) to uniform the age-colors;
2) the CMOD models;

Hiterto some screenshots of the partial map:

Image

Limit of level 3. For solving the first problem I adopt the color convention of US map for sedimentary, metamorphic, volcanic and intrusive ROCKS. Being the maps diversely colored amongst them, the US dataset is the only dataset whose DBF records shows the RGB colors values. The Canada's map is very detailled and exaustive but doesn't show colors. Bear in mind that Mapwindows doesn't import the ESRI exchange file .e00 and Qgis often fails in read it correctly. Nonetheless, I do manually assign them to the subshapes.

Image

level 4. By the manual assignment of color to the subshapes (shapes made from each eon/era/epoch/period) is possible the export also them as image for being successively "patternized" (I do not like the font choosed for "v" of volcanic. I may find another more sharpen. Instead the "+" of intrusives in near to be fine).

Image

Naturally, level 5 should be still better and in theory would be possible with yet another subdivision of the lat/long grid (4x4) and then to work on 16 pieces of the map... As experiment, on my machine GIMP does riassemble 16 x 8k pieces (32k) of black background in about 30 min. However, Photoshop isn't able to handle 32k images at all. Too bad for such an expensive software.

For what concern the "CMOD problem", is that one CMOD model only is BIGGER, even in binary, and the possible solution should be that of splitting the model in pieces, one for each continents; with the relevant scripts to make it active. CMOD models are yielded through an enhanced version of the Perl script shipped with CMODtracer and fitted to the Zonums http://www.zonums.com/shp2text.html shp2text converter.
So far, just an introduction.
Last edited by Fenerit on 25.08.2010, 09:19, edited 1 time in total.
Never at rest.
Massimo

Avatar
John Van Vliet
Posts: 2949
Joined: 28.08.2002
With us: 22 years 6 months

Re: Earth's geochronology map: a work in progress...

Post #2by John Van Vliet » 04.08.2010, 09:13

--- edit ---
Last edited by John Van Vliet on 19.10.2013, 23:13, edited 1 time in total.

Avatar
Topic author
Fenerit M
Posts: 1880
Joined: 26.03.2007
Age: 17
With us: 17 years 11 months
Location: Thyrrenian sea

Re: Earth's geochronology map: a work in progress...

Post #3by Fenerit » 04.08.2010, 21:52

john Van Vliet wrote:looks a lot like the USGS PIGWAD
http://webgis.wr.usgs.gov/pigwad/maps/index.html

Well, I do not find the Earth from the link above, :roll: but yes, I know several links hosted. (paywall, though). The map (as well as the others) are here: http://certmapper.cr.usgs.gov/rooms/we/index.jsp (mind the combo box below the interactive window :wink: there are the options for datasets not directly linkables from here).

Australia is here:
http://pubs.usgs.gov/of/1997/ofr-97-470/OF97-470F/aspac3.html#TOP but colored accordingly with U.S. geologic map. Apart the shapes, these maps are inside the PDF's.
Never at rest.
Massimo

Avatar
Topic author
Fenerit M
Posts: 1880
Joined: 26.03.2007
Age: 17
With us: 17 years 11 months
Location: Thyrrenian sea

Re: Earth's geochronology map: a work in progress...

Post #4by Fenerit » 04.08.2010, 21:56

One issue with the map. The image below show the shapefile of the former Soviet Union (CCCP) as imported in Mapwindow with its default equidistant conic projection.

Image

Now, when the shape is reprojected as WGS 1984, the lands which crosses the last meridian (180 deg. longitude) are stretched throughout their latitudes.

Image

Image

Diversely from SVG, whose vectors in an hypotetic conversion of the shape should be single polylines, here they are integrant parts of the polygons; and cannot be deleted without to delete yet the whole polygon. To delete the polygons makes holes in the far eastern/western parts of shape; as conseguence, several pieces of it will lacks both of geologic ages in graphical form and CMOD as well.
In such case, instead of to guess following the manual method of depict the likely boundaries of the rock beds, I think should be more appropriate to mark the areas like NaN (Not available Number). There are areas already marked as "undetermined age" for the various shapes, thus it shouldn't a great disrespectfulness. For the rest, the CCCP shapefile match very well with the neighborhood shapes, except for smallest areas along the boundaries with the South Asia, which can be either joined manually (CMOD and map) or left as well as "NaN".

Note that such a map has been scansioned and vectorized in the 1998 from an original soviet map of the 1966; while the areas outside it are somewhat recent works in Lambert Conformal Conic projection. Thus, summing up that soviets maybe were unawares of the precise locations of their boundaries along the mountains of the central/south Asia, with the paper scanning's method and with the different projections, such map is a good piece of the puzzle.
Never at rest.
Massimo

Avatar
Topic author
Fenerit M
Posts: 1880
Joined: 26.03.2007
Age: 17
With us: 17 years 11 months
Location: Thyrrenian sea

Re: Earth's geochronology map: a work in progress...

Post #5by Fenerit » 06.08.2010, 22:27

The "colors' issue".

Image

The sharp line amongst orange and reddish, is the rocks' boundary amongst Alaska and Canada (white areas are zones not yet colored). For the U.S shape, such intrusive granitic rocks have an age spanning Cretaceous-Tertiary (KTg) while for the Canadian shape such rocks are well identified as intrusive granitic of Paleogene (Tertiary as well); that is, they would have to be colored in orange (pgTg). Of course are the same rocks with the same ages, being ROCKS not patriotics, but the colors' assignment is different for the U.S legend. The Canada's map is beautiful in its rock's differentiation, but here will be "reduced" for uniformity purposes; as well as the remains world' geologic shapes have been enhanced accordingly with U.S.
Never at rest.
Massimo

Avatar
Topic author
Fenerit M
Posts: 1880
Joined: 26.03.2007
Age: 17
With us: 17 years 11 months
Location: Thyrrenian sea

Re: Earth's geochronology map: a work in progress...

Post #6by Fenerit » 08.08.2010, 22:41

Just two screenshots...

Far East...
Image

The origin of the psychedelic rock...
Image
Never at rest.
Massimo

Avatar
Topic author
Fenerit M
Posts: 1880
Joined: 26.03.2007
Age: 17
With us: 17 years 11 months
Location: Thyrrenian sea

Re: Earth's geochronology map: a work in progress...

Post #7by Fenerit » 10.08.2010, 00:08

The "CCCP stretching" problem has been solved for an half. One CCCP's WGS shape model have had the polygons deleted, in order to make the map without horizontal "stripes" and in which there are some blank areas (to be named like NaN, geology not shown, LOK etc.); while another CCCP's WGS shape model has been left stretched. When converted in CMOD the polygons' positions are restored and the model is complete.
Never at rest.
Massimo

Avatar
Topic author
Fenerit M
Posts: 1880
Joined: 26.03.2007
Age: 17
With us: 17 years 11 months
Location: Thyrrenian sea

Re: Earth's geochronology map: a work in progress...

Post #8by Fenerit » 10.08.2010, 21:45

Solved even the second half.
- The stretched shapes have been saved apart.
- Converted in SVG with an arrays/paragraphs version of the Perl script shipped with CMODtracer
- Imported the SVG in Inkscape and removed the straight lines.

Image


Image

TODO:
- To export the SVG as layer upon the map and to fill the areas with the proper colors seekables through the SHAPE__ID...

P.S.
Probably, isn't necessary to split furthermore the reference long/lat grid in order to achieve deep VT level; being the output a bitmap with color's blocks, the rescale without interpolations get always the same result. Then is matter to find the best antialias setting for the last level (at this point, possibly the fifth; beyond my machine goes overburn).
Never at rest.
Massimo

Avatar
Topic author
Fenerit M
Posts: 1880
Joined: 26.03.2007
Age: 17
With us: 17 years 11 months
Location: Thyrrenian sea

Re: Earth's geochronology map: a work in progress...

Post #9by Fenerit » 11.08.2010, 23:40

The images below shows some steps in georefencing's process of the geologic map of Mexico made with Quantum Gis (Qgis).

Image

The original image map in GIMP.

Image

The map polished, saturated and clipped to the layer.

Image

The map imported inside the georeferencing tool with shared sampled points drew up from the Mexico shape. The Mexico shape has been extracted from the world country dataset. The georeferencing tools is very friendly and allow such operation. Note that Qgis is an open source software available also for Linux and Mac.

Image

The image above show what come up with only 11 sampled points and spline interpolation.
Never at rest.
Massimo

Avatar
Topic author
Fenerit M
Posts: 1880
Joined: 26.03.2007
Age: 17
With us: 17 years 11 months
Location: Thyrrenian sea

Re: Earth's geochronology map: a work in progress...

Post #10by Fenerit » 12.08.2010, 15:40

The final GeoTIFF after spline transformation and interpolation...

Image

For a different test, the GeoTIFF how to look once imported in Mapwindow as missing piece amongst U.S. and Caribbean shapes.

Image
Never at rest.
Massimo

Avatar
Topic author
Fenerit M
Posts: 1880
Joined: 26.03.2007
Age: 17
With us: 17 years 11 months
Location: Thyrrenian sea

Re: Earth's geochronology map: a work in progress...

Post #11by Fenerit » 24.08.2010, 13:13

The code below is the Perl script with which I convert in CMOD the USGS and rest of the shapefiles after these latters before have been converted in TXT by the Zonums http://www.zonums.com/shp2text.html Shape2Text converter. Such (Windows) free converter get an output of the shape in this manner:

Code: Select all

shape2Text
Powered by Zonum Solutions
www.zonums.com

Entity type:Polygon

XCoo YCoo
94  Nodes
44.5877863223 39.7880631875
44.5923726216 39.7922647793
44.5969591894 39.7964705533
44.6040850954 39.8055222932
44.6103685864 39.8050787771
.
.
.

114  Nodes
44.5877863224 39.7880631874
44.6010694543 39.7793226069
44.6088577349 39.7760032622
44.6094029729 39.7702275558
44.6127868432 39.7607786752
.
.
.


The Perl script process such structure. Note that I'm not a Perl genius and the way in which I solved the CMOD structure's issues could be found primitive, expecially in the parts of the temporary reverses process in order to truncate the file. This can do the script a bit slow. The Canada shape, which its TXT conversion get a 60 Mb file, is processed on my machine in about 9 min. For such concerns, I've introduced some messages of procedure advancing in order to display that the script is working, otherwise could seem the script be "hunged up" (whether one doesn't has system meters of any sort). Nonetheless, neither Perl nor the CMODview utility and Celestia had fails in dealt with such resulting CMODs! And all the USGS and the remains shapes are now ready. (BTW, also the map, just some retouches).


Code: Select all

#    Author Fenerit <fenerit@interfree.it>
#    Version 1.0, 23.01.2010
#   Zonums <www.zonums.com> Shape2Text converter utility to CMOD
#    This script convert in ASCII CMOD format a set of geographic coordinates yielded in TXT form
#   through the Zonums Shape2Text converter.
#   The resulting 3d CMOD model is assumed to be used in Celestia, the 3D space simulator.
#   Acknowledgements: Dr. Fridger Schrempp <www.celestialmatters.org> for its basic routines.

#   USAGE:
#   1) Rename or save the TXT file yielded by the Zonums program as "shapefile.txt"
#   2) Script and coordinates text files must be in the same folder;
#    2) Double click on its icon, or type within console:
#   "shape2cmod.pl < shapefile.txt" (in Windows cmd: "perl shape2cmod.pl < shapefile.txt" then press ENTER)
#   A new file called "shape2cmod.cmod" will be created in the same folder;

my $radius = 1;
my $pi = 3.14159265359;
   
   print STDOUT "\n -- Text shapefile to CMOD converter. --\n\n -- Start polygons processing... Please wait...\n";
   
open(TXT, " < shapefile.txt") || die "Can not read the file!"; 
open(TMP, " > conv.tmp") || die "Can not write the file!";

{           

   @G_polygon = <TXT>;
   
}
   
   foreach $_(@G_polygon)

   {
   
       next if (/^$/);   
      next if (/^\s*shape2Text\s*$/);
      next if (/^\s*Powered by Zonum Solutions\s*$/);
      next if (/^\s*www.zonums.com\s*$/);
#      Change below "Polygon" to "Polyline" or "Point" accordingly with the SHP type (default: Polygon)
      next if (/^\s*Entity type:Polygon\s*$/);
      next if (/^\s*XCoo YCoo\s*$/);   
      
      chomp $_;
   
      ($num,$name) = split(" ",$_);
         
      my $vertexcount = $num;
         
      if($name == "Nodes") {
            
         print TMP "\nvertices $vertexcount\n";
               
      }
   
      else
             
      {
      
         $theta = $pi/180.0 * $name; # Polar angle $theta in radiants
         $phi = $pi/180.0 * $num;    # Azimuthal angle $phi in radiants

         $x = -($radius * cos($theta) * cos($phi));
         $y = $radius * sin($theta);
         $z = $radius * cos($theta) * sin($phi);
                  
         printf TMP "%8.6f %8.6f %8.6f\n", $x,$y,$z;
         
      }
   
   }
   
close TXT;
close TMP;

   sleep 1;
   
   print STDOUT "\n\n -- 25% done...\n";
   
   sleep 1;
   
open(TMP, " +< conv.tmp") || die "Can not read and write the file!";

{           

   local $/ = '';
   @S_polygon = <TMP>;
   
   print TMP "\nend__celmodel__ascii\n";

}

   foreach $_(@S_polygon)

   {
   
      next if (/^$/);   

      chomp $_;
   
      ($name,$num) = split(" ",$_);
      
      my $vertexcount = $num;
         
      if($name == "vertices")   {
      
         print TMP "\nmesh\n";
         print TMP "vertexdesc\n";
         print TMP "position f3\n";
         print TMP "end_vertexdesc\n\n";
         print TMP "$_";
         print TMP "\nlinestrip 1 $vertexcount\n";
      
         for ($index=0; $index < $vertexcount; $index++)
            
         {
            printf  TMP "%1d\n", $index;
         }
            
      print TMP "end_mesh\n\n";   
            
      }
   
   }

close TMP;

   sleep 1;
   
   print STDOUT "\n\n -- 50% done...\n";
   
   sleep 1;

open(TMP, " < conv.tmp") || die "Can not read the file!"; 
open(CTMP, " > shape2cmod.tmp") || die "Can not write the file!";

   @R_CMODascii = <TMP>;

    while ($_ = pop @R_CMODascii) {
       
      last if (/^end__celmodel__ascii$/);
      print CTMP $_;
   
   }
           
close TMP;
close CTMP;

   sleep 1;
   
   print STDOUT "\n\n -- 100% done. Parsing CMOD...\n";
   
   sleep 1;

open(CTMP, " < shape2cmod.tmp") || die "Can not read the file!"; 
open(CMOD, " > shape2cmod.cmod") || die "Can not write the file!";
   
   CMODplaceholder();
   
   sleep 1;
   
   @CMODascii = <CTMP>;

    while ($_ = pop @CMODascii) {
       
      print CMOD $_;
   
   }
   
close CTMP;
close CMOD;

   print STDOUT "\n *** Finished! ***\n\n", " -- Use CMODview for binary conversion.\n";

   sleep 2;
   
   unlink ("conv.tmp"); # delete the TMP file
   unlink ("shape2cmod.tmp"); # delete the TMP file

sub CMODplaceholder # placeholder constructor; set only where is ">" for your custom value
   
   {

   print CMOD "#celmodel__ascii\n";
   print CMOD "\n";
   
   print CMOD "material\n";
   print CMOD "diffuse 0 0 0\n";
   print CMOD "opacity 0\n";
   print CMOD "end_material\n";
   print CMOD "\n";
   
   print CMOD "material\n";
   print CMOD "diffuse 0 1 0\n"; # > for the color
#   print CMOD "emissive 1 1 1\n"; # > for its emissivity, uncomment first
   print CMOD "opacity 1\n"; # > for its opacity
   print CMOD "end_material\n";
   print CMOD "\n";
   
   print CMOD "mesh\n";
   print CMOD "vertexdesc\n";
   print CMOD "position f3\n";
   print CMOD "end_vertexdesc\n";
   print CMOD "\n";
   
   print CMOD "vertices 6\n";
   print CMOD "1 0 0\n";
   print CMOD "-1 0 0\n";
   print CMOD "0 1 0\n";
   print CMOD "0 -1 0\n";
   print CMOD "0 0 1\n";
   print CMOD "0 0 -1\n";
   print CMOD "\n";
   
   print CMOD "linestrip 0 6\n";
   print CMOD "0 1 2 3 4 5\n";
   print CMOD "end_mesh\n";
   print CMOD "\n";

   }


Copy the code in an empty document, and save the file as "yourname.pl".

P.S.
The Zonums Shape2Text converter although be a free program, is however limited as "expireware", in the sense that when one does launch the program, a message box show the expired warn and the utility doesn't start.
You must set the system clock year back to 2009 first, before to start the program. In such a way the utility is fully functional, even whether the date is restored forth after its loading. Of course is easy to write a "patch" which does perform the sequence: set the year back -> start the utility -> restore the date all in the same function. :wink: Probably this program is unfitted for whom deal either with real time stock options and financial stuffs or is an atomic clock maniac...
Never at rest.
Massimo

Avatar
Topic author
Fenerit M
Posts: 1880
Joined: 26.03.2007
Age: 17
With us: 17 years 11 months
Location: Thyrrenian sea

Re: Earth's geochronology map: a work in progress...

Post #12by Fenerit » 25.08.2010, 09:22

Updated the link for the Zonums shape2text, being no longer on the main page. Now is here: http://www.zonums.com/shp2text.html
Never at rest.
Massimo

Avatar
Topic author
Fenerit M
Posts: 1880
Joined: 26.03.2007
Age: 17
With us: 17 years 11 months
Location: Thyrrenian sea

Re: Earth's geochronology map: a work in progress...

Post #13by Fenerit » 26.08.2010, 14:28

As complement, the Perl script for converting in SVG the same TXT file got by the Zonum shape2text converter. The script is equivalent (for SVG) to that above.

Code: Select all

#    Author Fenerit <fenerit@interfree.it>
#    Version 2.0, 18.07.2010
#    This script built a 8192x4096 SVG (Scalable Vector Graphics) file
#   from a set of geographic coordinates yielded in TXT form through the Zonums Shape2Text converter.
#   The resulting file can be edited within graphic programs which support it.
#   Acknowledgements: Dr. Fridger Schrempp <www.celestialmatters.org> for its basic routines.

#   FOR EXTERNAL USAGE:
#    1) Rename or save the TXT file yielded by the Zonums progam as "shapefile.txt"
#   2) Script and dateset's text files must be in the same folder;
#    3) doubleclick on its icon or type within console:
#   "shptxt2svg.pl" < "shapefile.txt" (in Windows cmd: perl "shptxt2svg.pl" < "shapefile.txt" then ENTER)
#   A new file called "shapefile.svg" will be created in the same folder;

   print STDOUT "\n -- Text shapefile to SVG converter. --\n\n -- Start polygons processing... Please wait...\n";
   
##################################  SVG SIZE  #######################################
   
   my $constant = 32/45; # constant;
      
#   -------------------   C O N V E R S I O N    T A B L E:   -----------------------
#
#   0.25k  (256x128)  = 360 *   1 * (32/45) x 180 *   1 * (32/45) --> multiplier   1;
#    0.5k   (512x256)  = 360 *   2 * (32/45) x 180 *   2 * (32/45) --> multiplier   2;
#    1k    (1024x512)  = 360 *   4 * (32/45) x 180 *   4 * (32/45) --> multiplier   4;
#   2k   (2048x1024)  = 360 *   8 * (32/45) x 180 *   8 * (32/45) --> multiplier   8;
#   4k   (4096x2048)  = 360 *  16 * (32/45) x 180 *  16 * (32/45) --> multiplier  16;
#   8k   (8192x4096)  = 360 *  32 * (32/45) x 180 *  32 * (32/45) --> multiplier  32;
#  16k  (16384x8192)  = 360 *  64 * (32/45) x 180 *  64 * (32/45) --> multiplier  64;
#  32k (32768x16384)  = 360 * 128 * (32/45) x 180 * 128 * (32/45) --> multiplier 128;
#  ... thermoionic tubes were best!
   
   my $width = 8192; # set the SVG width;
   my $height = 4096; # set the SVG height;
   my $multiplier = 32; # scalar multiplier (see above);
   
#####################################################################################
   
   my $pixels = "M 0,0 L 0,$height L $width,$height L $width,0 z ";
   my $style = "fill:none;fill-rule:evenodd;stroke:#000000;stroke-width:0.1px;stroke-linecap:butt;stroke-linejoin:miter;stroke-opacity:1";
#   diminuish the stroke-width whether you choose smaller SVG sizes;

   open(TXT, " < shapefile.txt") || die "Can not read the file!\n";
   open(TMP, " > shptxtconv.tmp") || die "Can not write the file!\n";

   {           

   @G_polygon = <TXT>;
   
   }
   
   foreach $_(@G_polygon)

   {

      next if (/^$/);   
      next if (/^\s*shape2Text\s*$/);
      next if (/^\s*Powered by Zonum Solutions\s*$/);
      next if (/^\s*www.zonums.com\s*$/);
#      Change below "Polygon" to "Polyline" with the SHP type (default: Polygon)
      next if (/^\s*Entity type:Polygon\s*$/);
      next if (/^\s*XCoo YCoo\s*$/);   
      
      chomp $_;
   
      ($long,$lat) = split(" ",$_);
         
      if($lat == "Nodes") {
            
         print TMP "\n\n   <path\n";
         print TMP "      style=\"$style\"\n";
         print TMP "      d=\"M ";
               
      }
      
      else
             
      {
      
         $lat[$l] = (abs(90 - $lat)) * $constant * $multiplier;
         $long[$l] = (180 + $long) * $constant * $multiplier;
      
         printf TMP "%8.9f,%8.9f L ",$long[$index],$lat[$index];
            
      }
      
   }
   
close TXT;
close TMP;   

   sleep 2;
   
   print STDOUT "\n\n -- 25% done...\n";
   
   sleep 2;

   open(TMP, " +< shptxtconv.tmp") || die "Can not read and write the file!";

{           

   local $/ = '';
   @S_polygon = <TMP>;
   
   print TMP "\nend__svg__shape\n";

}
   
   foreach $_(@S_polygon)

   {
   
      next if (/^$/);   

      chomp $_;
   
      ($path) = split(" ",$_);
   
      if($path == "<path") {
      
         print TMP "$_";
         print TMP "      \"/>\n";

      }
   
   }

close TMP;

   sleep 2;
   
   print STDOUT "\n\n -- 50% done...\n";
   
   sleep 2;

   open(TMP, " < shptxtconv.tmp") || die "Can not read the file!";
   open(SVGTMP, " > shapefile.tmp") || die "Can not write the file!";

   @T_SVG = <TMP>;

    while ($_ = pop @T_SVG) {
       
      last if (/^end__svg__shape$/);
      print SVGTMP $_;
   
   }
           
close TMP;
close SVGTMP;

   sleep 2;
   
   print STDOUT "\n\n -- 100% done! Parsing SVG...\n";
   
   sleep 2;

   open(SVGTMP, " < shapefile.tmp") || die "Can not read the file!";
   open(SVG, " > shapefile.svg") || die "Can not write the file!";   
   
   print SVG "<svg\n";
   print SVG "   width=\"$width\"\n";
   print SVG "   height=\"$height\">\n";
   print SVG "   <g>\n";
   print SVG "   <path\n";
   print SVG "      style=\"$style\"\n";
   print SVG "      d=\"$pixels\"\n";
   print SVG "      />\n";
   
   @SVG = <SVGTMP>;

    while ($_ = pop @SVG) {
       
      print SVG $_;
   
   }
   
   print SVG "   </g>\n";
   print SVG "</svg>\n";
   
close SVGTMP;
close SVG;

   sleep 2;
   
   unlink ("shptxtconv.tmp");
   unlink ("shapefile.tmp");

   print STDOUT "\n -- shapefile.svg ready! --\n"


FYI, when the big shapes like the Canada shape or the Arabia shape are imported in Inkscape 0.47, I can handle it unless I've RAM, for example to zoom in it, then Inkscape crash. Same business for what concern GIMP: with the trace paths active it is slowest, as raster it work fine. Now, when I export in SVG the same shapes with Qgis, Inkscape neither open them. Same story GIMP. Thus, the crashing of Inkscape in the first case, probably is due to the bad managing of the memory threshold. Now is ready the newest 0.48. I will test soon. Below the U.S shape how look after conversion at 1024x512 px:

Image

Image
Never at rest.
Massimo

Avatar
Topic author
Fenerit M
Posts: 1880
Joined: 26.03.2007
Age: 17
With us: 17 years 11 months
Location: Thyrrenian sea

Re: Earth's geochronology map: a work in progress...

Post #14by Fenerit » 21.09.2010, 01:12

News from the map (level 5):

To Do: PDF/image for legend. Then a question: whether to add to the map also the age (geologic scale) of the ocean floor or not. With such solution my perplexity concern the "Christmas ball" effect of the whole...
Never at rest.
Massimo

Avatar
Topic author
Fenerit M
Posts: 1880
Joined: 26.03.2007
Age: 17
With us: 17 years 11 months
Location: Thyrrenian sea

Re: Earth's geochronology map: a work in progress...

Post #15by Fenerit » 23.12.2010, 12:54

...Hi! Vectorized Mexico and inserted the shape within the Caribbean shape...

Image
Never at rest.
Massimo

Avatar
Topic author
Fenerit M
Posts: 1880
Joined: 26.03.2007
Age: 17
With us: 17 years 11 months
Location: Thyrrenian sea

Re: Earth's geochronology map: a work in progress...

Post #16by Fenerit » 23.01.2011, 10:05

...Another issue.

Image

Seem that the Large Igneous Provinces (LIPs) as seen in the 3D geologic add-on doesn't matches with the areas ascribed to them by the geologic map; just some features at poles match. Such behaviour concern only the continental LIPs, while the oceanic ones matches perfectly with the seafloor reliefs as seen in the bathymetry map. Being in the same dataset, I've check their projections; but all is correct. Once finished the map, probably I will update the 3D LIPs by extracting the zones from such shapes; thus consider the existent LIPs model as "generic".
Never at rest.
Massimo

MiR
Posts: 247
Joined: 01.02.2010
With us: 15 years 1 month
Location: Germany

Re: Earth's geochronology map: a work in progress...

Post #17by MiR » 24.01.2011, 17:42

again: this looks very neat and accurate!

Good job, Massimo :)

Michael

Avatar
Topic author
Fenerit M
Posts: 1880
Joined: 26.03.2007
Age: 17
With us: 17 years 11 months
Location: Thyrrenian sea

Re: Earth's geochronology map: a work in progress...

Post #18by Fenerit » 24.01.2011, 20:25

Hi Michael, thank-you for supporting this add-on's enterprise and for having interrupt my monologue in these matters :wink: ; I do not have too much spare time now, but I'm facing with the most tedious work: to write the ages, which I do manually (with CMODtracer) mostly because I've under control the featured zones so that I can always check it against errors, and secondary because it's more fast this way than the required operations to get them through the centroids (USGS datesets either doesn't have exchange files or the text is polygonal). Bear in mind that there will be also the 3D faults for several world's areas.
Never at rest.
Massimo

MiR
Posts: 247
Joined: 01.02.2010
With us: 15 years 1 month
Location: Germany

Re: Earth's geochronology map: a work in progress...

Post #19by MiR » 25.01.2011, 17:30

Looks like an enormous work!
I'm looking forward to it. :)

Michael

Avatar
Topic author
Fenerit M
Posts: 1880
Joined: 26.03.2007
Age: 17
With us: 17 years 11 months
Location: Thyrrenian sea

Re: Earth's geochronology map: a work in progress...

Post #20by Fenerit » 25.01.2011, 17:47

MiR wrote:Looks like an enormous work!
I'm looking forward to it. :)

Michael

Lots of coffee around (usually I play with it in the night) and a personalized Celestia compilation without the right menu, otherwise during the right togglin the button's release often does show it causing interferences, then is fairly easy and fast (I collect points age by age and then do process them with Perl).
Never at rest.
Massimo


Return to “Textures”