Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1670 1331 (N17W50 L=198)   430    122        0 (1617/1314-1724/1360) -18.4 W45-W55 N15-N19
      2 1725 1060 (N00W48 L=197)   314     44        1 (1702/1036-1744/1093)  63.4 W46-W50 S01-N02
      3  908 1030 (N06E07 L=142)   401    130        1 ( 883/ 966- 931/1089) -82.9 E08-E05 N02-N10
      4 1195  944 (S01W10 L=158)   229     97        0 (1152/ 911-1229/ 962) -26.6 W07-W11 S03-N00
      5  855  840 (S05E12 L=137)   352    128        0 ( 828/ 802- 911/ 885) -45.0 E13-E08 S07-S01
      6 1583  699 (S21W36 L=185)  1669    283        3 (1478/ 612-1666/ 800)  39.8 W27-W42 S25-S14
      7  923  715 (S13E09 L=140)   306    117        1 ( 868/ 696- 958/ 751)   9.5 E12-E06 S14-S11
      8 1385  489 (S33W22 L=170)   487     67        0 (1350/ 463-1416/ 522)  31.0 W19-W24 S35-S31
      9 1203  293 (S48W08 L=156)   313     42        0 (1185/ 276-1227/ 315) -26.6 W06-W10 S49-S45
  

    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