Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1192 1469 (N24W08 L=218)   317     46        0 (1154/1448-1225/1485) -11.3 W05-W10 N22-N24
      2  290 1352 (N13E54 L=156)   380     65        1 ( 270/1325- 307/1385)  71.6 E56-E52 N11-N14
      3 1632 1154 (N08W38 L=248)   828    193        0 (1560/1078-1708/1240) -36.0 W32-W44 N03-N12
      4 1748 1215 (N13W48 L=258)   325     68        1 (1715/1191-1767/1240)   0.0 W45-W49 N11-N14
      5  826 1074 (S03E12 L=198)   526    136        1 ( 807/1001- 849/1141) -76.0 E13-E10 S07-N01
      6  810  979 (S08E13 L=197)   178     51        0 ( 780/ 961- 834/ 999)  33.7 E14-E11 S09-S07
      7 1379  919 (S08W22 L=232)   215     66        0 (1346/ 885-1411/ 953)  45.0 W20-W24 S10-S06
      8 1854  855 (S07W62 L=272)   690    100        1 (1822/ 818-1873/ 904)  63.4 W58-W64 S10-S04
      9 1262  862 (S13W16 L=225)   322    102        0 (1209/ 840-1309/ 876)  -9.5 W12-W18 S13-S11
     10 1552  730 (S18W37 L=247)  2002    356        1 (1389/ 622-1712/ 792)  28.6 W27-W49 S26-S13
     11  337  659 (S30E49 L=161)   664    134        0 ( 303/ 611- 368/ 711)  40.6 E53-E45 S32-S26
     12 1170  471 (S38W15 L=225)   225     76        0 (1127/ 450-1216/ 490)  33.7 W12-W18 S40-S36
  

    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