Re: [GRASS-user] Multidirectional hillshading scripts / addons?

[Carlos Henrique Grohmann de Carvalho via grass-user]

WhiteBox is another option (there is a QGIS plugin)

https://www.whiteboxgeo.com/manual/wbt_book/available_tools/geomorphometric_analysis.html#MultidirectionalHillshade

Carlos



On Sat, Feb 17, 2024 at 10:44 AM Paulo van Breugel via grass-user <
grass-user@lists.osgeo.org> wrote:

> Hi,
>
> Not sure about grass gis, but have you tried the multidirectional option
> when styling a raster as a hillshade in QGIS? See
> https://bnhr.xyz/2019/01/22/hillshade-in-qgis.html.
>
> Cheers,
>
> Paulo
>
>
>
>
>
> On Fri, Feb 16, 2024 at 11:20 PM B H via grass-user <
> grass-user@lists.osgeo.org> wrote:
>
>> Hi,
>> I am trying to create a hillshading effect using DEM rasters similar to
>> ESRI multidirectional hillshading, but I am unable to get anything  in.
>> grass to perform as good as ESRI (or even gdaldem multidirectional shading).
>> I am hoping someone in this group has some scripts /addon/tips that I can
>> try.
>>
>> So far I have tried the following on a few DEMs ( from 3m  to 100 m grid
>> sizes).
>> 1) gdaldem  multidirectional hillshading : It's not grass based, but so
>> far has performed better than other grass options :  with the parameter
>> zfactor between 1 and 10. But  it does not seem to be quite as good as ESRI.
>> 2) Following code is based on some old discussion on this email list.(Its
>> seems to exaggerate the hills and valleys more than I would expect)
>> https://gist.github.com/amcox/b54d65f08fe4448840462901fd6350d4
>> 3) Other variants of hillshading like r.skyview : Tried with several
>> parameters but it seems to be inferior to gdaldem (It seems to be a 'washed
>> out' version of the multidirectional. hillshading).
>>
>> Thanks,
>> HB
>> ___
>> grass-user mailing list
>> grass-user@lists.osgeo.org
>> https://lists.osgeo.org/mailman/listinfo/grass-user
>>
> ___
> grass-user mailing list
> grass-user@lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/grass-user
>


-- 
Prof. Carlos Henrique Grohmann
Institute of Energy and Environment - Univ. of São Paulo, Brazil
Editor-in-Chief - Brazilian Journal of Geology

http://carlosgrohmann.com
https://spamlab.github.io/
http://orcid.org/-0001-5073-5572

Can’t stop the signal.
___
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] single window gui in main branch

[Carlos Henrique Grohmann de Carvalho]

Thanks Martin. Works fine

C

On Mon, Jan 17, 2022 at 6:00 PM Martin Landa  wrote:

> Hi,
>
> pá 14. 1. 2022 v 14:02 odesílatel Carlos Henrique Grohmann de Carvalho
>  napsal:
> > And while it worked fine on linux, I couldn't get it to work on Windows.
> Is it not available yet? I tried the last executable from
> https://wingrass.fsv.cvut.cz/grass80/x86_64/
>
> try
>
> https://wingrass.fsv.cvut.cz/grass81/
>
> Ma
>
> --
> Martin Landa
> http://geo.fsv.cvut.cz/gwiki/Landa
> http://gismentors.cz/mentors/landa
>


-- 
Prof. Carlos Henrique Grohmann
Institute of Energy and Environment - Univ. of São Paulo, Brazil
Editor-in-Chief - Brazilian Journal of Geology

http://carlosgrohmann.com
https://spamlab.github.io/
http://orcid.org/-0001-5073-5572

Can’t stop the signal.
___
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] single window gui in main branch

[Carlos Henrique Grohmann de Carvalho]

Hello all

Just to add my experience. I compiled releasebranch_8 in Linux (Pop!OS
21.10). The single window is great! I have been using a lot o GRASS+Jupyer
lately (running things from Jupyter but with GRASS open for visualization)
and the ability to do just one alt-tab to GRASS is really nice (instead of
choosing the window to switch).

I did find it to be a bit slow to start, and also somewhat unresponsive
when resizing panels.

And while it worked fine on linux, I couldn't get it to work on Windows. Is
it not available yet? I tried the last executable from
https://wingrass.fsv.cvut.cz/grass80/x86_64/

Carlos



On Thu, Jan 13, 2022 at 4:14 PM Micha Silver  wrote:

> Wow, impressive!
>
> I especially like the search option in the right panel toolbox. Also the
> ability to switch the python and command shell to full screen is great.
>
>
> Best regards, and thanks for the nice work!
>
> Micha
>
>
> On 11/01/2022 20:16, Veronica Andreo wrote:
> > Hello everyone,
> >
> > Now that single window is merged into main branch I'd like to test :)
> >
> > Where's the switch? I just updated and compiled main branch and went
> > through all tabs and options in the settings?preferences menu from the
> > GUI, but I could not find anything saying single window something.
> > What am I doing wrong here?
> >
> > Can this be related to
> >
> https://github.com/OSGeo/grass/blob/f57940e8a096c593f23b2ce3f9653547f24b8e92/gui/wxpython/core/settings.py#L156
> > been set to False? Shall I change that to True to make it appear/work?
> >
> > Build info
> > GRASS version: 8.0.dev 
> > Code revision: 563797cd3
> > Build date: 2022-01-11
> >
> > Thanks much,
> > Vero
> > --
> > Dra. Verónica Andreo
> > Investigadora Asistente de CONICET
> > Instituto Gulich (CONAE - UNC)
> > Centro Espacial Teófilo Tabanera (CETT)
> > Falda del Cañete - Córdoba, Argentina
> > +54 3547 40 int. 1153
> > https://veroandreo.gitlab.io/
> >
> > ___
> > grass-user mailing list
> > grass-user@lists.osgeo.org
> > https://lists.osgeo.org/mailman/listinfo/grass-user
>
> --
> Micha Silver
> Ben Gurion Univ.
> Sde Boker, Remote Sensing Lab
> cell: +972-523-665918
>
> ___
> grass-user mailing list
> grass-user@lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/grass-user
>


-- 
Prof. Carlos Henrique Grohmann
Institute of Energy and Environment - Univ. of São Paulo, Brazil
Editor-in-Chief - Brazilian Journal of Geology

http://carlosgrohmann.com
https://spamlab.github.io/
http://orcid.org/-0001-5073-5572

Can’t stop the signal.
___
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] Imported .tif maps lose all data

[Carlos Henrique Grohmann de Carvalho]

Hi Micha

You're right. That's what happens when I try to answer quickly while
keeping an eye on a 5-year old. Thanks.

Carlos



On Sat, Sep 25, 2021 at 3:32 PM Micha Silver  wrote:

> Hi Carlos:
>
>
> On 9/25/21 9:16 PM, Carlos Henrique Grohmann de Carvalho wrote:
>
> Rich
>
> I believe r.in.gdal will require the region to match each of the tifs
> you're importing. Try using r.import with the option extent set to
>
> "input", so it won't do region cropping
>
>
> I don't think that is correct. r.in.gdal is one of the (very few) raster
> modules that does NOT honor the current region settings. You can import a
> raster from "anywhere", but then, as always, you must set the current
> region to match that raster for further analysis.
>
>
> What's more, in your previous post you mentioned that r.import "overrides
> region settings". IMHO this is also not correct. r.import will reproject
> an input dataset to current Location's coordinate system, if necessary.
> Regarding the region settings, by default with the extent parameter left at
> "input", the whole input raster will be imported. If you set the the extent
> parameter to "region" then only that subset of the input raster that falls
> in the current region will be imported. But in both cases, the current
> region is not altered.
>
>
> Best,
>
> Micha
>
>
>
> C
>
>
>
> On Sat, Sep 25, 2021 at 1:43 PM Rich Shepard 
> wrote:
>
>> On Sat, 25 Sep 2021, Carlos Henrique Grohmann de Carvalho wrote:
>>
>> > How is your region settings? Do they match the DEMs? Maybe importing
>> with
>> > r.import will help, as it can override region settings
>>
>> Carlos,
>>
>> Initially it was the default region. But, I changed the region to a test
>> map, dtm_35, and re-ran r.report and r.stats. There's still no data in any
>> of the cells.
>>
>> Thanks,
>>
>> Rich
>> ___
>> grass-user mailing list
>> grass-user@lists.osgeo.org
>> https://lists.osgeo.org/mailman/listinfo/grass-user
>>
>
>
> --
> Prof. Carlos Henrique Grohmann
> Institute of Energy and Environment - Univ. of São Paulo, Brazil
> - Digital Terrain Analysis | GIS | Remote Sensing -
>
> http://carlosgrohmann.com
> http://orcid.org/-0001-5073-5572
> 
> Can’t stop the signal.
>
> ___
> grass-user mailing 
> listgrass-user@lists.osgeo.orghttps://lists.osgeo.org/mailman/listinfo/grass-user
>
> --
> Micha Silver
> Ben Gurion Univ.
> Sde Boker, Remote Sensing Lab
> cell: +972-523-665918
>
>

-- 
Prof. Carlos Henrique Grohmann
Institute of Energy and Environment - Univ. of São Paulo, Brazil
- Digital Terrain Analysis | GIS | Remote Sensing -

http://carlosgrohmann.com
http://orcid.org/-0001-5073-5572

Can’t stop the signal.
___
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] Imported .tif maps lose all data

[Carlos Henrique Grohmann de Carvalho]

Rich

I believe r.in.gdal will require the region to match each of the tifs
you're importing. Try using r.import with the option extent set to "input",
so it won't do region cropping

C



On Sat, Sep 25, 2021 at 1:43 PM Rich Shepard 
wrote:

> On Sat, 25 Sep 2021, Carlos Henrique Grohmann de Carvalho wrote:
>
> > How is your region settings? Do they match the DEMs? Maybe importing with
> > r.import will help, as it can override region settings
>
> Carlos,
>
> Initially it was the default region. But, I changed the region to a test
> map, dtm_35, and re-ran r.report and r.stats. There's still no data in any
> of the cells.
>
> Thanks,
>
> Rich
> ___
> grass-user mailing list
> grass-user@lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/grass-user
>


-- 
Prof. Carlos Henrique Grohmann
Institute of Energy and Environment - Univ. of São Paulo, Brazil
- Digital Terrain Analysis | GIS | Remote Sensing -

http://carlosgrohmann.com
http://orcid.org/-0001-5073-5572

Can’t stop the signal.
___
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] Imported .tif maps lose all data

[Carlos Henrique Grohmann de Carvalho]

How is your region settings? Do they match the DEMs? Maybe importing with
r.import will help, as it can override region settings

Carlos




On Sat, 25 Sep 2021, 11:43 Rich Shepard,  wrote:

> I have 77 LiDAR topo files in .tif format. What gdalinfo tells me about a
> representative file:
> $ gdalinfo columbia_2010_e_dtm_35.tif
> Driver: GTiff/GeoTIFF
> Files: columbia_2010_e_dtm_35.tif
> Size is 10745, 15264
> Coordinate System is:
> PROJCRS["NAD_1983_HARN_StatePlane_Washington_South_FIPS_4602_Feet",
>  BASEGEOGCRS["GCS_North_American_1983_HARN",
>  DATUM["North_American_1983_HARN",
>  ELLIPSOID["GRS_1980",6378137,298.257222101,
>  LENGTHUNIT["metre",1,
>  ID["EPSG",9001,
>  PRIMEM["Greenwich",0,
>  ANGLEUNIT["degree",0.0174532925199433,
>  ID["EPSG",9122,
>  CONVERSION["Lambert Conic Conformal (2SP)",
>  METHOD["Lambert Conic Conformal (2SP)",
>  ID["EPSG",9802]],
>  PARAMETER["Latitude of false origin",45.3,
>  ANGLEUNIT["degree",0.0174532925199433],
>  ID["EPSG",8821]],
>  PARAMETER["Longitude of false origin",-120.5,
>  ANGLEUNIT["degree",0.0174532925199433],
>  ID["EPSG",8822]],
>  PARAMETER["Latitude of 1st standard parallel",45.8,
>  ANGLEUNIT["degree",0.0174532925199433],
>  ID["EPSG",8823]],
>  PARAMETER["Latitude of 2nd standard parallel",47.3,
>  ANGLEUNIT["degree",0.0174532925199433],
>  ID["EPSG",8824]],
>  PARAMETER["Easting at false origin",50,
>  LENGTHUNIT["metre",1],
>  ID["EPSG",8826]],
>  PARAMETER["Northing at false origin",0,
>  LENGTHUNIT["metre",1],
>  ID["EPSG",8827]]],
>  CS[Cartesian,2],
>  AXIS["easting",east,
>  ORDER[1],
>  LENGTHUNIT["US survey foot",0.304800609601219,
>  ID["EPSG",9003]]],
>  AXIS["northing",north,
>  ORDER[2],
>  LENGTHUNIT["US survey foot",0.304800609601219,
>  ID["EPSG",9003
> Data axis to CRS axis mapping: 1,2
> Origin = (1512164.000,152311.000)
> Pixel Size = (3.000,-3.000)
> Metadata:
>AREA_OR_POINT=Area
>DataType=Generic
> Image Structure Metadata:
>COMPRESSION=LZW
>INTERLEAVE=BAND
> Corner Coordinates:
> Upper Left  ( 1512164.000,  152311.000) (121d 0' 8.57"W, 45d44'59.58"N)
> Lower Left  ( 1512164.000,  106519.000) (121d 0' 4.46"W, 45d37'27.52"N)
> Upper Right ( 1544399.000,  152311.000) (120d52'34.01"W, 45d45' 1.35"N)
> Lower Right ( 1544399.000,  106519.000) (120d52'30.94"W, 45d37'29.29"N)
> Center  ( 1528281.500,  129415.000) (120d56'19.49"W, 45d41'14.50"N)
> Band 1 Block=512x256 Type=Float32, ColorInterp=Gray
>Min=159.316 Max=1037.520
>Minimum=159.316, Maximum=1037.520, Mean=262.041, StdDev=157.494
>NoData Value=-3.40282306073709653e+38
>Overviews: 5373x7632, 2687x3816, 1344x1908, 672x954, 336x477, 168x239
>Metadata:
>  BandName=Band_1
>  RepresentationType=ATHEMATIC
>  STATISTICS_COVARIANCES=24804.31297971113
>  STATISTICS_MAXIMUM=1037.5201416016
>  STATISTICS_MEAN=262.04051816209
>  STATISTICS_MINIMUM=159.31565856934
>  STATISTICS_SKIPFACTORX=1
>  STATISTICS_SKIPFACTORY=1
>  STATISTICS_STDDEV=157.49385060919
>
> I started GRASS with a new location:
> grass -c /path/to/project/data/topography/columbia_2010_e_dtm_35.tif
> /data/grassdata/columbia_river_dtms/topography
>
> and imported all raw data; e.g.,
>
> r.in.gdal
> in=/path/to/project/data/topography/columbia_2010_e_dtm_35.tif out=dtm_4
> --o
>
> The problem is the imported maps have no elevation data. For example,
> r.report dtm_35
>   100%
>
> +-+
> | RASTER MAP CATEGORY REPORT
> |
> |LOCATION: new_nevins_dockSat Sep 25 07:28:25
> 2021|
>
> |-|
> |  north: 106729east: 1544180
>|
> |REGIONsouth:  60934west: 1511876
>|
> |  res:3res:3
>|
>
> |-|
> |MASK: none
>|
>
> |-|
> |MAP: (untitled) (dtm_35 in topography)
>|
>
> |-|
> |   Category Information | cell|
>  %  |
> |#|description   |count|
> cover|
>
> |-|
> |*|no data. . . . . . . . . . . . . . . . . . . .

[GRASS-user] i.atcorr - range/rescale issues

[Carlos Henrique Grohmann de Carvalho]

 Hello all

I'm preparing a tutorial on remote sensing with GRASS (already have 7
videos - in portuguese - here: https://www.youtube.com/c/SPAMLabUSP)

What I'm planning now is to run i.landsat.toar and i.atcorr, and then
compare the results with the surface reflectance products from Landsat
collection 2.

I got it all running but the results are a bit different and it seems that
the "range" and "rescale" options are the ones to tweak.

checking the metadata:

L2 before scale factor: {'min': '3783', 'max': '65535'}
L2 after scale factor: {'min': '-0.0959675', 'max': '1.6022125'}

L1 radiance at sensor: {'min': '-6.2', 'max': '191.6'}

L1 reflectance dos1: {'min': '0', 'max': '0.398754597845553'}

L1 reflectance 6S: {'min': '8.560894e-11', 'max': '255'}

I can use the min/max from radiance-at-sensor as the "range" in i.atcorr,
but what would be best for the "rescale" option? If I don't set it, I end
up with 0-255
values. How to get the floating-point result?

thank you for any input


Carlos


-- 
Prof. Carlos Henrique Grohmann
Institute of Energy and Environment - Univ. of São Paulo, Brazil
- Digital Terrain Analysis | GIS | Remote Sensing -

http://carlosgrohmann.com
http://orcid.org/-0001-5073-5572

Can’t stop the signal.
___
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] surface roughness

[Carlos Henrique Grohmann de Carvalho]

Hello, Carlos Grohmann here (I wrote the first version of
r.surface.roughness)

Surface roughness is one of those subjects that can mean very different
things depending on the context. I believe that you already have seen these
papers, but I recommend them for a better understanding of surface
roughness in geosciences.

Berti, M.; Corsini, A. & Daehne, A. Comparative analysis of surface
roughness algorithms for the identification of active landslides
Geomorphology, 2013, 182, 1-18
http://dx.doi.org/10.1016/j.geomorph.2012.10.022

Grohmann, C. H., Smith, M. J., & Riccomini, C. 2011. Multiscale Analysis of
Topographic Surface Roughness in the Midland Valley, Scotland. Geoscience
and Remote Sensing, IEEE Transactions on, 49(4): 1200–1213.
http://dx.doi.org/10.1109/TGRS.2010.2053546

The thing is, there are many definitions of surface roughness. Vector
Strength and Fisher's K are two of those, so these results *are* maps of
surface roughness. Just like standard deviation of slope, which, in my
opinion, is very good for geomorphology and morphotectonics. One nice thing
about std-slope is that you don't need any add-on to calculate it, just
make a slope map (r.slope.aspect) and then you can play with r.param.scale
to calculate the standard deviation using several window sizes and see
which one gives the better answer for your scale of analysis.


And while we're here, I'll leave a few suggestions of papers on
morphotectonics. The isobase method is quite nice. If you don't have access
to the papers at the publishers, in my website (carlosgrohmann.com), I have
links for preprints/postprints which have the same content just different
typesetting.

Grohmann, C. H., Riccomini, C., & Chamani, M. A. C. 2011. Regional scale
analysis of landform configuration with base-level (isobase) maps.
Hydrology and Earth System Sciences, 15(5): 1493–1504.
http://dx.doi.org/10.5194/hess-15-1493-2011

Grohmann, C. H., & Riccomini, C. 2009. Comparison of roving-window and
search-window techniques for characterising landscape morphometry.
Computers & Geosciences, 35: 2164–2169.
http://dx.doi.org/10.1016/j.cageo.2008.12.014

Grohmann, C. H., Riccomini, C., & Alves, F. M. 2007. SRTM-based
morphotectonic analysis of the Poços de Caldas Alkaline Massif,
southeastern Brazil. Computers & Geosciences, 33: 10–19.
http://dx.doi.org/10.1016/j.cageo.2006.05.002

Grohmann, C. H. 2004. Morphometric analysis in Geographic Information
Systems: applications of free software GRASS and R. Computers &
Geosciences, 30: 1055–1067.
http://dx.doi.org/10.1016/j.cageo.2004.08.002


Hope this helps

Carlos







On Sun, Feb 7, 2021 at 10:42 AM sacit mutlu 
wrote:

> good day dear admin. What are the steps we need to take to measure surface 
> roughness in the grass gis program. Is there a list you can send me about 
> this? Because I ran r.roughness.vector code and got only vector strength and 
> Fisher's K parameters. how do i get the surface roughness map ...
>
> ___
> grass-user mailing list
> grass-user@lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/grass-user
>


-- 
Prof. Carlos Henrique Grohmann
Institute of Energy and Environment - Univ. of São Paulo, Brazil
- Digital Terrain Analysis | GIS | Remote Sensing -

http://carlosgrohmann.com
http://orcid.org/-0001-5073-5572

Can’t stop the signal.
___
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] Fwd: Update on Geomorphometry 2020 - updates

[Carlos Henrique Grohmann de Carvalho]

Apologies for cross-posting.

Carlos Grohmann

-- Forwarded message -
From: Massimiliano Alvioli 
Date: Wed, Feb 19, 2020 at 8:49 AM
Subject: Update on Geomorphometry 2020 - updates
To: Massimiliano Alvioli 



Dear colleague,

you are receiving this email because you submitted either an abstract, or a
workshop
proposal, or both - to Geomorphometry 2020.

Let us share with you an updated POSTER of the conference and the LEAFLET
with
the list and description of the accepted workshops, attached to this email.
Please
help us announcing the news, by sharing the information with colleagues
and/or
on social media (see the bottom of this email).

Additional important UPDATES:

due to many requests & announcement of a SPECIAL ISSUE for the Proceedings
of
the Sixth Geomorphometry conference (June 22-26, Perugia, Italy), the
DEADLINE
for abstract submission is now extended to FEBRUARY 29, 2020:
http://geomorphometry2020.org/submission-of-short-papers

Proceedings will be published as a special issue of TRANSACTIONS IN GIS
(Wiley):
https://onlinelibrary.wiley.com/journal/14679671

The first day of the conference will be devoted to a series of hands-on
WORKSHOPS:
http://geomorphometry2020.org/scientific-program/workshops

The list of KEYNOTE SPEAKERS include the Lifetime Achievement Awards
recipient:
http://geomorphometry2020.org/keynote-speakers

A GRANT will be awarded to a young geomorphologist by the International
Association
of Geomorphologists:
http://www.geomorph.org/2020/02/call-for-grants-geomorphometry-conference-2020/


We look forward to seeing you next June in Perugia!


The organizing committee:

Massimiliano Alvioli, CNR IRPI - Perugia
Ivan Marchesini, CNR IRPI, Perugia
Laura Melelli, University of Perugia
Peter Guth, US Naval Academy - USA
Scott Peckham, University of Colorado - USA & rivix.com


ResearchGate:
https://www.researchgate.net/project/Geomorphometry-2020

Facebook:
https://it-it.facebook.com/events/1795086087290744/

LinkedIn:
https://www.linkedin.com/events/geomorphometry2020/

Twitter:
https://twitter.com/GeomorphomOrg

Official Website:
http://geomorphometry2020.org

-- 
Massimiliano Alvioli
CNR IRPI, Perugia
tel 0039 339 3459008
http://www.irpi.cnr.it






This message is for the recipients only. If received in error, please
notify the sender and read http://www.irpi.cnr.it/en/disclaimer/
Questo messaggio e' per l'attenzione dei soli destinatari. Se lo avete
ricevuto per sbaglio, siete pregati di avvisare il mittente e leggere
http://www.irpi.cnr.it/disclaimer/




-- 
Prof. Carlos Henrique Grohmann
Institute of Energy and Environment - Univ. of São Paulo, Brazil
- Digital Terrain Analysis | GIS | Remote Sensing -

http://carlosgrohmann.com
http://orcid.org/-0001-5073-5572

Can’t stop the signal.
___
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] What's a GCP?

[Carlos Henrique Grohmann de Carvalho]

Should be Ground Control Point(s)

C

On Wed, Sep 18, 2019 at 12:34 PM Rich Shepard 
wrote:

> The r.in.gdal manual page mentions GCPs. This is a new acronym for me and I
> cannot find a relevant definition in my web searches. For what do those
> initials stand?
>
> Regards,
>
> Rich
> ___
> grass-user mailing list
> grass-user@lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/grass-user



-- 
Prof. Carlos Henrique Grohmann
Institute of Energy and Environment - Univ. of São Paulo, Brazil
- Digital Terrain Analysis | GIS | Remote Sensing -

http://carlosgrohmann.com
http://orcid.org/-0001-5073-5572

Can’t stop the signal.
___
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] 7.7.dev: No raster map displays [FIXED]

[Carlos Henrique Grohmann de Carvalho]

Glad to hear that you got it working!

best

Carlos

On Wed, Jul 24, 2019 at 6:03 PM Rich Shepard 
wrote:

> On Wed, 24 Jul 2019, Carlos Henrique Grohmann de Carvalho wrote:
>
> > Well, you do have data there.. could it be as simple as setting a new
> > colortable?
>
> Carlos,
>
> I suspect not. It originally displayed in one set of colors and I changed
> that to the elevation color table. Now it's nothing but blank.
>
> Created a new, replacement location and imported the hdr.adf. It displays
> with the default color set so I will copy it from PERMANENT to topography,
> remove the former, and change the colortable to 'elevation.' The region is
> set (as the default) to the DEM and the analytical area is the new MASK.
>
> Thanks and best regards,
>
> Rich
> ___
> grass-user mailing list
> grass-user@lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/grass-user



-- 
Prof. Carlos Henrique Grohmann
Institute of Energy and Environment - Univ. of São Paulo, Brazil
- Digital Terrain Analysis | GIS | Remote Sensing -

http://carlosgrohmann.com
http://orcid.org/-0001-5073-5572

Can’t stop the signal.
___
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] 7.7.dev: No raster map displays

[Carlos Henrique Grohmann de Carvalho]

Well, you do have data there.. could it be as simple as setting a new
colortable?

On Wed, Jul 24, 2019 at 3:54 PM Rich Shepard 
wrote:

> On Wed, 24 Jul 2019, Carlos Henrique Grohmann de Carvalho wrote:
>
> > what does the metadata of the raster says (r.info)?
>
> Redacted version:
>
>
> ++
>   | Map:  xxx_bedem  Date: Wed Jul 24 11:30:58
> 2019|
>   | Mapset:   topography Login of Creator: rshepard
> |
>   | Location: someplace
> |
>   | DataBase: /data/grassdata
> |
>   | Title:
>  |
>   | Timestamp: none
> |
>
> ||
>   |
> |
>   |   Type of Map:  raster   Number of Categories: 0
>  |
>   |   Data Type:FCELL
> |
>   |   Rows: 15497
> |
>   |   Columns:  10222
> |
>   |   Total Cells:  158410334
> |
>   |Projection: NAD_1983_HARN_Oregon_Statewide_Lambert_Feet_Intl
> |
>   |N:  1294726.5S:  1248235.5   Res: 3
>  |
>   |E:   443401.5W:   412735.5   Res: 3
>  |
>   |   Range of data:min = -3.39  max = 1505.69
>  |
>   |
> |
>   |   Data Description:
> |
>   |generated by r.in.gdal
> |
>   |
> |
>   |   Comments:
> |
>   |r.in.gdal --overwrite -o input="/home/rshepard/projects/xxx/xxx\
> |
>   |/data/LDQ-45123B8/2009_OLC_North_Coast/Bare_Earth/\
>  |
>   |be45123b8/hdr.adf" output="xxx_bedem" memory=300 offset=0
> nu\   |
>   |m_digits=0
> |
>   |
> |
>
> ++
>
> Rich
> ___
> grass-user mailing list
> grass-user@lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/grass-user



-- 
Prof. Carlos Henrique Grohmann
Institute of Energy and Environment - Univ. of São Paulo, Brazil
- Digital Terrain Analysis | GIS | Remote Sensing -

http://carlosgrohmann.com
http://orcid.org/-0001-5073-5572

Can’t stop the signal.
___
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] 7.7.dev: No raster map displays

[Carlos Henrique Grohmann de Carvalho]

what does the metadata of the raster says (r.info)?

On Wed, Jul 24, 2019 at 3:33 PM Rich Shepard 
wrote:

> On Wed, 24 Jul 2019, Rich Shepard wrote:
>
> > Will do. Good idea because the default region is probably much too small.
>
> Set the region to the DEM and re-imported the data file. Still no display.
> Very strange. Please send other ideas for finding why raster maps no longer
> are visible.
>
> Regards,
>
> Rich
> ___
> grass-user mailing list
> grass-user@lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/grass-user



-- 
Prof. Carlos Henrique Grohmann
Institute of Energy and Environment - Univ. of São Paulo, Brazil
- Digital Terrain Analysis | GIS | Remote Sensing -

http://carlosgrohmann.com
http://orcid.org/-0001-5073-5572

Can’t stop the signal.
___
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] 7.7.dev: No raster map displays

[Carlos Henrique Grohmann de Carvalho]

Hi Rich

regarding the region setting, to adjust the resolution to 3m, use a "-a"
flag ('g.region vect=analytical_area res=3 -pa')

also try to set the region, than re-run any previous steps (like re-import
the raster), to be sure that any operation is being run with the proper
region settings.

C


On Wed, Jul 24, 2019 at 11:44 AM Rich Shepard 
wrote:

> On Wed, 24 Jul 2019, Carlos Henrique Grohmann de Carvalho wrote:
>
> > This might sound silly but have you checked your region? Perhaps the
> > region is not on the same place of your raster? Didi you tried 'zoom to
> > layer' or 'zoom to region'?
>
> Carlos
>
> Not a silly question at all. I tried 'zoom to layer' with no results.
>
> However, there is a mask for the vector map, 'analytical_area', which is
> much smaller the area covered by the DEM. The region is also set to the
> analytical_area map and this might be the issue. Changing the region to the
> raster DEM does not allow the map to display.
>
> While waiting for a solution to the display issue I tried running r.flow on
> the DEM and aspact maps. Grass refused because of a region resolution
> discrepancy between the analytical_area and DEM maps:
>
> DEM resolution: 3 3
> analytical_area resolution: 2.99763360 3.00134608
>
> 'g.region vect=analytical_area res=3 -p' did not change the floating point
> values to integer 3s.
>
> Thanks,
>
> Rich
> ___
> grass-user mailing list
> grass-user@lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/grass-user



-- 
Prof. Carlos Henrique Grohmann
Institute of Energy and Environment - Univ. of São Paulo, Brazil
- Digital Terrain Analysis | GIS | Remote Sensing -

http://carlosgrohmann.com
http://orcid.org/-0001-5073-5572

Can’t stop the signal.
___
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] 7.7.dev: No raster map displays

[Carlos Henrique Grohmann de Carvalho]

This might sound silly but have you checked your region? Perhaps the region
is not on the same place of your raster? Didi you tried 'zoom to layer' or
'zoom to region'?

C

On Wed, Jul 24, 2019 at 10:44 AM Rich Shepard 
wrote:

> This is a major problem that I need to resolve ASAP: the monitor (and GUI
> map display window) will not display any raster map. The console reports
> 100% loaded and the icon in the GUI map layer window adds the raster to the
> list, but neither produces a visible output of the raster map.
>
> Running a new 7.7.dev from a fresh git repo clone.
>
> Help needed,
>
> Rich
> ___
> grass-user mailing list
> grass-user@lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/grass-user



-- 
Prof. Carlos Henrique Grohmann
Institute of Energy and Environment - Univ. of São Paulo, Brazil
- Digital Terrain Analysis | GIS | Remote Sensing -

http://carlosgrohmann.com
http://orcid.org/-0001-5073-5572

Can’t stop the signal.
___
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] question on Landsat 8 processing

[Carlos Henrique Grohmann de Carvalho]

Indeed it needs ALL bands.

I used this module in class a couple of weeks ago. You need to import and
use all bands, including the panchromatic.

Carlos



On Tue, Jun 18, 2019 at 9:32 AM Veronica Andreo 
wrote:

> Seems to me that the module needs all Landsat bands. Try importing all of
> them
>
> Vero
>
> El lun., 17 jun. 2019 a las 15:54, Gabriel Cotlier ()
> escribió:
>
>> Dear Vero, Nikos and Markus,
>>
>> I'm trying to get an NDVI from different preprocessing atmospheric
>> correction methods DOS1, DOS2, DOS 2b, DOS3, DOS4. However, unsuccessfully,
>> I'm not being able to run my supposed workflow on GRASS GIS which is:
>>
>> Step 1. DN to TOA Radiance ( *i.landsat.toar*)
>> Step 2. TOA Radiance to Land Surface Reflectance  (*i.atcorr*)
>>
>> that thought maybe is the correct workflow for Landdsat 8 Level 1
>> products, maybe is not?
>>
>> Since I need it only for Band 4 and Band 5 I proceeded as follows:
>>
>> 1. in a grassdata directory named with the correspondent Landsat 8 scene
>> "LC82270832018074LGN00" I put Band 4 and Band 5 together with the MTL
>> metadata txt file.
>>
>> 2. set a "Location" for this data set using Band 4 data by reading
>> projection and Datum from Band 4.
>>
>> 3.started grass gis and d.rast for displaying the already automatically
>> loaded Band 4
>>
>> 5. import with r.in.gdal Band 5 using Band 5 complete name as was Band 4
>> loaded "LC82270832018074LGN00_B5"
>>
>> [image: image.png]
>> 6. I changed working directory to where the data is:
>>
>> [image: image.png]
>> 7. in command line (the withe one in grass interface) I tried to
>> *i.landsat.toar* for obtaining TOA Radiance as follows :
>>
>> i.landsat.toar -r input=LC82270832018074LGN00_B
>> output=TOAR_LC82270832018074LGN00_B metfile=LC82270832018074LGN00_MTL.txt
>> sensor=oli8 date=2018-03-15 sun_elevation=45.06637424
>> product_date=2018-03-20
>>
>> and the resulting error was:
>>
>> [image: image.png]
>>
>> I do not really understand what could be happening related to the bands
>> names?
>>
>> Thanks a lot I will appreciate a way to solve this problem. May target is
>> to get finally both bands 4 and 5 radiometrically callibrated and further
>> atmospherically corrected with *i.atcorr*  using different DOS methos.
>>
>> Thanks a lot again for the help and guidance.
>>
>> Gabriel
>>
>> ___
> grass-user mailing list
> grass-user@lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/grass-user



-- 
Prof. Carlos Henrique Grohmann
Institute of Energy and Environment - Univ. of São Paulo, Brazil
- Digital Terrain Analysis | GIS | Remote Sensing -

http://carlosgrohmann.com
http://orcid.org/-0001-5073-5572

Can’t stop the signal.
___
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] 3D unstructured Grid

[Carlos Henrique Grohmann de Carvalho]

That seems very interesting! Any plans to support PLY?

Carlos

On Tue, Oct 2, 2018 at 9:32 AM, Saber Razmjooei <
saber.razmjo...@lutraconsulting.co.uk> wrote:

> Hi,
>
> Sorry to hijack the discussion related to GRASS. But, we are working a new
> library to handle mesh data. MDAL is currently supported in QGIS. But the
> idea is to make it generic so other OS GIS application can handle the mesh
> layer through a generic library similar to GDAL. Here is more info about
> MDAL:
> https://lutraconsulting.github.io/mdal-web/
>
> You can find some example of the formats supported there too.
>
> Regards
> Saber
>
>
> On Tue, 2 Oct 2018 at 03:01, Vaclav Petras  wrote:
>
>> On Sat, Sep 29, 2018 at 4:45 PM Francois Chartier 
>> wrote:
>>
>>> Can Grass GIS create a 3D unstructured grid based on a 3D point data set
>>> ?
>>>
>>
>> All data and visualizations at the following website were created in
>> GRASS GIS including a 3D regular grid (3D raster in GRASS GIS) interpolated
>> from a set of 3D points (vector which is 3D, i.e. with z, in GRASS GIS).
>> Points were soil properties measurement at different depths, 3D raster is
>> then continuous representation of a selected property.
>>
>> https://geospatial.ncsu.edu/geoforall/multidim-geovis.html
>>
>> Module used for interpolation was v.vol.rst:
>>
>> https://grass.osgeo.org/grass74/manuals/v.vol.rst.html
>>
>> Vaclav
>> ___
>> grass-user mailing list
>> grass-user@lists.osgeo.org
>> https://lists.osgeo.org/mailman/listinfo/grass-user
>
>
>
> --
> Saber Razmjooei
> www.lutraconsulting.co.uk
> +44 (0)7568 129733
>
> ___
> grass-user mailing list
> grass-user@lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/grass-user
>



-- 
Prof. Carlos Henrique Grohmann
Institute of Energy and Environment - Univ. of São Paulo, Brazil
- Digital Terrain Analysis | GIS | Remote Sensing -

http://carlosgrohmann.com
http://orcid.org/-0001-5073-5572

Can’t stop the signal.
___
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] 3D unstructured Grid

[Carlos Henrique Grohmann de Carvalho]

What do you mean? Like a 3D mesh?
You can use Meshlab, ClouCompare or Paraview to do that, and then exporting
it back to GRASS as 3D faces.

C

On Sat, Sep 29, 2018 at 5:44 PM, Francois Chartier 
wrote:

> Hi,
>
> Can Grass GIS create a 3D unstructured grid based on a 3D point data set ?
>
> thanks
> F
>
> ___
> grass-user mailing list
> grass-user@lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/grass-user
>



-- 
Prof. Carlos Henrique Grohmann
Institute of Energy and Environment - Univ. of São Paulo, Brazil
- Digital Terrain Analysis | GIS | Remote Sensing -

http://carlosgrohmann.com
http://orcid.org/-0001-5073-5572

Can’t stop the signal.
___
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] m.nviz.image - simulating aerial (nadir-looking) images

[Carlos Henrique Grohmann de Carvalho]

Hello all

Now that m.nviz.image is working (thanks devs!), I'm trying to simulate
aerial (UAV) images.

Is this possible?

Playing around with the parameters, I can't seem to get a nadir-looking
image. If I set the height of the 'camera' to something like 500 meters
(it's a coastal area, elevations are all below 40m), the resulting image
looks like a nadir observation, but as the height becomes smaller, say 120
meters, the surface is tilted (much like in Google Earth) and the
observation is an oblique view of the surface.

I tried setting the 'focus' parameter to the central coordinate of the
region, but the result is the same.

One issue that might be causing some confusion in my head is that some
parameters expect 'model coordinates' (0-1), while others expect map
units...

Another point is the 'perspective' param. What does it mean? Should I use
the field of view of my real camera here?

This is the region:

projection: 1 (UTM)
zone:   -22
datum:  wgs84
ellipsoid:  wgs84
north:  6900503
south:  6900383
west:   732294
east:   732474
nsres:  0.5
ewres:  0.5
rows:   240
cols:   360
cells:  86400


center point:
732384.0, 6900443.0

I want the view to be on the center of the image, so position=0.5,0.5
(right?)

set focus to center point (in map units, from SW corner - corner of what?
the raster? the region?)

focus=90.0,60.0,11.72

This is the full command (in python)

>>> grass.run_command('m.nviz.image', elevation_map=elev, color_map=shade,
mode='fine', resolution_fine='1', \
... style='surface', position='0.50,0.50', height='120', focus=foc,
perspective=cam_fov, size=img_size, \
... format='ppm', output='test', overwrite=True, verbose=True)


For testing, this are links to the 'elev' and 'shade' rasters:

https://www.dropbox.com/s/ou0ft6srs9w7e6t/elev.tif?dl=0
https://www.dropbox.com/s/urcf4m60rvtsi4f/shade.tif?dl=0


thanks

Carlos













-- 
Prof. Carlos Henrique Grohmann
Institute of Energy and Environment - Univ. of São Paulo, Brazil
- Digital Terrain Analysis | GIS | Remote Sensing -

http://carlosgrohmann.com
http://orcid.org/-0001-5073-5572

Can’t stop the signal.
___
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] help with m.nviz.image - no output

[Carlos Henrique Grohmann de Carvalho]

Hello all

I'm trying to create some images with m.nviz.image, but all I get as output
is a 11 byte ppm...

Some of the options are quite puzzling in terms of mixing units:

position=x,y
Viewpoint position (x,y model coordinates)

height=value
Viewpoint height (in map units)

focus=x,y,z
Focus to point on surface (from SW corner in map units)

So I have to convert my position coordinates from map units (say meters in
UTM) to 0-1, I can use height in meters, and also need to do some math to
get the position of the focus point in meters from the SW corner...

Also, regardless of the output size, all I get is an unreadable ppm with 11
bytes. With format=tif, I get an error creating the file
(TIFFScanlineSize64: Computed scanline size is zero.).

Opening the ppm file into a text editor, I see that all that was written
was the header, and that the image size was 0,0:

P6
0 0
255

navigating trough the code, I see that m.nviz.image calls GS_write_ppm to
write the file, which is in lib/ogsf/gsd_img_ppm.c:


int GS_write_ppm(const char *name)
{
unsigned int x;
int y;
unsigned int xsize, ysize;
FILE *fp;
unsigned char *pixbuf;

gsd_getimage(&pixbuf, &xsize, &ysize);

if (NULL == (fp = fopen(name, "w"))) {
G_warning(_("Unable to open file <%s> for writing"), name);
return (1);
}

fprintf(fp, "P6\n%d %d\n255\n", xsize, ysize);

for (y = ysize - 1; y >= 0; y--) {
for (x = 0; x < xsize; x++) {
   unsigned char r = pixbuf[(y * xsize + x) * 4 + 0];
   unsigned char g = pixbuf[(y * xsize + x) * 4 + 1];
   unsigned char b = pixbuf[(y * xsize + x) * 4 + 2];

   fputc((int)r, fp);
   fputc((int)g, fp);
   fputc((int)b, fp);
}

}
G_free(pixbuf);
fclose(fp);

return (0);
}


GS_write_ppm calls gsd_getimage to get the actual image data and size.

So it seems to me that gsd_getimage is not getting the image?




The goal is to reproduce a nadir-looking aerial view (airplane, UAV).

Is this even possible?

thanks





-- 
Prof. Carlos Henrique Grohmann
Institute of Energy and Environment - Univ. of São Paulo, Brazil
- Digital Terrain Analysis | GIS | Remote Sensing -

http://carlosgrohmann.com
http://orcid.org/-0001-5073-5572

Can’t stop the signal.
___
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user