<div dir="ltr">Dear Moritz and Helmut,<div><br></div><div>Thank you for letting me know that the GRASS GUI (and g.proj) have the ability to work directly with proj4 strings.</div><div><br></div><div>  -k.</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Oct 6, 2017 at 10:34 AM, Helmut Kudrnovsky <span dir="ltr"><<a href="mailto:hellik@web.de" target="_blank">hellik@web.de</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Ken Mankoff wrote<br>
<div><div class="h5">> Dear List,<br>
><br>
> Most of the time when I work with data in GRASS it is a provided on a<br>
> standard GIS grid (EPSG code) and I can work with it easily. I create a<br>
> location using the "-c" option when I start GRASS, or I import it into an<br>
> existing mapset with r.import.<br>
><br>
> I'm now working with data and all I know is this:<br>
><br>
> +proj=stere +lat_0=90 +lat_ts=71 +lon_0=-39 +k=1 +x_0=0 +y_0=0<br>
> +datum=WGS84<br>
> +units=m +no_defs”<br>
><br>
> In the past, when working with a few variables and this information, and<br>
> lat/lon grids in the NetCDF file, I would export the data (using Python)<br>
> to<br>
> lat,lon,data, and then r.in.ascii and project it that way. However, now<br>
> I'm<br>
> working with 14,000 high resolution daily records on the above projection.<br>
><br>
> It seems the only way to do this efficiently is to work in the projection<br>
> of the data, and r.import the few other variables to this projection. Or<br>
> perhaps cdo or gdal or some other 3rd party tool can quickly re-project<br>
> everything for me to a known EPSG code. I've solved this before with<br>
> vectors, ogr2ogr, and arbitrary proj4 strings.<br>
><br>
> Can anyone here suggest how to best get this data into GRASS?<br>
><br>
> I've tried using gdalwarp:<br>
><br>
> gdalwarp -of GTIFF -s_srs "+proj=stere +lat_0=90 +lat_ts=71 +lon_0=-39<br>
> +k=1<br>
> +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs" -t_srs EPSG:3413<br>
> NETCDF:"${file}":RU out.tif<br>
><br>
> But am concerned because when I then run gdalinfo on the out.tif file, it<br>
> reports strange bounds such as:<br>
><br>
> Origin = (-909.965548780931840,1546.<wbr>900601999830769)<br>
> Pixel Size = (7.477558376763198,-7.<wbr>477558376763198)<br>
> Corner Coordinates:<br>
> Upper Left  (    -909.966,    1546.901) (165d27'58.34"E, 89d59' 0.36"N)<br>
> Lower Left  (    -909.966,   -1264.661) ( 80d44'10.54"W, 89d59' 8.22"N)<br>
> Upper Right (     794.918,    1546.901) (107d48' 8.56"E, 89d59' 2.20"N)<br>
> Lower Right (     794.918,   -1264.661) ( 12d50'53.17"W, 89d59'10.36"N)<br>
> Center      (     -57.524,     141.120) (157d10'37.34"E, 89d59'54.94"N)<br>
><br>
> Which seems to be very close to the pole (both lower and upper at >89<br>
> deg).<br>
> This should be covering a much large portion of the Arctic.<br>
<br>
</div></div>e.g. you can use a proj4 string in the startup wizard to create a new<br>
location.<br>
<br>
just tested with your proj4 string,<br>
<br>
g.proj -j<br>
<span class="">+proj=stere<br>
+lat_0=90<br>
+lat_ts=71<br>
+lon_0=-39<br>
+k=1<br>
+x_0=0<br>
+y_0=0<br>
</span>+no_defs<br>
+a=6378137<br>
+rf=298.257223563<br>
+towgs84=0.000,0.000,0.000<br>
+to_meter=1<br>
<br>
HTH<br>
<br>
<br>
<br>
-----<br>
best regards<br>
Helmut<br>
--<br>
Sent from: <a href="http://osgeo-org.1560.x6.nabble.com/Grass-Users-f3884509.html" rel="noreferrer" target="_blank">http://osgeo-org.1560.x6.<wbr>nabble.com/Grass-Users-<wbr>f3884509.html</a><br>
______________________________<wbr>_________________<br>
grass-user mailing list<br>
<a href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/grass-user" rel="noreferrer" target="_blank">https://lists.osgeo.org/<wbr>mailman/listinfo/grass-user</a></blockquote></div><br></div>