Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1406 1784 (N38W62 L=299)   417     66        0 (1390/1753-1435/1804)   5.7 W56-W66 N37-N39
      2 1719 1488 (N12W62 L=300)  1503    196        4 (1616/1450-1805/1538) -34.5 W54-W71 N07-N18
      3 1030 1360 (N26W09 L=247)   659    202        2 ( 938/1315-1109/1414) -52.1 W04-W13 N21-N31
      4 1069 1264 (N19W09 L=247)   415    171        0 ( 986/1235-1145/1283) -11.3 W03-W13 N17-N20
      5 1619 1148 (S02W39 L=276)   751    155        0 (1593/1068-1645/1209)  45.0 W35-W41 S06-N01
      6 1083 1077 (N09W05 L=242)   283    141        0 (1048/1017-1117/1141)  35.5 W01-W08 N06-N11
      7  329 1069 (N25E46 L=192)   307     50        1 ( 298/1044- 351/1094) -36.9 E48-E43 N22-N26
      8  267  910 (N16E54 L=184)   573    108        0 ( 232/ 862- 298/ 946) -51.3 E57-E51 N13-N18
      9 1868  876 (S26W57 L=294)   827    112        3 (1822/ 848-1901/ 910)   7.1 W50-W62 S27-S23
     10  200  854 (N14E63 L=175)   703     98        3 ( 176/ 813- 221/ 902) -68.2 E66-E60 N11-N16
     11  471  851 (N09E38 L=200)   412     64        0 ( 444/ 812- 499/ 886) -68.2 E39-E36 N06-N11
     12  249  786 (N09E59 L=179)   577     65        1 ( 218/ 759- 292/ 807)   0.0 E62-E54 N07-N10
     13  444  541 (S08E50 L=187)   650    163        0 ( 395/ 474- 471/ 596) -90.0 E53-E48 S12-S04
     14 1402  401 (S42W08 L=246)  3367    666        1 (1140/ 244-1650/ 575)   9.3 E18-W31 S48-S36
     15 1197  522 (S27E03 L=235)   365     70        0 (1178/ 485-1219/ 569)  80.5 E04-E02 S29-S23
     16  654  478 (S16E38 L=200)   807    187        2 ( 612/ 410- 711/ 546)  84.3 E40-E34 S21-S11
     17  789  431 (S22E31 L=207)   510    110        0 ( 740/ 403- 820/ 473)   8.1 E34-E27 S23-S19
     18  380  365 (S19E72 L=166)   949     50        1 ( 359/ 336- 399/ 392)  16.4 E83-E66 S22-S17
     19  766  227 (S37E47 L=191)   459     67        0 ( 749/ 195- 782/ 265)  33.7 E51-E42 S39-S33
  

    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