Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1767 1369 (N18W56 L=276)  1433    310        5 (1670/1254-1853/1484) -47.1 W49-W62 N11-N25
      2  347 1271 (N14E46 L=174)   678    180        1 ( 279/1238- 421/1290)  -8.7 E52-E39 N11-N15
      3  224 1154 (N07E57 L=164)   389     59        2 ( 206/1122- 243/1191) -45.0 E58-E54 N05-N09
      4  169 1021 (S00E63 L=157)   731    100        0 ( 146/ 977- 189/1066)  38.7 E66-E60 S02-N02
      5  517  818 (S14E33 L=187)   869    264        1 ( 494/ 700- 547/ 940)  78.7 E36-E31 S21-S06
      6 1076  863 (S13W03 L=223)   220     82        1 (1037/ 836-1102/ 893) -36.9 W00-W04 S14-S10
      7 1049  793 (S17W01 L=222)   274    138        0 ( 984/ 777-1109/ 831) -26.6 E02-W04 S17-S14
      8 1707  562 (S32W55 L=276)   634    111        0 (1672/ 514-1751/ 610)  63.4 W52-W58 S34-S28
      9 1577  457 (S39W47 L=267)   903    137        3 (1504/ 435-1654/ 481)   3.8 W39-W55 S40-S37
  

    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