Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1317 1660 (N42W35 L=287)   326     62        1 (1284/1644-1348/1674)   0.0 W32-W38 N40-N43
      2 1537 1545 (N30W46 L=298)   524    111        2 (1496/1506-1566/1573) -68.2 W43-W48 N27-N32
      3 1844 1164 (N02W63 L=315)  4024    625       10 (1777/ 944-1881/1361)  80.5 W58-W71 S11-N13
      4 1010 1334 (N25W03 L=255)   214     50        1 ( 985/1310-1039/1358) -45.0 W01-W04 N23-N26
      5  541 1158 (N19E30 L=222)  1871    403        7 ( 371/1065- 691/1273) -32.0 E44-E20 N12-N27
      6 1326 1169 (N11W21 L=272)   365    120        0 (1268/1110-1363/1218)  40.6 W16-W23 N07-N13
      7  377 1025 (N12E44 L=208)   468    100        0 ( 336/ 982- 399/1056) -31.0 E47-E42 N09-N13
      8 1900  971 (S11W67 L=319)   686    134        2 (1878/ 928-1918/1013) -63.4 W63-W69 S14-S08
      9 1618  745 (S20W37 L=289)  2345    414        5 (1496/ 607-1701/ 900)  45.0 W27-W43 S27-S11
     10  921  777 (S08E09 L=243)   232     37        1 ( 900/ 755- 951/ 798)  18.4 E10-E07 S09-S06
     11  610  758 (S06E29 L=222)   283     86        2 ( 582/ 728- 641/ 794)  31.0 E31-E26 S08-S03
     12  986  693 (S14E06 L=246)   955    326        1 ( 859/ 613-1090/ 793)  28.1 E15-W00 S17-S08
     13  816  727 (S10E16 L=235)   355    105        2 ( 795/ 681- 836/ 788)  63.4 E18-E15 S12-S06
     14  691  691 (S11E25 L=227)   647    157        2 ( 631/ 632- 762/ 749) -71.6 E28-E20 S15-S06
     15 1378  630 (S24W18 L=270)   503    186        2 (1323/ 559-1406/ 686) -78.7 W14-W20 S28-S18
     16 1285  598 (S24W12 L=264)   236     52        2 (1272/ 568-1302/ 626) -90.0 W10-W13 S26-S22
     17  557  584 (S17E37 L=215)   534    128        1 ( 502/ 547- 600/ 619) -51.3 E40-E34 S19-S14
     18  544  498 (S23E41 L=211)   313     53        0 ( 518/ 481- 567/ 516)  21.8 E43-E38 S23-S21
  

    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