Where communities thrive


  • Join over 1.5M+ people
  • Join over 100K+ communities
  • Free without limits
  • Create your own community
People
Repo info
Activity
  • Sep 30 18:51
    dopplershift commented #1940
  • Sep 30 09:57
    JoMaHoPo commented #1940
  • Sep 30 00:58
    AlxndrLhr closed #2086
  • Sep 29 18:47
    dopplershift commented #2080
  • Sep 29 11:39
    jorisvandenbossche commented #2080
  • Sep 28 08:19
    jorisvandenbossche synchronize #2080
  • Sep 27 20:51
    jorisvandenbossche commented #2080
  • Sep 27 20:50
    jorisvandenbossche synchronize #2080
  • Sep 27 18:58
    dopplershift closed #879
  • Sep 27 18:58
    dopplershift commented #879
  • Sep 27 18:54
    dopplershift closed #2090
  • Sep 27 18:54
    dopplershift commented #2090
  • Sep 27 16:16
    cameronbracken commented #879
  • Sep 27 16:12
    cameronbracken commented #879
  • Sep 27 14:09
    4goettma opened #2090
  • Sep 27 07:35
    QuLogic closed #2088
  • Sep 27 07:30
    PBrockmann commented #2088
  • Sep 27 05:46
    jrenaud90 opened #2089
  • Sep 26 23:28
    dopplershift commented #2087
  • Sep 26 23:27
    dopplershift commented #2088
Ryan May
@dopplershift
@Pheosics I'm not sure what your radius represents? The math in NearsidePerspective is calculating based on tangent lines for a given satellite height.
Miles Engel
@Pheosics
By radius I mean the radius of the circle that is formed when you draw a cone from the satellite down to a tangent plane which is centered right below the satellite. Of course if the satellite is far away or has a large enough look angle, the circle formed will be larger than the cross section of the Earth, so you have to limit it. In an ideal world NearsidePerspective would handle that, but it seems like it doesn't and so if you create a circle with too large of radius and try and use NearsidePerspective to display it it will just not display anything.
Ryan May
@dopplershift
So Nearside perspective is calculating purely where the line from the satellite is tangent to the surface of the earth—so the physical limit one can see without looking through rock. It’s not using any kind of field of view limit from the satellite.
George Duffy
@CowboyGeorge

Hello. I just installed Cartopy on a new MacOS Pro 13", running Monterey. I just tried to copy the tutorial script "https://scitools.org.uk/cartopy/docs/latest/gallery/lines_and_polygons/always_circular_stereo.html", and I'm getting an error at the line where i'm supposed to apply cfeature.LAND,
"PermissionError: [Errno 13] Permission denied: '~/.local/share/cartopy'"

I don't have "cartopy" directory in /.local/share/, which seems like problem. I've tried reinstalling to no avail (using the conda forge), and I haven't seen any similar problems on stack overflow or anything similar. Any ideas on what could be going wrong?

Ryan May
@dopplershift
Cartopy should create the directory cartopy inside ~/.local/share, where ~ is your user home directory. I'd check the permissions on ~/.local/share.
George Duffy
@CowboyGeorge

Thanks Ryan,

I did find that my ~/.local/share directory had root permissions, which is also strange since according to the person who set up my permission there wasn't even a ~/.local directory on my computer when received it. The directory was empty, so I deleted the directory and made a new one (~./local/share) with user permission, but when I installed it with conda it just tells me that "All requested packages already installed", with no new files created inside ~./local/share.

Greg Lucas
@greglucas
I'm not sure whether conda installs anything to that directory when you install, but if you run that example again, it should give you some "DownloadWarnings" that populate that directory with the needed files.
George Duffy
@CowboyGeorge
Tutorial is working! Thank you!
George Duffy
@CowboyGeorge
image.png
image.png
anuappu
@anuappu
Hello: I am a a beginner in Python/cartopy etc so pardon my ignorance. I am encountering an error when I utilize the cointourf utility. Toplology exception side conflict error.
Snippet of my code
 ax = plt.axes(projection=ccrs.LambertConformal(central_longitude=10.0, central_latitude=40.0, standard_parallels=(40.0, 40.0)))
I am a bit confused on what the ax.extent transform should be:
ax.set_extent([-30, 40, 25, 70],crs=ccrs.???())
and finally sc=plt.contourf(xnew, ynew, perc, clevs,
transform=ccrs.PlateCarree(),colors=colours,extend='max')
where xnew is lon and lat in degrees and perc is probability values that I am computing.
anuappu
@anuappu
The code works 95% of the time but once in while fails with a an error "Toplology exception side conflict error at .."
The input data and the calculated probability values seem correct.
Appreciate the feedback.
Greg Lucas
@greglucas
Yes, it looks like there are some bad transformations going on creating invalid geometries at some point in the chain. No one has addressed that issue yet. You can try "transform_first=True" to the contouring algorithm which may lead to nicer geometries... See: https://scitools.org.uk/cartopy/docs/latest/gallery/scalar_data/contour_transforms.html
anuappu
@anuappu
Thank you @greglucas for the quick response. The "transform_first=True" did not work. However, may I ask, you mention that "bad transformations going on creating invalid geometries at some point in the chain" could it be in the way the ax.extent is defined?
Greg Lucas
@greglucas
Yes, it certainly could be, you can try adjusting the limits to see if that helps.
anuappu
@anuappu
Thanks, I will try that and keep you posted.
anuappu
@anuappu
@greglucas I have not have had any success with the issue by altering the ax.extent limits. Is it fair to say the it could be the conflict between the projection of LambertConformal and the transformation to PlateCaree? Also, the area where I am trying to project crosses the 0 longitude (domain includes both the hemispheres). Any suggestions would really be appreciated.
Greg Lucas
@greglucas
Yes, it is very likely bad transformation with some points in a Polygon, and that object being cut/stitched back together incorrectly. No one has found the cause of this yet in that issue, so unfortunately there isn't much we can do until the root cause is determined. One other thought is to try adjusting your contouring levels, which will produce different sets of coordinates.
anuappu
@anuappu
Thanks @greglucas
Alexandre ROULEAU
@arouleau56

Hello everyone, I trie to apply the second cell of this link https://www.geophysique.be/2014/02/25/shaded-relief-map-in-python/

See the code here :

Subplot 2: STRM Map

plt.figure(figsize=(15, 10))
ax = plt.subplot(111, projection=ccrs.PlateCarree())

elev, crs, extent = srtm_composite(12, 47, 1, 1)
elev = np.ma.masked_less_equal(elev,0,copy=False)
plt.imshow(elev, extent=extent, transform=crs,
cmap='gist_earth', origin='lower')

However I got these message : "HTTP Error 401: Unauthorized".
Do you have an idea howto solve it?
Thank you :)

Ryan May
@dopplershift
I'm not even sure how that code could run that far? srtm_composite was removed in CartoPy 0.12 according to the changelog
The SRTM example was removed in #847 due to login issues, but if you have access maybe that version will work for you?
Alex Crosby
@acrosby
is there a problem with the recent versions' packaging for pypi/pip
getting pip complaints about the version metadata not matching the stated version, eventually rolls back to an earlier version that it is happy with
Greg Lucas
@greglucas
What version are you trying? This should be fixed on v0.20.3 that was just released. SciTools/cartopy#2053
Alex Crosby
@acrosby
oh I see, i think i saw it going through 0.20.x<0.20.3 but this was just prior to 0.20.3, thanks!
Elliott Sales de Andrade
@QuLogic
lots of image comparison failures with mpl 3.6, after patching the removed API
haven't checked yet if main is better with the switch to pytest-mpl
Elliott Sales de Andrade
@QuLogic
Ryan May
@dopplershift
oof
Houssam Halaby
@hhalaby
Hello,
Any plans to release cartopy without proj dependency (#2069) ? would be very useful for our installation process.
edit: I see it's in the 0.21 milestone. Any rough idea when that might be ?