Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  894 1855 (N61W05 L=157)   683     86        1 ( 864/1817- 930/1883)  33.7 W00-W10 N57-N64
      2 1002 1541 (N32W06 L=158)   467    147        0 ( 948/1486-1054/1588) -53.1 W02-W08 N27-N35
      3 1876 1307 (N06W70 L=222)   602     60        1 (1860/1284-1897/1344)  45.0 W66-W72 N04-N08
      4 1765 1285 (N06W55 L=207)   486     69        0 (1741/1249-1783/1327)  33.7 W51-W57 N04-N08
      5  372 1244 (N21E42 L=110)  1186    196        5 ( 308/1188- 460/1291)  13.0 E47-E34 N18-N23
      6 1572  885 (S15W33 L=185)  1469    376        2 (1527/ 714-1605/1043)  81.5 W30-W35 S25-S05
      7 1873  906 (S18W62 L=214)  4272    367        7 (1781/ 765-1930/1030) -67.6 W51-W71 S26-S09
      8 1460  768 (S21W25 L=177)   220    100        0 (1429/ 730-1484/ 816)  45.0 W22-W26 S22-S18
      9 1707  601 (S35W48 L=200)   797    114        0 (1662/ 563-1765/ 643)  21.8 W43-W54 S37-S33
     10  575  575 (S22E37 L=115)   316     96        0 ( 558/ 526- 605/ 612)  39.8 E39-E33 S24-S19
     11  317  569 (S19E60 L= 92)   980    104        5 ( 280/ 537- 360/ 602) -38.7 E64-E54 S20-S16
  

    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