"Fossies" - the Fresh Open Source Software Archive  

Source code changes of the file "Demos/Cartography_demo.pm" between
PDL-2.076.tar.gz and PDL-2.077.tar.gz

About: PDL (Perl Data Language) aims to turn perl into an efficient numerical language for scientific computing (similar to IDL and MatLab).

Cartography_demo.pm  (PDL-2.076)	:	Cartography_demo.pm  (PDL-2.077)
#
package PDL::Demos::Cartography_demo;		package PDL::Demos::Cartography_demo;
use PDL;		sub init {'
use PDL::Graphics::PGPLOT::Window;		use PDL::Graphics::PGPLOT::Window;
use PDL::Transform::Cartography;		use PDL::Transform::Cartography;
		'}
		sub done {'
		$w->close;
		undef $w;
		'}
use File::Spec;		sub info {('cartography','Cartographic projections (Req.: PGPLOT)')}
PDL::Demos::Routines->import();
sub comment($);
sub act($);
sub output;
sub run {
local($PDL::debug) = 0;
local($PDL::verbose) = 0;
##$ENV{PGPLOT_XW_WIDTH}=0.6;		##$ENV{PGPLOT_XW_WIDTH}=0.6;
$ENV{PGPLOT_DEV} = $^O =~ /MSWin32/ ? '/GW' :		$ENV{PGPLOT_DEV} = $^O =~ /MSWin32/ ? '/GW' :
defined($ENV{PGPLOT_DEV}) ? $ENV{PGPLOT_DEV} : "/XWIN";		defined($ENV{PGPLOT_DEV}) ? $ENV{PGPLOT_DEV} : "/XWIN";
unless( PDL->rpiccan('JPEG') ) {		my @demos = (
comment q|		PDL->rpiccan('JPEG') ? () :
		[comment => q|
This demo illustrates the PDL::Transform::Cartography module.		This demo illustrates the PDL::Transform::Cartography module.
It requires PGPLOT and also the ability to read/write JPEG images.		It requires PGPLOT and also the ability to read/write JPEG images.
You don't seem to have that ability at the moment -- this is likely		You don't seem to have that ability at the moment -- this is likely
because you do not have NetPBM installed. See the man page for PDL::IO::Pic.		because you do not have NetPBM installed. See the man page for PDL::IO::Pic.
I'll continue with the demo anyway, but it will likely crash on the		I'll continue with the demo anyway, but it will likely crash on the
earth_image('day') call on the next screen.		earth_image('day') call on the next screen.
|		|],
}
comment q|		[comment => q|
This demo illustrates the PDL::Transform::Cartography module.		This demo illustrates the PDL::Transform::Cartography module.
It also requires PGPLOT support: you must have PGPLOT installed to run it.		It also requires PGPLOT support: you must have PGPLOT installed to run it.
PDL::Transform::Cartography includes a global earth vector coastline map		PDL::Transform::Cartography includes a global earth vector coastline map
and night and day world image maps, as well as the infrastructure for		and night and day world image maps, as well as the infrastructure for
transforming them to different coordinate systems.		transforming them to different coordinate systems.
If you are using an RGB-enabled version of PGPLOT, then the map		If you are using an RGB-enabled version of PGPLOT, then the map
images will appear in color on your screen. Otherwise they will		images will appear in color on your screen. Otherwise they will
appear in greyscale. (The vector data will appear in color		appear in greyscale. (The vector data will appear in color
regardless).		regardless).
|;		|],
act q|		[act => q|
### Load the necessary modules		### Load the necessary modules
use PDL::Graphics::PGPLOT::Window;		use PDL::Graphics::PGPLOT::Window;
use PDL::Transform::Cartography;		use PDL::Transform::Cartography;
### Get the vector coastline map (and a lon/lat grid), and load the Earth		### Get the vector coastline map (and a lon/lat grid), and load the Earth
### RGB daytime image -- both of these are built-in to the module. The		### RGB daytime image -- both of these are built-in to the module. The
### coastline map is a set of (X,Y,Pen) vectors.		### coastline map is a set of (X,Y,Pen) vectors.
$coast = earth_coast() -> glue( 1, scalar graticule(15,1) );		$coast = earth_coast() -> glue( 1, scalar graticule(15,1) );
print "Coastline data are a collection of vectors: ",		print "Coastline data are a collection of vectors: ",
join("x",$coast->dims),"\n";		join("x",$coast->dims),"\n";
$map = earth_image('day');		$map = earth_image('day');
print "Map data are RGB: ",join("x",$map->dims),"\n\n";		print "Map data are RGB: ",join("x",$map->dims),"\n\n";
|;		|],
act q&		[act => q&
### Map data are stored natively in Plate Caree format.		### Map data are stored natively in Plate Caree format.
### The image contains a FITS header that contains coordinate system info.		### The image contains a FITS header that contains coordinate system info.
print "FITS HEADER INFORMATION:\n";		print "FITS HEADER INFORMATION:\n";
for $_(sort keys %{$map->hdr}){		for $_(sort keys %{$map->hdr}){
next if(m/SIMPLE/ || m/HISTORY/ || m/COMMENT/);		next if(m/SIMPLE/ || m/HISTORY/ || m/COMMENT/);
printf (" %8s: %10s%s", $_, $map->hdr->{$_}, (++$i%3) ? " " : "\n");		printf (" %8s: %10s%s", $_, $map->hdr->{$_}, (++$i%3) ? " " : "\n");
}		}
print "\n";		print "\n";
$dev = $^O =~ /MSWin32/ ? '/GW' :		$dev = $^O =~ /MSWin32/ ? '/GW' :
defined($ENV{PGPLOT_DEV}) ? $ENV{PGPLOT_DEV} : "/XW";		defined($ENV{PGPLOT_DEV}) ? $ENV{PGPLOT_DEV} : "/XW";
# $dev = $^O =~ /MSWin/i ? '/GW' : '/xw';		# $dev = $^O =~ /MSWin/i ? '/GW' : '/xw';
$w = pgwin(Dev=> $dev, size=>[8,6]);		$w = pgwin(Dev=> $dev, size=>[8,6]);
$w->fits_imag($map, {Title=>"NASA/MODIS Earth Map (Plate Caree)",J=>0});		$w->fits_imag($map, {Title=>"NASA/MODIS Earth Map (Plate Caree)",J=>0});
&;		&],
act q&		[act => q&
### The map data are co-aligned with the vector data, which can be drawn		### The map data are co-aligned with the vector data, which can be drawn
### on top of the window with the "->lines" PGPLOT method. The		### on top of the window with the "->lines" PGPLOT method. The
### clean_lines method breaks lines that pass over the map's singularity		### clean_lines method breaks lines that pass over the map's singularity
### at the 180th parallel.		### at the 180th parallel.
$w->hold;		$w->hold;
$w->lines( $coast -> clean_lines );		$w->lines( $coast -> clean_lines );
$w->release;		$w->release;
&;		&],
act q&		[act => q&
### There are a large number of map projections -- to list them all,		### There are a large number of map projections -- to list them all,
### say "??cartography" in the perldl or pdl2 shell. Here are four		### say "??cartography" in the perldl or pdl2 shell. Here are four
### of them:		### of them:
$w->close; # Close old window		$w->close; # Close old window
undef $w;		undef $w;
#$dev = $^O =~ /MSWin/i ? '/GW' : '/xw';		#$dev = $^O =~ /MSWin/i ? '/GW' : '/xw';
$dev = $^O =~ /MSWin32/ ? '/GW' :		$dev = $^O =~ /MSWin32/ ? '/GW' :
defined($ENV{PGPLOT_DEV}) ? $ENV{PGPLOT_DEV} : "/XW";		defined($ENV{PGPLOT_DEV}) ? $ENV{PGPLOT_DEV} : "/XW";
skipping to change at line 123		skipping to change at line 118
($tx, $t, $px, $opt ) = @_;		($tx, $t, $px, $opt ) = @_;
$w->fits_imag( $map->map( $tx, $px, $opt ),{Title=>$t, CharSize=>1.5, J=>1} );		$w->fits_imag( $map->map( $tx, $px, $opt ),{Title=>$t, CharSize=>1.5, J=>1} );
$w->hold; $w->lines( $coast -> apply( $tx ) -> clean_lines ); $w->release;		$w->hold; $w->lines( $coast -> apply( $tx ) -> clean_lines ); $w->release;
}		}
## (The "or" option specifies the output range of the mapping)		## (The "or" option specifies the output range of the mapping)
draw( t_mercator, "Mercator Projection", [400,300] );		draw( t_mercator, "Mercator Projection", [400,300] );
draw( t_aitoff, "Aitoff / Hammer", [400,300] );		draw( t_aitoff, "Aitoff / Hammer", [400,300] );
draw( t_gnomonic, "Gnomonic", [400,300],{or=>[[-3,3],[-2,2]]} );		draw( t_gnomonic, "Gnomonic", [400,300],{or=>[[-3,3],[-2,2]]} );
draw( t_lambert, "Lambert Conformal Conic",[400,300],{or=>[[-3,3],[-2,2]]} );		draw( t_lambert, "Lambert Conformal Conic",[400,300],{or=>[[-3,3],[-2,2]]} );
		&],
&;		[act => q|
act q|
### You can create oblique projections by feeding in a different origin.		### You can create oblique projections by feeding in a different origin.
### Here, the origin is centered over North America.		### Here, the origin is centered over North America.
draw( t_mercator( o=>[-90,40] ), "Mercator Projection", [400,300] );		draw( t_mercator( o=>[-90,40] ), "Mercator Projection", [400,300] );
draw( t_aitoff ( o=>[-90,40] ), "Aitoff / Hammer", [400,300] );		draw( t_aitoff ( o=>[-90,40] ), "Aitoff / Hammer", [400,300] );
draw( t_gnomonic( o=>[-90,40] ), "Gnomonic",[400,300],{or=>[[-3,3],[-2,2]]} );		draw( t_gnomonic( o=>[-90,40] ), "Gnomonic",[400,300],{or=>[[-3,3],[-2,2]]} );
draw( t_lambert( o=>[-90,40] ), "Lambert ",[400,300],{or=>[[-3,3],[-2,2]]} );		draw( t_lambert( o=>[-90,40] ), "Lambert ",[400,300],{or=>[[-3,3],[-2,2]]} );
|;		|],
act q|		[act => q|
### There are three main perspective projections (in addition to special		### There are three main perspective projections (in addition to special
### cases like stereographic and gnomonic projection): orthographic,		### cases like stereographic and gnomonic projection): orthographic,
### vertical, and true perspective. The true perspective has options for		### vertical, and true perspective. The true perspective has options for
### both downward-looking and aerial-view coordinate systems.		### both downward-looking and aerial-view coordinate systems.
draw( t_orthographic( o=>[-90,40] ),		draw( t_orthographic( o=>[-90,40] ),
"Orthographic", [400,300]);		"Orthographic", [400,300]);
draw( t_vertical( r0=> (2 + 1), o=>[-90,40] ),		draw( t_vertical( r0=> (2 + 1), o=>[-90,40] ),
"Vertical (Altitude = 2 r\\\\de\\\\u)", [400,300]);		"Vertical (Altitude = 2 r\\\\de\\\\u)", [400,300]);
draw( t_perspective( r0=> (2 + 1), o=>[-90,40] ),		draw( t_perspective( r0=> (2 + 1), o=>[-90,40] ),
"True Perspective (Altitude= 2 r\\\\de\\\\u)", [400,300]);		"True Perspective (Altitude= 2 r\\\\de\\\\u)", [400,300]);
# Observer is 0.1 earth-radii above surface, lon 117W, lat 31N (over Tijuana).		# Observer is 0.1 earth-radii above surface, lon 117W, lat 31N (over Tijuana).
# view is 45 degrees below horizontal, azimuth -22 (338) degrees.		# view is 45 degrees below horizontal, azimuth -22 (338) degrees.
draw( t_perspective( r0=> 1.1, o=>[-117,31], cam=>[-22,-45,0] ),		draw( t_perspective( r0=> 1.1, o=>[-117,31], cam=>[-22,-45,0] ),
"Aerial view of West Coast of USA", [400,300],		"Aerial view of West Coast of USA", [400,300],
{or=>[[-60,60],[-45,45]], method=>'linear'});		{or=>[[-60,60],[-45,45]], method=>'linear'});
|;		|],
comment q|		[comment => q|
That concludes the basic cartography demo. Numerous other transforms		That concludes the basic cartography demo. Numerous other transforms
are available.		are available.
Because PDL's cartographic transforms work with in the Transform module		Because PDL's cartographic transforms work with in the Transform module
and are invertible, it's easy to use them both forwards and backwards.		and are invertible, it's easy to use them both forwards and backwards.
In particular, the perspective transformation is useful for ingesting		In particular, the perspective transformation is useful for ingesting
scientific image data of the Earth or other planets, and converting to		scientific image data of the Earth or other planets, and converting to
a map of the imaged body.		a map of the imaged body.
Similarly, scanned images of map data can easily be converted into		Similarly, scanned images of map data can easily be converted into
lat/lon coordinates or reprojected to make other projections.		lat/lon coordinates or reprojected to make other projections.
|;		|],
		);
$w->close;		sub demo { @demos }
undef $w;
}
1;		1;
End of changes. 23 change blocks.
37 lines changed or deleted		29 lines changed or added

---