Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1326 1680 (N35W38 L=348)   389     66        0 (1295/1654-1363/1698)   0.0 W35-W41 N33-N36
      2 1534 1629 (N27W52 L=  1)  2177    302        3 (1370/1568-1635/1666) -34.5 W39-W60 N21-N32
      3 1746 1518 (N16W65 L= 15)   485     51        0 (1726/1499-1768/1540) -71.6 W62-W68 N14-N17
      4 1444 1470 (N19W36 L=346)   287     55        0 (1415/1456-1481/1486) -26.6 W33-W38 N17-N20
      5  557 1282 (N25E25 L=285)  2403    550        3 ( 332/1218- 826/1343)  -9.0 E44-E06 N21-N29
      6 1852 1268 (S00W64 L= 14)   880    104        1 (1811/1235-1878/1317)   9.5 W58-W70 S02-N01
      8 1882  871 (S24W62 L= 12)  1464    194        1 (1853/ 765-1907/ 950)  84.8 W58-W65 S30-S19
      9 1298  773 (S18W12 L=321)   415    150        0 (1265/ 714-1321/ 855)  63.4 W08-W14 S21-S13
     10 1131  801 (S13W02 L=312)   221     71        0 (1111/ 768-1157/ 831)  36.9 W00-W04 S14-S11
     11  465  409 (S25E56 L=254)   464     65        0 ( 448/ 376- 481/ 445)  29.7 E60-E52 S26-S22
  

    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