Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  689 1443 (N08E31 L=284)  1051    313        1 ( 653/1305- 722/1587) -62.1 E36-E26 N00-N17
      2 1794 1464 (N47W46 L=  2)  1029    193        0 (1761/1402-1810/1514)  78.7 W42-W49 N41-N51
      3 1428 1464 (N31W12 L=327)   260     47        0 (1399/1450-1459/1478)  36.9 W09-W13 N29-N32
      4 1672 1337 (N31W32 L=348)   496    105        0 (1632/1292-1702/1378) -23.2 W28-W35 N29-N33
      5  869 1307 (N05E16 L=299)   244     74        0 ( 829/1286- 911/1324)  36.9 E18-E13 N02-N05
      6  385 1157 (S15E43 L=272)  1914    498        3 ( 199/1113- 533/1185)  16.5 E59-E32 S19-S11
      7 1868  754 (N06W68 L= 24)   848    135        2 (1840/ 733-1902/ 780)   5.2 W63-W74 N04-N07
      8 1205  508 (S30W28 L=344)   695    121        2 (1144/ 475-1252/ 559)  76.0 W25-W30 S33-S25
      9 1704  506 (S13W64 L= 20)  1362    182        3 (1647/ 466-1754/ 537)  26.6 W57-W69 S15-S09
     10 1075  475 (S36W22 L=337)   264     28        0 (1057/ 452-1095/ 500) -33.7 W19-W23 S37-S34
     11 1510  355 (S27W57 L= 13)  1943    386        1 (1425/ 294-1589/ 408)   2.7 W47-W68 S29-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