Re: [Qgis-user] vector layer only displays when zoomed out to layer extent
So for the first question I asked below, I was referring to your statement that "When you cross 360 degrees, some projections have difficulty dealing with that.” My lines don’t technically cross 360 deg, but they cross 0 deg, as the data values for the line positions are all > -360 and < 360 deg. So they cross at 0 deg with values like -0.345, 33.00 to 0.231, 32.94. I was able to save the densification layer to a file as you suggested, but when I do that the lines in the densified file will only display when the zoom is for the full extent of the layer like the original file. Only if the densified lines are in a temporary scratch layer will they display at other zoom levels. Regarding the other suggestion "It's likely a corrupt spatial index on the file. Try removing and recreating it.”. I don’t know how to find, remove and recreate a corrupt spatial index. Any instructions on how to do this would be appreciated. Thanks, Lee > On Dec 17, 2019, at 1:54 PM, Nicolas Cadieux > wrote: > > Hi, > > Not sure for question one. Have you tried creating a spatial index as others > have suggested? When you run the densify tool, If you don't specify the name, > the output will be a temporary layer. You should be able to do a save as for > the layer. Remember to close the temporary file after exporting it or else > you will get that message. > > Nicolas > > On 2019-12-17 12:51 p.m., Lee Eddington wrote: >> If my layer extent is > -360 and < 360 am I really crossing 360? I’m >> crossing 0 from negative values to positive values. Is that the same? >> >> Also, regarding not being able to save the layer I created with >> densification, when I try to close QGIS it says I have a ’scratch’ layer >> that will be lost, which would be this layer created by densification. How >> do you save a scratch layer? >> >> Best, >> Lee >> >> >>> On Dec 16, 2019, at 7:43 PM, Nicolas Cadieux >>> wrote: >>> >>> Hi, >>> When you cross 360 degrees, some projections have difficulty dealing with >>> that. Try not crossing that or try splitting the shapes that cross it. >>> Nicolas >>> >>>> Le 16 déc. 2019 à 22:32, Lee Eddington a écrit >>>> : >>>> >>>> I’m not sure that’s the reason as I have plenty of points in my lines. >>>> Also, I saved the densified layer to disk, but the saved layer has the >>>> same display problem as the original layer. >>>> >>>> Lee >>>> >>>> >>>>> On Dec 16, 2019, at 8:35 AM, Nicolas Cadieux >>>>> wrote: >>>>> >>>>> Hi, >>>>> >>>>> If you have a very long straight line represented by only two point, >>>>> reprojection operates only on those two points (called vertices or >>>>> nodes). The lines just represent a straight connections between two >>>>> point. The problem is that the earth is round (contrary to some popular >>>>> belief :) and those long straight lines need to be bent depending on the >>>>> projection being used. Since reprojection and zooming only operates on >>>>> vertices or nodes, densification adds nodes to the line giving more >>>>> points to calculate and thus, less errors. >>>>> Cheer! >>>>> Nicolas >>>>> >>>>> >>>>> >>>>>>> Le 16 déc. 2019 à 10:23, Lee Eddington a >>>>>>> écrit : >>>>>> Hi Nicolas, >>>>>> >>>>>> Vector densification worked. I went with the default 1 vertices to add >>>>>> and it created a new layer that now displays at any zoom level. Can you >>>>>> explain why my original layer doesn’t display? >>>>>> >>>>>> Thanks, >>>>>> Lee >>>>>> >>>>>> >>>>>>> On Dec 14, 2019, at 3:13 PM, Nicolas Cadieux >>>>>>> wrote: >> . >> ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] vector layer only displays when zoomed out to layer extent
If my layer extent is > -360 and < 360 am I really crossing 360? I’m crossing 0 from negative values to positive values. Is that the same? Also, regarding not being able to save the layer I created with densification, when I try to close QGIS it says I have a ’scratch’ layer that will be lost, which would be this layer created by densification. How do you save a scratch layer? Best, Lee > On Dec 16, 2019, at 7:43 PM, Nicolas Cadieux > wrote: > > Hi, > When you cross 360 degrees, some projections have difficulty dealing with > that. Try not crossing that or try splitting the shapes that cross it. > Nicolas > >> Le 16 déc. 2019 à 22:32, Lee Eddington a écrit : >> >> I’m not sure that’s the reason as I have plenty of points in my lines. >> Also, I saved the densified layer to disk, but the saved layer has the same >> display problem as the original layer. >> >> Lee >> >> >>> On Dec 16, 2019, at 8:35 AM, Nicolas Cadieux >>> wrote: >>> >>> Hi, >>> >>> If you have a very long straight line represented by only two point, >>> reprojection operates only on those two points (called vertices or nodes). >>> The lines just represent a straight connections between two point. The >>> problem is that the earth is round (contrary to some popular belief :) and >>> those long straight lines need to be bent depending on the projection being >>> used. Since reprojection and zooming only operates on vertices or nodes, >>> densification adds nodes to the line giving more points to calculate and >>> thus, less errors. >>> Cheer! >>> Nicolas >>> >>> >>> >>>>> Le 16 déc. 2019 à 10:23, Lee Eddington a >>>>> écrit : >>>> >>>> Hi Nicolas, >>>> >>>> Vector densification worked. I went with the default 1 vertices to add >>>> and it created a new layer that now displays at any zoom level. Can you >>>> explain why my original layer doesn’t display? >>>> >>>> Thanks, >>>> Lee >>>> >>>> >>>>> On Dec 14, 2019, at 3:13 PM, Nicolas Cadieux >>>>> wrote: ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] vector layer only displays when zoomed out to layer extent
I’m not sure that’s the reason as I have plenty of points in my lines. Also, I saved the densified layer to disk, but the saved layer has the same display problem as the original layer. Lee > On Dec 16, 2019, at 8:35 AM, Nicolas Cadieux > wrote: > > Hi, > > If you have a very long straight line represented by only two point, > reprojection operates only on those two points (called vertices or nodes). > The lines just represent a straight connections between two point. The > problem is that the earth is round (contrary to some popular belief :) and > those long straight lines need to be bent depending on the projection being > used. Since reprojection and zooming only operates on vertices or nodes, > densification adds nodes to the line giving more points to calculate and > thus, less errors. > Cheer! > Nicolas > > > >> Le 16 déc. 2019 à 10:23, Lee Eddington a écrit : >> >> Hi Nicolas, >> >> Vector densification worked. I went with the default 1 vertices to add and >> it created a new layer that now displays at any zoom level. Can you explain >> why my original layer doesn’t display? >> >> Thanks, >> Lee >> >> >>> On Dec 14, 2019, at 3:13 PM, Nicolas Cadieux >>> wrote: >>> >>> Hi, >>> >>> Try matching the project CRS to the layer CRS. Make sure you don’t have a >>> scale rendering option. If it’s not that, try vector densification. Could >>> work. >>> >>> Nicolas >>> >>>>> Le 14 déc. 2019 à 16:38, Lee Eddington a >>>>> écrit : >>>> >>>> I have a vector layer of lines that span -360 to 360 deg. They display >>>> fine when I’m zoomed out full, but if I zoom in to any degree they >>>> disappear. >>>> >>>> I’m using 3.4.14 on a Mac. >>>> >>>> Lee >>>> >>>> ___ >>>> Qgis-user mailing list >>>> Qgis-user@lists.osgeo.org >>>> List info: https://lists.osgeo.org/mailman/listinfo/qgis-user >>>> Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user >> ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] vector layer only displays when zoomed out to layer extent
Hi Nicolas, Vector densification worked. I went with the default 1 vertices to add and it created a new layer that now displays at any zoom level. Can you explain why my original layer doesn’t display? Thanks, Lee > On Dec 14, 2019, at 3:13 PM, Nicolas Cadieux > wrote: > > Hi, > > Try matching the project CRS to the layer CRS. Make sure you don’t have a > scale rendering option. If it’s not that, try vector densification. Could > work. > > Nicolas > >> Le 14 déc. 2019 à 16:38, Lee Eddington a écrit : >> >> I have a vector layer of lines that span -360 to 360 deg. They display >> fine when I’m zoomed out full, but if I zoom in to any degree they disappear. >> >> I’m using 3.4.14 on a Mac. >> >> Lee >> >> ___ >> Qgis-user mailing list >> Qgis-user@lists.osgeo.org >> List info: https://lists.osgeo.org/mailman/listinfo/qgis-user >> Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
[Qgis-user] vector layer only displays when zoomed out to layer extent
I have a vector layer of lines that span -360 to 360 deg. They display fine when I’m zoomed out full, but if I zoom in to any degree they disappear. I’m using 3.4.14 on a Mac. Lee ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] writing vector data for QGIS input from a python script using osgeo ogr
Etienne, Being a newer user I don’t know about the QGIS API. I will take a look. Thank you, Lee > On Jan 21, 2019, at 9:03 PM, Etienne Trimaille > wrote: > > QGIS is using OGR in the background to read vector data. So every format > supported by OGR is fine for QGIS. You can write geojson, geopackage, sqlite, > shapefile, it's fine for OGR API(write) and QGIS Desktop (read). > But are you sure you want to use the OGR API to write your file? IMHO, QGIS > API is easier. > > Le lun. 21 janv. 2019 à 23:49, Nicolas Cadieux <mailto:nicolas.cadi...@archeotec.ca>> a écrit : > QGIS can read wkt. This could help. > https://en.m.wikipedia.org/wiki/Well-known_text > <https://en.m.wikipedia.org/wiki/Well-known_text> > Nicolas > > Le 21 janv. 2019 à 20:26, Lee Eddington <mailto:lee.w.edding...@gmail.com>> a écrit : > >> I’d like to write vector line data in lat/lon coordinates to a file that >> QGIS can read and display. I have the osgeo ogr module, but can’t see to >> find the right documentation on how to use it for this simple task. >> >> Thanks, >> Lee >> >> ___ >> Qgis-user mailing list >> Qgis-user@lists.osgeo.org <mailto:Qgis-user@lists.osgeo.org> >> List info: https://lists.osgeo.org/mailman/listinfo/qgis-user >> <https://lists.osgeo.org/mailman/listinfo/qgis-user> >> Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user >> <https://lists.osgeo.org/mailman/listinfo/qgis-user>___ > Qgis-user mailing list > Qgis-user@lists.osgeo.org <mailto:Qgis-user@lists.osgeo.org> > List info: https://lists.osgeo.org/mailman/listinfo/qgis-user > <https://lists.osgeo.org/mailman/listinfo/qgis-user> > Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user > <https://lists.osgeo.org/mailman/listinfo/qgis-user> ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] writing vector data for QGIS input from a python script using osgeo ogr
Nicolas, Thank you for the link. I was able to use the info to create a wkt file that worked. But I also had to use a blog post to help me to figure out how to add the proper header: https://gis.stackexchange.com/questions/107152/how-to-create-wkt-file-readable-by-qgis-2-4-0 <https://gis.stackexchange.com/questions/107152/how-to-create-wkt-file-readable-by-qgis-2-4-0> Do you know where the official documentation is that explains that? I’m sure it’s out there, but I couldn’t find it. Thanks, Lee > On Jan 21, 2019, at 8:48 PM, Nicolas Cadieux > wrote: > > QGIS can read wkt. This could help. > https://en.m.wikipedia.org/wiki/Well-known_text > <https://en.m.wikipedia.org/wiki/Well-known_text> > Nicolas > > Le 21 janv. 2019 à 20:26, Lee Eddington <mailto:lee.w.edding...@gmail.com>> a écrit : > >> I’d like to write vector line data in lat/lon coordinates to a file that >> QGIS can read and display. I have the osgeo ogr module, but can’t see to >> find the right documentation on how to use it for this simple task. >> >> Thanks, >> Lee >> >> ___ >> Qgis-user mailing list >> Qgis-user@lists.osgeo.org <mailto:Qgis-user@lists.osgeo.org> >> List info: https://lists.osgeo.org/mailman/listinfo/qgis-user >> <https://lists.osgeo.org/mailman/listinfo/qgis-user> >> Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user >> <https://lists.osgeo.org/mailman/listinfo/qgis-user> ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
[Qgis-user] writing vector data for QGIS input from a python script using osgeo ogr
I’d like to write vector line data in lat/lon coordinates to a file that QGIS can read and display. I have the osgeo ogr module, but can’t see to find the right documentation on how to use it for this simple task. Thanks, Lee ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
[Qgis-user] compute spatial derivatives from raster bands
I’m fairly new to QGIS and would like to compute spatial derivatives (d/dx, d/dy) of raster band data. Is there a tool, computational interface or plugin that does this already? Thanks, Lee ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] how to properly reference longitude values in different datasets
Coming back to this…I did some reading on MDAL and tried the instructions found at: https://www.lutraconsulting.co.uk/blog/2018/10/18/mdal/ <https://www.lutraconsulting.co.uk/blog/2018/10/18/mdal/> trying to load the GFS dataset as a mesh instead of a raster, but I get the following error message: Tools::BufferedFileWriter::open: Cannot open file.: unspecified upstream_category error I also tried loading the ECMWF_ERA-40_subset.nc dataset linked to the blog post and get an error message saying it is an invalid data source. I’m using QGIS version 3.4.1-Madeira on MacOSX 10.14.2 Lee > On Dec 21, 2018, at 2:46 PM, Lee Eddington wrote: > > As a raster. > > Sent from my iPhone > > On Dec 21, 2018, at 12:09 AM, Saber Razmjooei > <mailto:saber.razmjo...@lutraconsulting.co.uk>> wrote: > >> Did you open your layer as Mesh layer or Raster? >> >> Cheers >> Saber >> >> On Fri, 21 Dec 2018 at 05:59, Lee Eddington > <mailto:lee.w.edding...@gmail.com>> wrote: >> That worked. >> >> Thanks, >> Lee >> >> >>> On Dec 20, 2018, at 6:06 AM, Saber Razmjooei >>> >> <mailto:saber.razmjo...@lutraconsulting.co.uk>> wrote: >>> >>> Climate Data Operator (CDO): >>> https://code.mpimet.mpg.de/projects/cdo/ >>> <https://code.mpimet.mpg.de/projects/cdo/> >>> >>> cdo sellonlatbox,-180,180,-90,90 input.grib output.grib >>> >>> Regards >>> Saber >>> >>> On Thu, 20 Dec 2018 at 13:00, Lee Eddington >> <mailto:lee.w.edding...@gmail.com>> wrote: >>> Can you point me to the external tools? >>> >>> Thanks, >>> Lee >>> >>> >>>> On Dec 20, 2018, at 12:46 AM, Saber Razmjooei >>>> >>> <mailto:saber.razmjo...@lutraconsulting.co.uk>> wrote: >>>> >>>> This is a known issue and can be fixed with some external tools. >>>> I have filed a bug: >>>> https://issues.qgis.org/issues/20853 <https://issues.qgis.org/issues/20853> >>>> And hopefully get around to fix it in 3.6 >>>> >>>> Regards >>>> Saber >>>> >>>> On Thu, 20 Dec 2018 at 01:44, Lee Eddington >>> <mailto:lee.w.edding...@gmail.com>> wrote: >>>> I’m pretty new to QGIS and am trying to figure out how to modify longitude >>>> values of layers. I was able to add a raster layer in GRIB2 format (NCEP >>>> global meteorological data). The longitude values range from 0 to 360 >>>> degrees. I’m able to view and do calculations on the data with no >>>> problem. I also added shaded relief and coastline data from Natural >>>> Earth. This data has longitude values from -180 to 180 degrees. So from >>>> -180 to 0 degrees I only see the Natural Earth data, from 0 to 180 degrees >>>> I see both datasets and from 180 to 360 degrees I only see the GRIB2 data. >>>> I want both datasets to be referenced the same with regards to longitude >>>> (-180 to 180, 0 to 360 or whatever). It seems like this should be a >>>> fairly easy problem to fix by either adding 180 degrees to longitudes from >>>> the Natural Earth data or subtracting 180 from the GRIB2 longitudes, but I >>>> can’t figure out how to do it. >>>> >>>> I’m using QGIS version 3.4.1-Madeira. The CRS I’m working in is what QGIS >>>> put the GRIB2 data into when I added it: * Generated CRS (+proj=longlat >>>> +a=6371229 +b=6371229 +no_defs) >>>> >>>> Thanks, >>>> Lee >>>> >>>> ___ >>>> Qgis-user mailing list >>>> Qgis-user@lists.osgeo.org <mailto:Qgis-user@lists.osgeo.org> >>>> List info: https://lists.osgeo.org/mailman/listinfo/qgis-user >>>> <https://lists.osgeo.org/mailman/listinfo/qgis-user> >>>> Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user >>>> <https://lists.osgeo.org/mailman/listinfo/qgis-user> >>>> >>>> -- >>>> Saber Razmjooei >>>> www.lutraconsulting.co.uk <http://www.lutraconsulting.co.uk/> >>>> +44 (0)7568 129733 >>> >>> >>> >>> -- >>> Saber Razmjooei >>> www.lutraconsulting.co.uk <http://www.lutraconsulting.co.uk/> >>> +44 (0)7568 129733 >> >> >> >> -- >> Saber Razmjooei >> www.lutraconsulting.co.uk <http://www.lutraconsulting.co.uk/> >> +44 (0)7568 129733 ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] how to properly reference longitude values in different datasets
As a raster. Sent from my iPhone > On Dec 21, 2018, at 12:09 AM, Saber Razmjooei > wrote: > > Did you open your layer as Mesh layer or Raster? > > Cheers > Saber > >> On Fri, 21 Dec 2018 at 05:59, Lee Eddington >> wrote: >> That worked. >> >> Thanks, >> Lee >> >> >>> On Dec 20, 2018, at 6:06 AM, Saber Razmjooei >>> wrote: >>> >>> Climate Data Operator (CDO): >>> https://code.mpimet.mpg.de/projects/cdo/ >>> >>> cdo sellonlatbox,-180,180,-90,90 input.grib output.grib >>> >>> Regards >>> Saber >>> >>>> On Thu, 20 Dec 2018 at 13:00, Lee Eddington >>>> wrote: >>>> Can you point me to the external tools? >>>> >>>> Thanks, >>>> Lee >>>> >>>> >>>>> On Dec 20, 2018, at 12:46 AM, Saber Razmjooei >>>>> wrote: >>>>> >>>>> This is a known issue and can be fixed with some external tools. >>>>> I have filed a bug: >>>>> https://issues.qgis.org/issues/20853 >>>>> And hopefully get around to fix it in 3.6 >>>>> >>>>> Regards >>>>> Saber >>>>> >>>>>> On Thu, 20 Dec 2018 at 01:44, Lee Eddington >>>>>> wrote: >>>>>> I’m pretty new to QGIS and am trying to figure out how to modify >>>>>> longitude values of layers. I was able to add a raster layer in GRIB2 >>>>>> format (NCEP global meteorological data). The longitude values range >>>>>> from 0 to 360 degrees. I’m able to view and do calculations on the data >>>>>> with no problem. I also added shaded relief and coastline data from >>>>>> Natural Earth. This data has longitude values from -180 to 180 degrees. >>>>>> So from -180 to 0 degrees I only see the Natural Earth data, from 0 to >>>>>> 180 degrees I see both datasets and from 180 to 360 degrees I only see >>>>>> the GRIB2 data. I want both datasets to be referenced the same with >>>>>> regards to longitude (-180 to 180, 0 to 360 or whatever). It seems like >>>>>> this should be a fairly easy problem to fix by either adding 180 degrees >>>>>> to longitudes from the Natural Earth data or subtracting 180 from the >>>>>> GRIB2 longitudes, but I can’t figure out how to do it. >>>>>> >>>>>> I’m using QGIS version 3.4.1-Madeira. The CRS I’m working in is what >>>>>> QGIS put the GRIB2 data into when I added it: * Generated CRS >>>>>> (+proj=longlat +a=6371229 +b=6371229 +no_defs) >>>>>> >>>>>> Thanks, >>>>>> Lee >>>>>> >>>>>> ___ >>>>>> Qgis-user mailing list >>>>>> Qgis-user@lists.osgeo.org >>>>>> List info: https://lists.osgeo.org/mailman/listinfo/qgis-user >>>>>> Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user >>>>> >>>>> >>>>> -- >>>>> Saber Razmjooei >>>>> www.lutraconsulting.co.uk >>>>> +44 (0)7568 129733 >>>> >>> >>> >>> -- >>> Saber Razmjooei >>> www.lutraconsulting.co.uk >>> +44 (0)7568 129733 >> > > > -- > Saber Razmjooei > www.lutraconsulting.co.uk > +44 (0)7568 129733 ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] how to properly reference longitude values in different datasets
That worked. Thanks, Lee > On Dec 20, 2018, at 6:06 AM, Saber Razmjooei > wrote: > > Climate Data Operator (CDO): > https://code.mpimet.mpg.de/projects/cdo/ > <https://code.mpimet.mpg.de/projects/cdo/> > > cdo sellonlatbox,-180,180,-90,90 input.grib output.grib > > Regards > Saber > > On Thu, 20 Dec 2018 at 13:00, Lee Eddington <mailto:lee.w.edding...@gmail.com>> wrote: > Can you point me to the external tools? > > Thanks, > Lee > > >> On Dec 20, 2018, at 12:46 AM, Saber Razmjooei >> > <mailto:saber.razmjo...@lutraconsulting.co.uk>> wrote: >> >> This is a known issue and can be fixed with some external tools. >> I have filed a bug: >> https://issues.qgis.org/issues/20853 <https://issues.qgis.org/issues/20853> >> And hopefully get around to fix it in 3.6 >> >> Regards >> Saber >> >> On Thu, 20 Dec 2018 at 01:44, Lee Eddington > <mailto:lee.w.edding...@gmail.com>> wrote: >> I’m pretty new to QGIS and am trying to figure out how to modify longitude >> values of layers. I was able to add a raster layer in GRIB2 format (NCEP >> global meteorological data). The longitude values range from 0 to 360 >> degrees. I’m able to view and do calculations on the data with no problem. >> I also added shaded relief and coastline data from Natural Earth. This data >> has longitude values from -180 to 180 degrees. So from -180 to 0 degrees I >> only see the Natural Earth data, from 0 to 180 degrees I see both datasets >> and from 180 to 360 degrees I only see the GRIB2 data. I want both datasets >> to be referenced the same with regards to longitude (-180 to 180, 0 to 360 >> or whatever). It seems like this should be a fairly easy problem to fix by >> either adding 180 degrees to longitudes from the Natural Earth data or >> subtracting 180 from the GRIB2 longitudes, but I can’t figure out how to do >> it. >> >> I’m using QGIS version 3.4.1-Madeira. The CRS I’m working in is what QGIS >> put the GRIB2 data into when I added it: * Generated CRS (+proj=longlat >> +a=6371229 +b=6371229 +no_defs) >> >> Thanks, >> Lee >> >> ___ >> Qgis-user mailing list >> Qgis-user@lists.osgeo.org <mailto:Qgis-user@lists.osgeo.org> >> List info: https://lists.osgeo.org/mailman/listinfo/qgis-user >> <https://lists.osgeo.org/mailman/listinfo/qgis-user> >> Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user >> <https://lists.osgeo.org/mailman/listinfo/qgis-user> >> >> -- >> Saber Razmjooei >> www.lutraconsulting.co.uk <http://www.lutraconsulting.co.uk/> >> +44 (0)7568 129733 > > > > -- > Saber Razmjooei > www.lutraconsulting.co.uk <http://www.lutraconsulting.co.uk/> > +44 (0)7568 129733 ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] how to properly reference longitude values in different datasets
Can you point me to the external tools? Thanks, Lee > On Dec 20, 2018, at 12:46 AM, Saber Razmjooei > wrote: > > This is a known issue and can be fixed with some external tools. > I have filed a bug: > https://issues.qgis.org/issues/20853 <https://issues.qgis.org/issues/20853> > And hopefully get around to fix it in 3.6 > > Regards > Saber > > On Thu, 20 Dec 2018 at 01:44, Lee Eddington <mailto:lee.w.edding...@gmail.com>> wrote: > I’m pretty new to QGIS and am trying to figure out how to modify longitude > values of layers. I was able to add a raster layer in GRIB2 format (NCEP > global meteorological data). The longitude values range from 0 to 360 > degrees. I’m able to view and do calculations on the data with no problem. > I also added shaded relief and coastline data from Natural Earth. This data > has longitude values from -180 to 180 degrees. So from -180 to 0 degrees I > only see the Natural Earth data, from 0 to 180 degrees I see both datasets > and from 180 to 360 degrees I only see the GRIB2 data. I want both datasets > to be referenced the same with regards to longitude (-180 to 180, 0 to 360 or > whatever). It seems like this should be a fairly easy problem to fix by > either adding 180 degrees to longitudes from the Natural Earth data or > subtracting 180 from the GRIB2 longitudes, but I can’t figure out how to do > it. > > I’m using QGIS version 3.4.1-Madeira. The CRS I’m working in is what QGIS > put the GRIB2 data into when I added it: * Generated CRS (+proj=longlat > +a=6371229 +b=6371229 +no_defs) > > Thanks, > Lee > > ___ > Qgis-user mailing list > Qgis-user@lists.osgeo.org <mailto:Qgis-user@lists.osgeo.org> > List info: https://lists.osgeo.org/mailman/listinfo/qgis-user > <https://lists.osgeo.org/mailman/listinfo/qgis-user> > Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user > <https://lists.osgeo.org/mailman/listinfo/qgis-user> > > -- > Saber Razmjooei > www.lutraconsulting.co.uk <http://www.lutraconsulting.co.uk/> > +44 (0)7568 129733 ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
[Qgis-user] how to properly reference longitude values in different datasets
I’m pretty new to QGIS and am trying to figure out how to modify longitude values of layers. I was able to add a raster layer in GRIB2 format (NCEP global meteorological data). The longitude values range from 0 to 360 degrees. I’m able to view and do calculations on the data with no problem. I also added shaded relief and coastline data from Natural Earth. This data has longitude values from -180 to 180 degrees. So from -180 to 0 degrees I only see the Natural Earth data, from 0 to 180 degrees I see both datasets and from 180 to 360 degrees I only see the GRIB2 data. I want both datasets to be referenced the same with regards to longitude (-180 to 180, 0 to 360 or whatever). It seems like this should be a fairly easy problem to fix by either adding 180 degrees to longitudes from the Natural Earth data or subtracting 180 from the GRIB2 longitudes, but I can’t figure out how to do it. I’m using QGIS version 3.4.1-Madeira. The CRS I’m working in is what QGIS put the GRIB2 data into when I added it: * Generated CRS (+proj=longlat +a=6371229 +b=6371229 +no_defs) Thanks, Lee ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
Re: [Qgis-user] /bin/sh: gdal_merge.py: command not found
Not sure how to remove this post, but I figured it out the problem. gdal_merge.py was in /Library/Frameworks/Python.framework/Versions/3.6/lib/python3.6/site-packages/GDAL-2.3.2.dist-info/scripts. Previously when I searched for gdal_merge.py I couldn’t find it. I was able to find it when searching as root. Reran the merge command and it wouldn’t work until I changed the permissions on the scripts to executable. Lee > On Nov 20, 2018, at 3:57 PM, Lee Eddington wrote: > > Hi, > > I’m new to QGIS and to this list. I’m going through the Training Manual and > have a problem doing a raster merge: > > > > > > I’m not sure how to get gdal_merge.py (or any other scripts I might be > missing) installed properly. I added custom environment variables for PATH > to get around problems creating a virtual raster and reprojecting (warping) a > raster. I have gdal in > /Library/Frameworks/Python.framework/Versions/3.6/lib/python3.6/site-packages > (2.3.2). > > Thanks, > Lee > ___ Qgis-user mailing list Qgis-user@lists.osgeo.org List info: https://lists.osgeo.org/mailman/listinfo/qgis-user Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user
