Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  722 1875 (N50E50 L= 76)   529     60        0 ( 698/1846- 749/1899) -54.5 E53-E47 N46-N53
      2 1293  850 (S12W19 L=145)   246     84        0 (1263/ 814-1334/ 876)  53.1 W17-W20 S14-S10
      3 1486  782 (S13W33 L=158)   349    131        0 (1450/ 719-1520/ 824)  74.1 W31-W33 S17-S10
      4 1720  647 (S17W56 L=182)   388    106        2 (1685/ 624-1755/ 670)   0.0 W52-W60 S18-S15
      5 1562  498 (S29W49 L=175)  1543    359        2 (1402/ 442-1703/ 588)  24.4 W37-W60 S34-S22
      6 1439  522 (S30W38 L=164)   346    109        0 (1407/ 484-1476/ 560)  80.5 W36-W39 S32-S26
      7 1645  417 (S32W64 L=190)   708     93        1 (1591/ 395-1692/ 442)  29.1 W58-W70 S34-S29
      8 1328  171 (S57W64 L=190)   645     53        1 (1305/ 143-1348/ 202)  -8.1 W56-W77 S58-S53
  

    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