Re: [GRASS-user] Reprojecting to x,y grid
Hi Markus, On 2020-11-17 at 14:40 -08, Markus Neteler wrote... > On Tue, Nov 17, 2020 at 7:06 PM Ken Mankoff wrote: >> I have a data set provided on an (x,y) grid, where each coordinate >> does have a paired (lon,lat). I'd like to work on the native (x,y) >> grid, and bring some other geo-referenced data there. How do I do >> this? > > Here 2.5 wild guesses: > > If it is about "pushing around" a raster map (no squeezing): > - https://grass.osgeo.org/grass78/manuals/r.region.html > > If a rotation is needed: > - https://grass.osgeo.org/grass7/manuals/addons/i.rotate.html > > If it is about a real transform: > - convert the raster matrix to its vector center point representation; then > - https://grass.osgeo.org/grass78/manuals/v.transform.html > - rasterize again Thanks for getting me thinking down this path. This isn't about the raster. I want to work in the raster XY location. I want to reproject a vector there. I know the (lon,lat) of every (x,y). I think maybe I need to solve this outside of GRASS - I can take the (lon,lat) of each point of the geo-referenced vector, and look up the nearest (lon,lat) in the XY location, and set those as the new coordinates for the vector. I'm not yet sure what the difference is between that method and using v.rectify. I would build the control points using the same method: Take the (lon,lat) of each point of the geo-referenced vector, look up the nearest (lon,lat) in the XY location, from there jump to the (x,y) coord in the XY location, and use those as my control points. I then rectify (un-project?) the vector to the XY coords. Is the only difference that the second method would "fill in" points that I don't map between (lon,lat) and (x,y), while the first method requires that I update every point? Thanks, -k. ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Reprojecting to x,y grid
On Tue, Nov 17, 2020 at 7:06 PM Ken Mankoff wrote: > > Hi GRASS list, > > I have a data set provided on an (x,y) grid, where each coordinate does have > a paired (lon,lat). I'd like to work on the native (x,y) grid, and bring some > other geo-referenced data there. How do I do this? Here 2.5 wild guesses: If it is about "pushing around" a raster map (no squeezing): - https://grass.osgeo.org/grass78/manuals/r.region.html If a rotation is needed: - https://grass.osgeo.org/grass7/manuals/addons/i.rotate.html If it is about a real transform: - convert the raster matrix to its vector center point representation; then - https://grass.osgeo.org/grass78/manuals/v.transform.html - rasterize again But again, just some almost random ideas. best, Markus ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] [GRASS-dev] let's set a meeting time during OSGeo Code Sprint!
Hi, út 17. 11. 2020 v 20:43 odesílatel Veronica Andreo napsal: > https://www.timeanddate.com/worldclock/converter.html?iso=20201118T20=1440=204=48=207=485=671 > > Let's use the dedicated channel: https://meet.jit.si/GRASS great, maybe one hour is not so much. What about two hours: 19 - 21 UTC? Ma -- Martin Landa http://geo.fsv.cvut.cz/gwiki/Landa http://gismentors.cz/mentors/landa ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] let's set a meeting time during OSGeo Code Sprint!
Seems we have a winner so far: https://www.timeanddate.com/worldclock/converter.html?iso=20201118T20=1440=204=48=207=485=671 Let's use the dedicated channel: https://meet.jit.si/GRASS See you tomorrow! Vero El mar., 17 nov. 2020 a las 13:39, Veronica Andreo () escribió: > Hi everyone, > > I just created a doodle to see if we can meet during these sprint days. > Considering our previous meetings this year I have proposed time slots > between 7 and 10 pm (Amsterdam time zone). I'm not sure how it will appear > for you, though :) > > Here's the doodle: > https://doodle.com/poll/rzek6s7vd7tivta7?utm_source=poll_medium=link > > Cheers, > Vero > ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] grass python problem
Hi Everyone Thanks for your response, for some reason it works now. Thanks Ming Micha Silver 于2020年11月17日周二 上午10:44写道: > > On 11/17/20 6:03 AM, ming han wrote: > > Hi Everyone > > Hope this email finds you well. > I can run r.cross in console with following command: > "r.cross -z --overwrite --quiet input=Connect_Lake@PERMANENT > ,str_grass_r@PERMANENT output=test" > > But I got error by run r.cross in python with following command: > grass.run_command('r.cross', input = ['str_grass_r@PERMANENT > ','Connect_Lake@PERMANENT'],output = 'test', quiet = True) > > What is the correct format to use r.cross function in Python? and how > to obtain the output table generated by r.cross? > > > This worked fine for me also, using GRASS 7.8.4 in the nc_basin_spm > LOCATION: > > > In [2]: import grass.script as gscript > > In [3]: gscript.run_command("r.cross", input=["landuse@PERMANENT > ","basins@PERMANENT"], output="lu_bas", overwrite=True) > r.cross: STEP 1 ... > 100% > r.cross: STEP 2 ... > r.cross: STEP 3 ... > 100% > Creating support files for ... > 89 categories > Out[3]: 0 > > > Just to be clear, you do not get a table output, rather a raster. In your > case you should see a new raster "test". > > > > Thanks > Ming > > ___ > 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 > > ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Reprojecting to x,y grid
Hi GRASS list, I have a data set provided on an (x,y) grid, where each coordinate does have a paired (lon,lat). I'd like to work on the native (x,y) grid, and bring some other geo-referenced data there. How do I do this? I've successfully worked in rotated pole coordinates (e.g., [1]) but because I don't know what (x,y) coord contains the pole, I haven't figured out how to adapt this to the rotated pole situation as in [1]. Anyway, this isn't a rotated pole... The coordinates for the data look like this: # x, y, lon, lat 1,1,-46.930,59.105 2,1,-46.843,59.127 3,1,-46.757,59.148 4,1,-46.670,59.169 5,1,-46.584,59.191 6,1,-46.497,59.212 7,1,-46.410,59.233 293,441,-10.886,84.128 294,441,-10.415,84.120 295,441,-9.944,84.112 296,441,-9.476,84.104 297,441,-9.008,84.095 298,441,-8.542,84.086 299,441,-8.077,84.076 300,441,-7.614,84.066 Meaning there are 300 columns and 441 rows, each with a (lon,lat). I can easily import this into any projection I want with m.proj, but I'd like to work in this XY coord: grass -c XY ./G g.region s=0 n=441 w=0 e=300 res=1 -pas The problem is I have one other vector file (in EPSG:4326 or EPSG:3413) that I need to use with this grid-coordinate data. How do I reproject or import that? Thanks, -k. [1] https://lists.osgeo.org/pipermail/grass-user/2011-October/062180.html ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] grass python problem
On 11/17/20 6:03 AM, ming han wrote: Hi Everyone Hope this email finds you well. I can run r.cross in console with following command: "r.cross -z --overwrite --quiet input=Connect_Lake@PERMANENT,str_grass_r@PERMANENT output=test" But I got error by run r.cross in python with following command: grass.run_command('r.cross', input = ['str_grass_r@PERMANENT','Connect_Lake@PERMANENT'],output = 'test', quiet = True) What is the correct format to use r.cross function in Python? and how to obtain the output table generated by r.cross? This worked fine for me also, using GRASS 7.8.4 in the nc_basin_spm LOCATION: In [2]: import grass.script as gscript In [3]: gscript.run_command("r.cross", input=["landuse@PERMANENT","basins@PERMANENT"], output="lu_bas", overwrite=True) r.cross: STEP 1 ... 100% r.cross: STEP 2 ... r.cross: STEP 3 ... 100% Creating support files for ... 89 categories Out[3]: 0 Just to be clear, you do not get a table output, rather a raster. In your case you should see a new raster "test". Thanks Ming ___ 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
Re: [GRASS-user] [OSGeo-Discuss] Reminder: OSGeo Virtual Community Sprint 2020 (17-19 November 2020)
Hi, I have created a draft board in GitHub: OSGeo Virtual Community Sprint 2020 board https://github.com/OSGeo/grass/projects/4 and added a series of open pull requests which we may tackle in the sprint. Feel free to modify! cheers Markus On Tue, Nov 17, 2020 at 10:29 AM Markus Neteler wrote: > > Starting today - please join - open to anyone, no programming > knowledge required! > > - Main zoom channel for kick off, presentations, free discussion and > daily reports > https://zoom.us/j/95817956179?pwd=VGROK0NVK2lTaDg4eHJkdW1IdWNlUT09 > - GRASS GIS related channel: https://meet.jit.si/GRASS > - IRC chat during the event: #tosprint on Server: irc.freenode.net >(https://webchat.freenode.net/) > ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] let's set a meeting time during OSGeo Code Sprint!
Hi everyone, I just created a doodle to see if we can meet during these sprint days. Considering our previous meetings this year I have proposed time slots between 7 and 10 pm (Amsterdam time zone). I'm not sure how it will appear for you, though :) Here's the doodle: https://doodle.com/poll/rzek6s7vd7tivta7?utm_source=poll_medium=link Cheers, Vero ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Fwd: [OSGeo-Discuss] Reminder: OSGeo Virtual Community Sprint 2020 (17-19 November 2020)
Starting today - please join - open to anyone, no programming knowledge required! - Main zoom channel for kick off, presentations, free discussion and daily reports https://zoom.us/j/95817956179?pwd=VGROK0NVK2lTaDg4eHJkdW1IdWNlUT09 - GRASS GIS related channel: https://meet.jit.si/GRASS - IRC chat during the event: #tosprint on Server: irc.freenode.net (https://webchat.freenode.net/) -- Forwarded message - From: Victoria Rautenbach Date: Mon, Nov 16, 2020 at 2:25 PM Subject: [OSGeo-Discuss] Reminder: OSGeo Virtual Community Sprint 2020 (17-19 November 2020) To: osgeolist , Dear All We would like to remind everyone of the OSGeo Virtual Community Sprint 2020 that will take place 17-19 November 2020. Information about schedule and meetings have been added to the wiki: https://wiki.osgeo.org/wiki/OSGeo_Virtual_Community_Sprint_2020 The sprint will be open every day from 8am UTC to midnight UTC, with a couple of special sessions scheduled (meeting links available in the wiki): A short kick off meeting will start at 8am UTC every day in the main room A welcome session is planned for Tue 17 Nov 1pm UTC in the main room At the end of each day, a daily report session is planned around 10pm UTC The code sprint is also open to anyone! If you are interested in getting involved in any of the projects, join one of the daily sessions to get started. Any questions can be posted also to any of the available communication channels, https://wiki.osgeo.org/wiki/OSGeo_Virtual_Community_Sprint_2020#Communication Happy coding, OSGeo Board ___ Discuss mailing list disc...@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/discuss ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] [GRASS-dev] grass python problem
Please provide more information: GRASS version Python version exact error message I did not observe any errors while running command you provided (of course I used different map names to match map names I had on my system). Māris. otrd., 2020. g. 17. nov., plkst. 06:04 — lietotājs ming han () rakstīja: > > Hi Everyone > > Hope this email finds you well. > I can run r.cross in console with following command: > "r.cross -z --overwrite --quiet > input=Connect_Lake@PERMANENT,str_grass_r@PERMANENT output=test" > > But I got error by run r.cross in python with following command: > grass.run_command('r.cross', input = > ['str_grass_r@PERMANENT','Connect_Lake@PERMANENT'],output = 'test', quiet = > True) > > What is the correct format to use r.cross function in Python? and how to > obtain the output table generated by r.cross? > > Thanks > Ming ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user