Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  759 1569 (N36E13 L= 60)   349     58        1 ( 723/1550- 788/1588) -21.8 E16-E11 N34-N37
      2  263 1357 (N27E56 L= 18)  2988    363        7 ( 228/1202- 325/1543) -42.7 E77-E46 N17-N41
      3 1368 1456 (N22W27 L=101)   864    218        0 (1261/1429-1480/1478) -19.7 W20-W34 N19-N24
      4 1863 1276 (N07W66 L=140)  2855    220        5 (1793/1214-1939/1350) -19.3 W58-W78 N02-N11
      5  220 1023 (N07E57 L= 17)   361     30        0 ( 202/1003- 237/1048)  71.6 E58-E55 N05-N08
      6 1744  914 (S14W48 L=122)   330     22        1 (1721/ 891-1761/ 935)  18.4 W46-W50 S15-S12
      7  257  620 (S18E64 L=  9)   548     49        1 ( 234/ 588- 283/ 641)  23.2 E67-E59 S19-S16
      8  934  618 (S25E10 L= 64)   292     60        0 ( 904/ 601- 971/ 639) -26.6 E11-E07 S26-S23
      9 1326  446 (S41W17 L= 90)   667     99        1 (1276/ 427-1376/ 475)   7.1 W12-W21 S42-S38
     10 1218  414 (S42W08 L= 81)   455     74        1 (1194/ 390-1253/ 451) -59.0 W05-W10 S44-S39
     11 1131  381 (S44W00 L= 74)   392     49        0 (1106/ 362-1164/ 410) -14.0 E01-W03 S45-S41
  

    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