Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  588 1706 (N58E36 L= 63)  1752    308        4 ( 422/1679- 721/1735)   0.0 E67-E17 N55-N60
      2 1292 1482 (N30W26 L=125)   235     43        1 (1269/1462-1318/1501) -33.7 W24-W27 N28-N31
      3  500 1364 (N33E34 L= 65)  1239    333        2 ( 411/1224- 571/1468) -46.7 E45-E27 N23-N40
      4  381 1271 (N28E44 L= 55)   376    104        0 ( 356/1231- 418/1323) -49.4 E47-E40 N24-N31
      5 1038 1163 (N14W03 L=102)   332    135        0 (1000/1110-1062/1228)  56.3 W00-W04 N11-N17
      6 1897  788 (S24W70 L=168)  1347    172        2 (1854/ 742-1921/ 878)  69.4 W62-W76 S26-S18
      7 1555  824 (S14W32 L=131)   377     67        1 (1529/ 791-1579/ 860)  76.0 W29-W33 S15-S11
      8 1418  536 (S31W21 L=120)  2530    569        4 (1176/ 395-1629/ 645)  17.0 W01-W37 S37-S26
      9  550  604 (S15E37 L= 62)   683    155        0 ( 478/ 572- 628/ 639)   8.7 E43-E30 S16-S12
     10  832  275 (S42E29 L= 70)   362     52        1 ( 802/ 255- 860/ 302) -21.8 E31-E26 S44-S39
     11  691  206 (S48E51 L= 48)   916     74        0 ( 655/ 186- 735/ 224)   7.6 E59-E44 S49-S45
  

    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