[GRASS-user] ITRF92 Not found.

Markus Neteler neteler at osgeo.org
Thu Mar 18 08:16:51 PDT 2021


On Tue, Mar 16, 2021 at 9:51 AM Markus Neteler <neteler at osgeo.org> wrote:
>
> Hi,
>
> On Sun, Mar 14, 2021 at 3:32 PM javier Garcia Prieto
> <fjgarciaprieto at hotmail.com> wrote:
> >
> > Thank you!
> > You are really great! This is going to be so much useful to Mexican people, You are very kind.
>
> Glad to help. Now it is merged, i.e. it is already available in the
> current source code of GRASS GIS 7.8.git and will be part of GRASS GIS
> 7.8.6 and later.

BTW, I made a test which you indicated to me off-list, with the
udpated GRASS GIS 7.8.git version:

gdalinfo republica2017_r15m.bil
Driver: EHdr/ESRI .hdr Labelled
Files: republica2017_r15m.bil
       republica2017_r15m.bil.ovr
       republica2017_r15m.bil.aux.xml
       republica2017_r15m.hdr
       republica2017_r15m.stx
       republica2017_r15m.prj
Size is 218763, 130943
Coordinate System is:
GEOGCRS["ITRF92",
    DATUM["International Terrestrial Reference Frame 1992",
        ELLIPSOID["GRS 1980",6378137,298.257222101,
            LENGTHUNIT["metre",1]],
        ID["EPSG",6651]],
    PRIMEM["Greenwich",0,
        ANGLEUNIT["Degree",0.0174532925199433]],
    CS[ellipsoidal,2],
        AXIS["longitude",east,
            ORDER[1],
            ANGLEUNIT["Degree",0.0174532925199433]],
        AXIS["latitude",north,
            ORDER[2],
            ANGLEUNIT["Degree",0.0174532925199433]]]
Data axis to CRS axis mapping: 1,2
Origin = (-117.124268882302502,32.718626154577500)
Pixel Size = (0.000138888889000,-0.000138888889000)
Corner Coordinates:
Upper Left  (-117.1242689,  32.7186262) (117d 7'27.37"W, 32d43' 7.05"N)
Lower Left  (-117.1242689,  14.5320984) (117d 7'27.37"W, 14d31'55.55"N)
Upper Right ( -86.7405189,  32.7186262) ( 86d44'25.87"W, 32d43' 7.05"N)
Lower Right ( -86.7405189,  14.5320984) ( 86d44'25.87"W, 14d31'55.55"N)
Center      (-101.9323939,  23.6253623) (101d55'56.62"W, 23d37'31.30"N)
Band 1 Block=218763x1 Type=Int16, ColorInterp=Undefined
  Min=-32.000 Max=5610.000
  Minimum=-32.000, Maximum=5610.000, Mean=1041.610, StdDev=820.476
  NoData Value=-32768
  Overviews: 109382x65472, 54691x32736, 27346x16368, 13673x8184,
6837x4092, 3419x2046, 1710x1023, 855x512, 428x256, 214x128
  Metadata:
    SourceBandIndex=0
    STATISTICS_SKIPFACTORX=1
    STATISTICS_SKIPFACTORY=1

# generate location from dataset
grass78 -c republica2017_r15m.bil /mnt/data/grassdata/mexico

# verify projection
GRASS 7.8.6dev (mexico):~ > g.proj -w
GEOGCRS["ITRF92",
    DATUM["International Terrestrial Reference Frame 1992",
        ELLIPSOID["GRS 1980",6378137,298.257222101,
            LENGTHUNIT["metre",1]],
        ID["EPSG",6651]],
    PRIMEM["Greenwich",0,
        ANGLEUNIT["Degree",0.0174532925199433]],
    CS[ellipsoidal,2],
        AXIS["longitude",east,
            ORDER[1],
            ANGLEUNIT["Degree",0.0174532925199433]],
        AXIS["latitude",north,
            ORDER[2],
            ANGLEUNIT["Degree",0.0174532925199433]]]

# import DEM
r.import republica2017_r15m.bil output=republica2017_r15m

# set color table
r.colors republica2017_r15m color=srtm_plus

The result looks good to me.

best,
Markus


More information about the grass-user mailing list