Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  615 1469 (N25E34 L= 29)  1580    401        2 ( 523/1391- 678/1619) -79.4 E40-E29 N19-N35
      2 1501 1411 (N30W31 L= 94)  2412    577        2 (1275/1231-1641/1530) -26.6 W13-W41 N19-N37
      3 1716 1232 (N20W48 L=112)   709    190        1 (1667/1166-1777/1291) -45.0 W44-W54 N16-N23
      4 1384 1190 (N15W22 L= 85)   379    152        1 (1353/1136-1425/1268) -58.0 W19-W24 N11-N19
      5 1561 1200 (N17W34 L= 98)   343    126        0 (1535/1141-1575/1252) -74.1 W32-W35 N13-N20
      6  146  963 (S11E70 L=353)   675     53        0 ( 120/ 946- 178/ 983)   0.0 E76-E64 S12-S10
      7  468  771 (S20E36 L= 27)   658    185        0 ( 447/ 694- 504/ 840)  73.3 E37-E32 S24-S14
      8 1497  669 (S16W36 L= 99)   601    140        1 (1455/ 618-1548/ 711)  56.3 W32-W38 S19-S13
      9 1743  569 (S21W64 L=127)  1343    178        0 (1683/ 517-1799/ 612)  36.9 W58-W72 S24-S17
     10 1277  494 (S30W24 L= 87)   285     95        1 (1227/ 474-1324/ 513)  15.9 W20-W27 S31-S28
  

    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