Hi Johannes, On 2019-10-23 at 10:49 -07, Johannes Radinger <[email protected]> wrote... > So maybe another way (if rotated pole is not possible with GRASS/gdal)
Rotated pole works fine in GRASS. Here are my notes from when I did this a few years ago. -k. https://lists.osgeo.org/pipermail/grass-user/2011-October/062180.html #+BEGIN_SRC bash :results verbatim # Normal and rotated GRASS locations trash Gnorm Grot grass74 -c EPSG:4326 Gnorm # Import something into the normal location v.import input=~/data/Zwally_2012/sectors output=Z cat << EOF > ./Gnorm/PERMANENT/PROJ_INFO name: General Oblique Transformation datum: wgs84 towgs84: 0.000,0.000,0.000 proj: ob_tran o_proj: latlon ellps: wgs84 a: 6378137.0000000000 es: 0.0066943800 f: 298.2572235630 lat_0: 0.0000000000 lon_0: 180.0000000000 o_lat_p: 18.0 o_lon_p: -200.0 EOF # rotated_pole:grid_north_pole_latitude = 18. ; # rotated_pole:grid_north_pole_longitude = -200. cat << EOF > ./Gnorm/PERMANENT/PROJ_UNITS unit: degree units: degrees meters: .0174532925 EOF grass74 -e -c EPSG:4326 Grot grass74 ./Grot/PERMANENT v.proj location=Gnorm input=Z g.region vector=Z d.mon start=wx0 d.erase d.vect Z d.grid 1:0:0 #+END_SRC _______________________________________________ grass-user mailing list [email protected] https://lists.osgeo.org/mailman/listinfo/grass-user
