Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  549 1779 (N26E60 L= 90)   641    156        2 ( 525/1733- 573/1816) -54.5 E62-E57 N22-N29
      2  864 1689 (N28E31 L=118)   422    135        0 ( 803/1655- 921/1708)   0.0 E35-E26 N26-N29
      3  990 1504 (N20E16 L=134)   210     93        0 ( 972/1462-1015/1547) -38.7 E18-E13 N17-N21
      4 1676 1095 (N17W37 L=187)   422    126        0 (1635/1064-1717/1120)  18.4 W34-W40 N14-N18
      5 1061  487 (S35W20 L=170)   227     48        1 (1039/ 467-1089/ 513) -11.3 W18-W23 S36-S33
      6 1271  446 (S31W35 L=185)   237     80        1 (1252/ 424-1296/ 479) -21.8 W32-W37 S31-S29
      7 1634  410 (S20W63 L=213)   652    207        2 (1600/ 374-1666/ 439)   0.0 W58-W69 S20-S18
  

    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