Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  466 1641 (N16E56 L=278)  1026    182        3 ( 407/1571- 518/1705)  90.0 E57-E54 N10-N20
      2 1137 1263 (N09E00 L=333)   483    131        0 (1086/1191-1179/1313) -33.7 E04-W04 N05-N10
      3  404 1200 (S12E43 L=291)   361     54        1 ( 376/1180- 436/1221)  36.9 E44-E40 S13-S09
      4 1531 1101 (N11W27 L=  0)   442    122        0 (1515/1036-1545/1167) -60.3 W25-W29 N07-N14
      5 1380  999 (N01W20 L=354)   689    249        0 (1299/ 882-1419/1100) -33.7 W13-W25 S04-N04
      6  680  818 (S27E15 L=318)   681    278        0 ( 556/ 791- 821/ 843)  11.9 E24-E05 S30-S25
      7  398  762 (S38E36 L=298)   302     91        0 ( 364/ 751- 443/ 780)   7.1 E39-E31 S38-S36
      8  875  556 (S38W05 L=338)   300     86        0 ( 836/ 531- 910/ 590)  71.6 W02-W06 S40-S34
      9 1542  390 (S25W56 L= 30)   486    113        0 (1512/ 375-1581/ 412)   7.1 W52-W60 S26-S24
  

    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