Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  845 1674 (N33E24 L=137)   310     36        1 ( 821/1655- 870/1690) -26.6 E25-E21 N30-N33
      2 1666 1303 (N21W39 L=199)  1615    289        3 (1591/1186-1740/1407) -37.6 W32-W45 N15-N27
      3  888 1289 (N08E12 L=149)   272     98        1 ( 846/1252- 924/1333) -31.0 E14-E09 N05-N09
      4 1102 1135 (N01W03 L=163)   369    112        1 (1060/1097-1154/1179) -26.6 W00-W06 N00-N03
      5 1118 1015 (S05W06 L=166)   333    198        0 (1015/ 986-1197/1037)   0.0 W00-W10 S06-S04
      6 1194  916 (S10W12 L=172)   241    132        0 (1162/ 853-1232/ 965)  76.0 W10-W13 S14-S06
      7  206  887 (S23E60 L=100)   635     72        2 ( 183/ 847- 228/ 924) -59.0 E62-E57 S25-S20
      8 1839  752 (S07W64 L=225)   644     68        2 (1813/ 724-1864/ 781) -76.0 W61-W68 S09-S05
      9 1218  689 (S23W17 L=178)   281     78        0 (1173/ 673-1259/ 705)   0.0 W14-W20 S23-S21
     10 1728  613 (S18W59 L=219)   832    146        1 (1654/ 584-1799/ 638)  19.7 W52-W66 S20-S15
  

    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