Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  625 1679 (N58E25 L=  7)  2312    319       11 ( 515/1599- 770/1746)   9.5 E40-E07 N52-N62
      2  998 1607 (N44W11 L= 44)   729    184        1 ( 912/1580-1076/1634) -30.3 W04-W16 N40-N47
      3 1221 1587 (N38W27 L= 60)   675    146        2 (1162/1564-1300/1611) -10.3 W22-W34 N36-N39
      4  442 1165 (N23E38 L=355)   536    146        4 ( 422/1095- 468/1230) -82.9 E39-E35 N19-N27
      5  970 1184 (N17E01 L= 32)   346    113        1 ( 931/1132-1000/1221) -71.6 E02-W00 N13-N19
      6 1372  655 (S22W16 L= 49)  1443    323        7 (1304/ 548-1461/ 784) -77.5 W11-W21 S31-S13
      7 1786  713 (S29W53 L= 85)   563     90        1 (1747/ 690-1817/ 742)   7.1 W48-W56 S29-S26
      8 1580  639 (S28W32 L= 64)   273     67        0 (1557/ 611-1603/ 667)  45.0 W29-W33 S29-S26
  

    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