The gridliner handles some of that to turn off the visibility on some ticks based upon that criteria. There could be some useful code in there to build from: https://github.com/SciTools/cartopy/blob/019f0130ec9320cad8ac2532697afb51db0f97a9/lib/cartopy/mpl/gridliner.py#L1276-L1280
4 replies
Hello!
When I plot a high resolution satellite image (tiff format) in imshow and save as jpeg with very high dpi (1000) the quality of the image is retained.
But when I plot the same image with cartopy the image is resampled (I assume onto the resolution of the lat-lon grid) and I loose alot of image quality.
Please could some one advice me how to plot a high resolution sat. image in cartopy retaining the original image quality?
Thank you very much!
There is a
regrid_shape argument that is likely what you want to look into: https://scitools.org.uk/cartopy/docs/latest/reference/generated/cartopy.mpl.geoaxes.GeoAxes.html#cartopy.mpl.geoaxes.GeoAxes.imshow
4 replies
Hi there, I have a question regarding the natural earth resources. The Naciscdn server has been down for days now. I've seen that the fix is to use the aws link instead but I like to avoid amazon whenever I can so I was wondering whether someone was in touch with nacis and would now if the server would get back up some day or whether it's been taken down or something
Sorry, those were supposed to be in that thread @greglucas. Weird behavior here. I'm pretty sure image 2 just results in a blank image because when I take off the cartopy features and set_extent, there's still nothing
Tried with pcolormesh (passing rgba tuples to color kwarg) and the transform works, but all these dots appear. Frustrating day!
Hm... My initial guess is that your image 2 may have to do with you using img1.size as the regrid_shape argument. That will be nm4 for an (n, m) RGBA image, not just the shape of the output, so I'm guessing that is producing a really large output array that may cause issues with the queries...? Just a stab in the dark there. I'd try just setting it to a large value at first (2000) and see if that works and increase appropriately to your problem. For the second pcolormesh idea, I'm not sure that pcolormesh takes RGBA tuples...? So it could be doing some odd resampling based on the strides in your array to produce those dots, or reordering the data in an odd way.
1 reply
Hi all, i want to test some modification on cartopy, so i forked the project and did python setup.py install to install it localy.
I got a error ModuleNotFoundError: No module named 'cartopy._version'
I saw that in v0.18 there was a file _version.py in the cartopy project, but there isn"t anymore. Do i have to do extra step to use it locally ? Thanks
There is an environment.yml file listing a bunch of dependencies in the base repo now that you can use to create a new conda environment that will be able to compile/install Cartopy for testing. You can also just use that as a reference for which packages you'd need. There are also some install_requires and requirements.txt files for different dependencies
I'm experiencing some weirdness with pcolormesh. I'm trying to pass the kwarg shading='nearest', which per the matplotlib docs will work when C has the same dimensions of X and Y. It seems like the cartopy transformation adds an element to X and Y because I get the error:
Even though I double checked and the shape of C (300, 400) is the same as my input X and Y. I've tested and the same line that works with plain matplotlib fails with cartopy.
TypeError: Dimensions of C (300, 400) are incompatible with X (401) and/or Y (301); see help(pcolormesh)Even though I double checked and the shape of C (300, 400) is the same as my input X and Y. I've tested and the same line that works with plain matplotlib fails with cartopy.
I don't think this makes much difference for anyone, except the GitHub UI
Hi all, im trying to replicate the following srtm shading example https://scitools.org.uk/cartopy/docs/v0.14/examples/srtm_shading.html however when i run the code i get this error
UserWarning: SRTM requires an account set up and log in to access. Use of this Downloader is likely to fail with HTTP 401 errors.
warnings.warn('SRTM requires an account set up and log in to access.. I have looked at the documentation thoroughly and i can't find anything regarding this. Is there a new version of the srtm_shading example? Thanks.
There is no direct API for setting that, but it uses regular
urlopen, so you can set a password override using that
Elliot, any interest in PRing your branch that removes the final PROJ dependencies in trace? QuLogic/cartopy@d109aac
I still think the benefits outweigh the decreased speed of spherical interpolators. Especially with all of the recent PROJ installation issues being opened.
I still think the benefits outweigh the decreased speed of spherical interpolators. Especially with all of the recent PROJ installation issues being opened.
Hello, someone how can I solve this problem?
Python 3.8.10 (default, Nov 26 2021, 20:14:08)
Type 'copyright', 'credits' or 'license' for more information
IPython 7.24.0 -- An enhanced Interactive Python. Type '?' for help.
In [1]: import cartopy.crs as ccrs
ImportError Traceback (most recent call last)
<ipython-input-1-139fa5ad2567> in <module>
----> 1 import cartopy.crs as ccrs
~/.local/lib/python3.8/site-packages/cartopy/init.py in <module>
102 # Commonly used sub-modules. Imported here to provide end-user
103 # convenience.
--> 104 import cartopy.crs
105 import cartopy.feature # noqa: F401 (flake8 = unused import)
~/.local/lib/python3.8/site-packages/cartopy/crs.py in <module>
20 from shapely.prepared import prep
21
---> 22 from cartopy._crs import (CRS, Geodetic, Globe, PROJ4_VERSION,
23 WGS84_SEMIMAJOR_AXIS, WGS84_SEMIMINOR_AXIS)
24 from cartopy._crs import Geocentric # noqa: F401 (flake8 = unused import)
ImportError: libproj.so.19: cannot open shared object file: No such file or directory
Hi, How can I install Cartopy? I get the error:
pip3 install cartopy==0.20.2 --no-warn-script-location
Defaulting to user installation because normal site-packages is not writeable
Collecting cartopy==0.20.2
Using cached Cartopy-0.20.2.tar.gz (10.8 MB)
Installing build dependencies ... done
Getting requirements to build wheel ... error
ERROR: Command errored out with exit status 1:
command: /usr/bin/python3 /home/leo/.local/lib/python3.8/site-packages/pip/_vendor/pep517/in_process/_in_process.py get_requires_for_build_wheel /tmp/tmpxc5rsxyg
cwd: /tmp/pip-install-edx6ys1f/cartopy_faed090c1dd0401bb4dbbd520024780d
Complete output (3 lines):
setup.py:117: UserWarning: Unable to determine GEOS version. Ensure you have 3.7.2 or later installed, or installation may fail.
warnings.warn(
pip3 install cartopy==0.20.2 --no-warn-script-location
Defaulting to user installation because normal site-packages is not writeable
Collecting cartopy==0.20.2
Using cached Cartopy-0.20.2.tar.gz (10.8 MB)
Installing build dependencies ... done
Getting requirements to build wheel ... error
ERROR: Command errored out with exit status 1:
command: /usr/bin/python3 /home/leo/.local/lib/python3.8/site-packages/pip/_vendor/pep517/in_process/_in_process.py get_requires_for_build_wheel /tmp/tmpxc5rsxyg
cwd: /tmp/pip-install-edx6ys1f/cartopy_faed090c1dd0401bb4dbbd520024780d
Complete output (3 lines):
setup.py:117: UserWarning: Unable to determine GEOS version. Ensure you have 3.7.2 or later installed, or installation may fail.
warnings.warn(
Proj 8.0.0 must be installed.
Make sure you use your Linux package manager to install libproj-dev and libgeos-dev (or their equivalents) and that libproj-dev is at least version 8. If that's not possible, I'd recommend installing using conda, following the official install instructions: https://scitools.org.uk/cartopy/docs/latest/installing.html#conda-pre-built-binaries
Yes, I work with plotting satellite data and passes. Have you binned your data into a grid and are trying to plot that? Or trying to plot the ground tracks of satellites? My guess is that if you are not able to transform anything is that you're using the incorrect transform for your data coordinates. If you have lat/lon ground tracks use the PlateCarree/Geodetic input to the function call, and then whatever projection you want in the Axes() creation calls. But, remember the
transform=data for the plotting function calls.
What I have is the TLE of a satellite (or the actual ephemeris data, depending on the satellite) and would like to plot a field of view coverage area (assuming some fixed field of view), the circle gets plotted, but it doesn't get properly transformed to account for the geodetic Earth
You may want to look at the Nightshade implementation which follows the subsolar point and handles some curvatures in that way. https://scitools.org.uk/cartopy/docs/latest/gallery/lines_and_polygons/nightshade.html https://github.com/SciTools/cartopy/blob/main/lib/cartopy/feature/nightshade.py
Thanks, I will look at that. I also found this example, https://scitools.org.uk/cartopy/docs/v0.10/examples/geostationary.html, which seems to take a circular image and uses the Geostationary reference frame to transform it into a more expected shape. I'm guessing this won't work for an arbitrary satellite location, but maybe I can figure out something by looking at that code
Hmm, NearsidePerspective does seem to be what I want... it looks like it can take the lat/lon/alt of the satellite and do a transformation. The projected coordinates are supposed to be an x,y measured from the origin of the plane tangent to the Earth directly below the satellite. That should be easy enough to calculate. I'm guessing thats x,y in meters?
I'm confused now... I think it does happen for the Circle() patch.
import matplotlib.pyplot as plt
import cartopy.crs as ccrs
fig, ax = plt.subplots(subplot_kw={"projection": ccrs.Robinson()})
ax.coastlines()
patch = plt.Circle((140, 20), 20, color='g',
alpha=0.5, transform=ccrs.PlateCarree().as_geodetic())
ax.add_patch(patch)
ax.set_global()
plt.show()
So it seems like my issue was that I was using the PlateCarree projection with a Geodetic transformation which didn't alter the shape of the circle. Using Robinson and Geodetic did alter the circle in a way consistent with the coordinate transformation, but also made me realize that I will either need to figure out NearsidePerspective or perform a more detailed calculation for determining the view coverage as a circle doesn't just become a rounded rectangle by coordinate transformation, but instead by field of view mapping
looks like NearsidePerspective gives me what I need (or at least something that looks a lot like what I need). The only tricky thing is that if your radius is too large it just won't plot anything. Luckily it seems my calculation of the circle from before is enough to limit the radius to something it will correctly plot