Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1199 1792 (N49W01 L= 10)   621    102        2 (1164/1755-1232/1826) -56.3 E02-W03 N45-N52
      2 1069 1806 (N48E11 L=358)   330     44        0 (1045/1783-1097/1825)  53.1 E12-E08 N45-N50
      3 1508 1590 (N38W27 L= 36)   511     90        0 (1483/1554-1542/1613)  50.2 W24-W30 N34-N40
      4  584 1528 (N20E36 L=333)   276     60        2 ( 559/1504- 612/1554)  18.4 E37-E34 N18-N21
      5  230 1469 (N14E67 L=302)   746    245        3 ( 200/1437- 261/1511) -33.7 E73-E62 N11-N17
      6  507 1431 (N13E39 L=330)   286     61        2 ( 482/1403- 529/1458) -45.0 E41-E37 N11-N14
      7  378 1405 (N10E49 L=320)   348     99        1 ( 343/1388- 427/1424)  15.9 E51-E44 N09-N11
      8  441 1270 (N02E40 L=329)   291     86        3 ( 412/1247- 468/1292)  14.0 E42-E38 N01-N03
      9 1891 1170 (N17W61 L= 70)  2226    393        6 (1844/1092-1930/1236) -37.9 W55-W66 N13-N21
     10 1391 1155 (N07W20 L= 29)   204     53        1 (1374/1121-1408/1179) -76.0 W19-W21 N04-N08
     11  359 1136 (S06E45 L=324)   361     48        1 ( 337/1107- 380/1166) -36.9 E46-E42 S07-S04
     12  748 1098 (S05E17 L=352)   370    133        0 ( 685/1044- 791/1149) -35.5 E21-E14 S07-S02
     13  718  942 (S14E18 L=351)   505    127        0 ( 670/ 901- 794/ 984)  35.5 E20-E13 S17-S11
     14  147  892 (S21E64 L=305)  1521    200        2 ( 111/ 845- 198/ 948)  29.1 E68-E58 S24-S18
     15 1344  881 (S10W21 L= 30)   440    134        1 (1294/ 827-1398/ 934)  66.8 W18-W23 S13-S06
     16 1852  874 (S01W61 L= 70)   555    138        0 (1811/ 816-1887/ 924)  49.4 W57-W64 S05-N02
     17 1028  818 (S18W03 L= 12)   206     71        1 ( 987/ 807-1070/ 832)   0.0 W00-W05 S18-S17
     18  500  703 (S31E32 L=337)   472     93        2 ( 464/ 674- 545/ 737) -31.0 E35-E28 S33-S29
     19 1802  577 (S19W66 L= 75)   973    152        5 (1767/ 534-1838/ 630)  56.3 W63-W70 S21-S15
  

    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