Differences between revisions 9 and 10
⇤ ← Revision 9 as of 2024-10-28 12:58:29 →
Size: 3586
Editor: EricPalmer
Comment:
← Revision 10 as of 2024-10-31 11:21:29 → ⇥
Size: 3641
Editor: EricPalmer
Comment:
Deletions are marked like this. Additions are marked like this.
Line 3: Line 3:
----
Line 9: Line 10:
----
 * INTEGER NTMP
  NTMP is a global compile parameter that defines the maximum size of maplets (it allocates memory)
   Standard value is 128
Line 10: Line 15:
* INTEGER NTMP
  NTMP is a global compile parameter that defines the maximum size of maplets (it allocates memory)
  Standard value is 128		  * INTEGER QSZ
  QSZ is the half size of maplet. width = QSZ * 2 + 1
   It is derived from [[READ_MAP]]
Line 14: Line 19:
* INTEGER QSZ
  QSZ is the half size of maplet. width = QSZ * 2 + 1
  It is derived from [[READ_MAP]]		  * DOUBLE PRECISION SCALE
  The scale (both horizontally x and y, as well as verticlally). The values in the [[MAPLET]] should be multiplied by this value to get its actual value. Units are in km.
   It is derived from [[READ_MAP]]
Line 18: Line 23:
* DOUBLE PRECISION SCALE
  The scale (both horizontally x and y, as well as verticlally). The values in the [[MAPLET]] should be multiplied by this value to get its actual value. Units are in km.
  It is derived from [[READ_MAP]]

* CHARACTER*12 PICNM		  * CHARACTER*12 PICNM
Line 25: Line 26:
* DOUBLE PRECISION IPL(2)  * DOUBLE PRECISION IPL(2)
Line 27: Line 28:
  Comes from [[V2IMGPLV]]   Comes from [[V2IMGPLV]]
Line 31: Line 32:
* DOUBLE PRECISION V(3)  * DOUBLE PRECISION V(3)
Line 33: Line 34:
  It is derived from [[READ_MAP]]   It is derived from [[READ_MAP]]
Line 35: Line 36:
* DOUBLE PRECISION UX(3)  * DOUBLE PRECISION UX(3)
 * DOUBLE PRECISION UY(3)
 * DOUBLE PRECISION UZ(3)
  Normal plane of the maplet in X, Y and Z space
   It is derived from [[READ_MAP]]
Line 37: Line 42:
* DOUBLE PRECISION UY(3)

* DOUBLE PRECISION UZ(3)
  Normal plane of the maplet in X, Y and Z space
  It is derived from [[READ_MAP]]

* LOGICAL HUSE(-NTMP:NTMP,-NTMP:NTMP)		  * LOGICAL HUSE(-NTMP:NTMP,-NTMP:NTMP)
Line 45: Line 44:
  Defined in [[HGT2SLP]]   Defined in [[HGT2SLP]]
Line 48: Line 47:
* REAL*4 HT(-NTMP:NTMP,-NTMP:NTMP)  * REAL*4 HT(-NTMP:NTMP,-NTMP:NTMP)
Line 52: Line 51:
* DOUBLE PRECISION CP(3)  * DOUBLE PRECISION CP(3)
Line 54: Line 53:
  Calculated in this subroutine   Calculated in this subroutine
Line 56: Line 55:
* DOUBLE PRECISION SP(3)  * DOUBLE PRECISION SP(3)
Line 58: Line 57:
  Calculated in this subroutine   Calculated in this subroutine
Line 61: Line 60:
* LOGICAL EUSE  * LOGICAL EUSE
Line 64: Line 63:
  Uses the ETOL variable   Uses the ETOL variable
Line 66: Line 65:
* DOUBLE PRECISION ETOL  * DOUBLE PRECISION ETOL
Line 69: Line 68:
  Typically set to .025 (standard within the scripts for LITHOSP)   Typically set to .025 (standard within the scripts for LITHOSP)
Line 72: Line 71:
* DOUBLE PRECISION FINV  * DOUBLE PRECISION FINV
Line 75: Line 74:
* DOUBLE PRECISION DPIC_DLOC(2,3)  * DOUBLE PRECISION DPIC_DLOC(2,3)
Line 77: Line 76:
   Changed within [[LOC2PIX]] that is typically called from this subroutine.   Changed within [[LOC2PIX]] that is typically called from this subroutine.
Line 81: Line 80:
* DOUBLE PRECISION DIDH  * DOUBLE PRECISION DIDH
 * DOUBLE PRECISION DJDH
  Rate of change for the position I and J in the picture
   Calcualted in this routine and reported back to who called it
Line 83: Line 85:
* DOUBLE PRECISION DJDH
  Rate of change for the position I and J in the picture
  Calcualted in this routine and reported back to who called it		  * INTEGER Z0(2)
  Center point to extract the data from the image.
   Defined in [[lithos]] and derived from IPL. Used in [[PICINPT]] and [[BRT1]]
Line 87: Line 89:
* INTEGER Z0(2)
  Center point to extract the data from the image.
  Defined in [[lithos]] and derived from IPL. Used in [[PICINPT]] and [[BRT1]]		  * INTEGER KK, ESZ
  * KK is a spacing parameter, which I think is based upon the maplet.
  * ESZ is the size of the array.
   Defined in [[lithos]]. Used in [[PICINPT]] and [[BRT1]]
Line 91: Line 94:
* INTEGER KK, ESZ
  KK is a spacing parameter, which I think is based upon the maplet.
  ESZ is the size of the array.
  Defined in [[lithos]]. Used in [[PICINPT]] and [[BRT1]]		  * INTEGER DN(-1300:1300,-1300:1300)
  The workspace of displayed DN for both the image and the maplet (template)
   Set in this subroutine.
Line 96: Line 98:
*      INTEGER DN(-1300:1300,-1300:1300)
  The workspace of displayed DN for both the image and the maplet (template)
  Set in this subroutine.

* REAL*4 DNX(-NTMP:NTMP,-NTMP:NTMP)
  The filled out data that we really want. Set in this routine.		  * REAL*4 DNX(-NTMP:NTMP,-NTMP:NTMP)
   '''The filled out data that we really want.''' Set in this routine.

Describe EXTRACT_DATA_PIC here. 

      SUBROUTINE EXTRACT_DATA_PIC(NTMP,QSZ,SCALE,PICNM,IPL,
     .                     V,UX,UY,UZ,HUSE,HT,CP,SP,EUSE,ETOL,
     .                     FINV,DPIC_DLOC,DIDH,DJDH,Z0,KK,ESZ,DN, DNX)

  • INTEGER NTMP
    • NTMP is a global compile parameter that defines the maximum size of maplets (it allocates memory)
      • Standard value is 128
  • INTEGER QSZ
    • QSZ is the half size of maplet. width = QSZ * 2 + 1
  • DOUBLE PRECISION SCALE

    • The scale (both horizontally x and y, as well as verticlally). The values in the MAPLET should be multiplied by this value to get its actual value. Units are in km.

  • CHARACTER*12 PICNM
    • The name of the image that you want to use. The file is read in (again) in this program)
  • DOUBLE PRECISION IPL(2)
    • IPL = Image Pixel Line. Scaling factor to translate from physical to image frame.
  • DOUBLE PRECISION V(3)
    • Central vector of the maplet
  • DOUBLE PRECISION UX(3)
  • DOUBLE PRECISION UY(3)
  • DOUBLE PRECISION UZ(3)
    • Normal plane of the maplet in X, Y and Z space

  • LOGICAL HUSE(-NTMP:NTMP,-NTMP:NTMP)

    • Tracks each pixel of the maplet to define if its value should be considered

  • REAL*4 HT(-NTMP:NTMP,-NTMP:NTMP)

    • HT - Height. The height of the template at that location of the maplet. Used for trying to identifing shadow.
  • DOUBLE PRECISION CP(3)
    • Camera vector
      • Calculated in this subroutine
  • DOUBLE PRECISION SP(3)
    • Sun vector
      • Calculated in this subroutine
  • LOGICAL EUSE
    • "Expanded T1" Flag for when to use the expanded threshold for DN.

      Is a flag set within lithos or lithosp using the 'x' option

      • Uses the ETOL variable
  • DOUBLE PRECISION ETOL
    • Expanded T1 flag variable.

      Is a flag set within lithos or lithosp using the 'x' option

      • Typically set to .025 (standard within the scripts for LITHOSP)
  • DOUBLE PRECISION FINV
    • Reports the fraction of invisible pixels within the maplet/image intersection
  • DOUBLE PRECISION DPIC_DLOC(2,3)
    • This subroutine computes the change in the image-space location (p,l) of a point (x,y,z) on a maplet as the point is moved in the maplet frame.

      • Changed within LOC2PIX that is typically called from this subroutine.

  • DOUBLE PRECISION DIDH
  • DOUBLE PRECISION DJDH
    • Rate of change for the position I and J in the picture
      • Calcualted in this routine and reported back to who called it
  • INTEGER Z0(2)
    • Center point to extract the data from the image.
  • INTEGER KK, ESZ
    • KK is a spacing parameter, which I think is based upon the maplet.
    • ESZ is the size of the array.
  • INTEGER DN(-1300:1300,-1300:1300)
    • The workspace of displayed DN for both the image and the maplet (template)
      • Set in this subroutine.

  • REAL*4 DNX(-NTMP:NTMP,-NTMP:NTMP)

    • The filled out data that we really want. Set in this routine.

categorySubroutine

EXTRACT_DATA_PIC (last edited 2024-10-31 11:21:44 by EricPalmer)