Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  651 1230 (N26E17 L=147)   762    203        3 ( 625/1140- 681/1332)  76.0 E19-E14 N20-N32
      2  807 1228 (N22E07 L=156)   252    100        0 ( 767/1202- 845/1243)   0.0 E10-E04 N20-N22
      3  253 1060 (N25E53 L=111)   353     85        2 ( 232/1015- 270/1098) -90.0 E55-E50 N22-N27
      4  169 1049 (N26E64 L=100)  1200    126        6 ( 129/1000- 217/1092)   8.1 E70-E56 N22-N28
      5 1558  842 (S20W27 L=191)   410    154        2 (1501/ 785-1620/ 886)  11.3 W22-W32 S21-S18
      6  316  505 (S09E66 L= 97)   511     68        2 ( 300/ 479- 345/ 543)  33.7 E70-E63 S10-S06
      7 1235  403 (S38E05 L=159)   543    153        1 (1167/ 368-1303/ 449)   5.2 E10-W01 S39-S36
      8  536  265 (S29E66 L= 98)  1220     90        3 ( 508/ 222- 557/ 311)  20.2 E78-E58 S32-S25
  

    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