Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1147 1821 (N57E16 L= 36)  1688    235        3 (1049/1771-1235/1852)  49.4 E23-E07 N48-N63
      2 1779 1097 (N19W51 L=103)  1989    340        5 (1726/ 984-1835/1228)  73.6 W46-W57 N10-N27
      3 1703 1156 (N20W43 L= 95)   517     98        1 (1664/1121-1738/1194) -15.9 W39-W46 N18-N22
      4  231 1138 (S10E59 L=353)   933    152        1 ( 186/1068- 280/1192)  53.1 E63-E53 S14-S06
      5 1653 1115 (N17W39 L= 91)   908    203        1 (1591/1021-1704/1172) -28.6 W33-W44 N12-N20
      6 1574 1098 (N14W33 L= 85)   276     64        1 (1554/1068-1590/1134)  78.7 W31-W34 N11-N16
      7 1211  621 (S22W21 L= 72)   301    130        0 (1154/ 591-1266/ 657)  45.0 W17-W23 S24-S18
      8 1701  509 (S17W65 L=117)   960    126        2 (1663/ 451-1738/ 562) -69.4 W61-W68 S21-S13
      9 1476  446 (S27W47 L= 99)  1030    159        2 (1402/ 423-1549/ 478)  33.7 W41-W52 S29-S22
     10 1283  379 (S36W37 L= 88)   302     56        0 (1258/ 365-1315/ 399)  14.0 W34-W39 S36-S34
     11 1154  238 (S50W39 L= 91)   477     56        0 (1122/ 212-1176/ 259)  -6.3 W33-W42 S52-S48
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y