{"id":240,"date":"2016-01-15T05:32:11","date_gmt":"2016-01-15T05:32:11","guid":{"rendered":"https:\/\/existencia.org\/pro\/?p=240"},"modified":"2016-01-15T05:35:35","modified_gmt":"2016-01-15T05:35:35","slug":"tropict-a-clearer-depiction-of-the-tropics","status":"publish","type":"post","link":"https:\/\/existencia.org\/pro\/?p=240","title":{"rendered":"Tropict: A clearer depiction of the tropics"},"content":{"rendered":"

Tropict is a set of python and R scripts that adjust the globe to make land masses in the tropics fill up more visual real estate. It does this by exploiting the ways continents naturally “fit into” each other, splicing out wide areas of empty ocean and nestling the continents closer together.<\/p>\n

All Tropict scripts are designed to show the region between 30°S and 30°N. In an equirectangular projection, that looks like this:<\/p>\n

\"original\"<\/a><\/p>\n

It is almost impossible to see what is happening on land: the oceans dominate. By removing open ocean and applying the Gall-Peters projection, we get a clearer picture:<\/p>\n

\"version4\"<\/a><\/p>\n

There’s even a nice spot for a legend in the lower-left! Whether for convenience or lack of time, the tools I’ve made to allow you to make these maps are divided between R and Python. Here’s a handy guide for which tool to use:<\/p>\n

\"decisions\"<\/a>
\n
\n(1)<\/i> Supported image formats are listed in the Pillow<\/a> documentation.
\n(2)<\/i> A TSR file is a Tropict Shapefile Reinterpretation file, and includes the longitudinal shifts for each hemisphere.
\n<\/font><\/p>\n

Let’s say you find yourself with a NetCDF file in need of Tropiction, called bio-2.nc4<\/tt>. It’s already clipped to between 30°S and 30°N. The first step is to call splice_grid.py<\/tt> to create a Tropicted NetCDF:<\/p>\n

\r\npython ..\/splice_grid.py subjects\/bio-2.nc4 ..\/bio-2b.nc4\r\n<\/pre>\n
But that NetCDF doesn’t show country boundaries.  To show country boundaries, you can follow the example<\/a> for using draw_map.R<\/tt>:<\/p>\n
\r\nlibrary(ncdf4)\r\nlibrary(RColorBrewer)\r\n\r\n## Open the Tropicted NetCDF\r\ndatabase <- nc_open("bio-2b.nc4")\r\n## Extract one variable\r\nmap <- ncvar_get(database, "change")\r\n\r\n## Identify the range of values there\r\nmaxmap <- max(abs(map), na.rm=T)\r\n\r\n## Set up colors centered on 0\r\ncolors <- rev(brewer.pal(11,"RdYlBu"))\r\nbreaks <- seq(-maxmap, maxmap, length.out=12)\r\n\r\n## Draw the NetCDF image as a background\r\nsplicerImage(map, colors, breaks=breaks)\r\n## Add country boundaries\r\naddMap(border="#00000060")\r\n## Add seams where Tropict knits the map together\r\naddSeams(col="#00000040")\r\n<\/pre>\n
Here’s an example of the final result, for a bit of my coffee work:<\/p>\n
\"arabica-futureb\"<\/a><\/p>\n
For more details, check out the documentation at the GitHub page<\/a>!<\/p>\n
And just for fun, here were two previous attempts of re-hashing the globe:<\/p>\n
\"version1\"<\/a><\/p>\n
I admit that moving Australia and Hawaii into the India Ocean was over-zealous, but they fill up the space so well!<\/p>\n
\"version3\"<\/a><\/p>\n
Here I can still use the slick division between Indonesian and Papua New Guinea and Hawaii fits right on the edge, but Australia gets split in two.<\/p>\n
Enjoy the tropics!<\/p>\n","protected":false},"excerpt":{"rendered":"
Tropict is a set of python and R scripts that adjust the globe to make land masses in the tropics fill up more visual real estate. It does this by exploiting the ways continents naturally “fit into” each other, splicing out wide areas of empty ocean and nestling the continents closer together. All Tropict scripts […]<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"jetpack_post_was_ever_published":false,"_jetpack_newsletter_access":"","_jetpack_dont_email_post_to_subs":false,"_jetpack_newsletter_tier_id":0,"_jetpack_memberships_contains_paywalled_content":false,"_jetpack_memberships_contains_paid_content":false,"footnotes":"","jetpack_publicize_message":"","jetpack_publicize_feature_enabled":true,"jetpack_social_post_already_shared":true,"jetpack_social_options":{"image_generator_settings":{"template":"highway","default_image_id":0,"font":"","enabled":false},"version":2}},"categories":[4],"tags":[],"class_list":["post-240","post","type-post","status-publish","format-standard","hentry","category-software"],"jetpack_publicize_connections":[],"jetpack_featured_media_url":"","jetpack_sharing_enabled":true,"jetpack_shortlink":"https:\/\/wp.me\/p4Zh9E-3S","_links":{"self":[{"href":"https:\/\/existencia.org\/pro\/index.php?rest_route=\/wp\/v2\/posts\/240","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/existencia.org\/pro\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/existencia.org\/pro\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/existencia.org\/pro\/index.php?rest_route=\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/existencia.org\/pro\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=240"}],"version-history":[{"count":0,"href":"https:\/\/existencia.org\/pro\/index.php?rest_route=\/wp\/v2\/posts\/240\/revisions"}],"wp:attachment":[{"href":"https:\/\/existencia.org\/pro\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=240"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/existencia.org\/pro\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=240"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/existencia.org\/pro\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=240"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}