Join GitHub today
GitHub is home to over 28 million developers working together to host and review code, manage projects, and build software together.Sign up
The geos projection pictures how a geostationary satellite scans the earth at regular scanning angle intervals.
In order to project using the geos projection you can do the following:
proj +proj=geos +h=35785831.0
The required argument
h is the viewing point (satellite position) height above the earth.
Other arguments include:
lon_0: the subsatellite longitude point
sweep: the sweep angle axis of the viewing instrument. Can be
yis the default). (new in 4.8.0)
The projection coordinate relate to the scanning angle by the following simple relation:
scanning_angle (radians) = projection_coordinate / h
The viewing instrument on-board geostationary satellites described by this projection have a two-axis gimbal viewing geometry. This means that the different scanning positions are obtained by rotating the gimbal along a N/S axis (or ''y'') and a E/W axis (or ''x'').
In the image above, the outer-gimbal axis, or sweep-angle axis, is the N/S axis (
y) while the inner-gimbal axis, or fixed-angle axis, is the E/W axis (
This example represents the scanning geometry of the Meteosat series satellite. However, the GOES satellite series use the opposite scanning geometry, with the E/W axis (
x) as the sweep-angle axis, and the N/S (
y) as the fixed-angle axis.
sweep argument is used to tell
proj.4 which on which axis the outer-gimbal is rotating. The possible values are
y being the default. Thus, the scanning geometry of the Meteosat series satellite should take
x, and GOES should take