Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1145 1821 (N57E16 L= 36)  1732    228        3 (1047/1772-1234/1852)  53.7 E24-E07 N48-N63
      2 1760 1116 (N19W49 L=101)  2248    349        9 (1664/1023-1833/1195) -28.3 W39-W57 N13-N24
      3  226 1136 (S10E59 L=353)   971    148        3 ( 181/1068- 276/1194)  56.3 E64-E53 S15-S06
      4 1651 1114 (N17W39 L= 91)   978    208        1 (1581/1025-1702/1170) -26.6 W32-W44 N12-N20
      5 1570 1098 (N14W33 L= 85)   322     68        2 (1550/1068-1587/1133)  78.7 W31-W34 N11-N16
      6 1215  628 (S21W21 L= 73)   335    135        0 (1155/ 594-1267/ 665)  54.5 W17-W22 S24-S17
      7 1699  509 (S18W65 L=117)   888    122        1 (1667/ 453-1735/ 559) -69.4 W61-W68 S21-S13
      8 1473  446 (S27W47 L= 99)  1046    155        2 (1402/ 421-1546/ 480)  37.9 W41-W51 S29-S22
      9 1260  382 (S36W35 L= 87)   495    121        1 (1202/ 367-1313/ 399)  12.5 W29-W38 S37-S33
     10 1153  240 (S50W39 L= 91)   527     71        1 (1113/ 213-1183/ 264)   0.0 W32-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