Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  576 1699 (N27E52 L=246)  1334    406        3 ( 526/1607- 617/1789) -30.1 E62-E43 N21-N32
      2  802 1743 (N35E38 L=259)   565     94        1 ( 762/1711- 851/1769)   0.0 E42-E33 N33-N37
      3  487 1503 (N12E48 L=250)   513     84        1 ( 453/1462- 525/1546)  71.6 E49-E46 N09-N15
      4  261 1482 (N07E71 L=227)   732    170        1 ( 248/1436- 272/1520) -33.7 E76-E66 N03-N09
      5  514 1419 (N08E42 L=255)   535     95        2 ( 478/1387- 549/1446)   0.0 E45-E39 N05-N09
      6  214 1163 (S13E60 L=237)   688     73        2 ( 197/1128- 232/1207) -78.7 E62-E57 S15-S10
      7  167 1124 (S16E65 L=233)   696     99        0 ( 151/1066- 184/1161) -80.5 E66-E62 S19-S13
      8  615  999 (S14E24 L=274)   405    101        1 ( 587/ 943- 654/1033)  68.2 E25-E21 S17-S12
      9  592  770 (S28E21 L=276)   309     65        0 ( 570/ 736- 621/ 795)  59.0 E22-E19 S31-S26
     10 1183  667 (S20W19 L=316)   278     63        0 (1144/ 648-1219/ 685)  36.9 W16-W20 S21-S18
     11  444  631 (S42E32 L=265)   710    140        3 ( 398/ 597- 496/ 669)  29.1 E36-E27 S45-S38
     12 1544  374 (S27W57 L=355)   520    127        1 (1509/ 357-1573/ 395)   0.0 W53-W61 S27-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