Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1074 1853 (N58W44 L= 94)   449     69        0 (1041/1826-1107/1875) -14.0 W38-W50 N56-N59
      2 1585 1581 (N24W56 L=106)   699    111        0 (1539/1536-1631/1613) -63.4 W53-W58 N19-N26
      3 1256 1549 (N32W30 L= 80)   423     83        2 (1232/1513-1294/1582) -90.0 W28-W31 N28-N34
      4 1702 1385 (N10W55 L=105)  2720    462        5 (1593/1324-1841/1485) -39.8 W45-W69 N02-N17
      5 1401 1445 (N22W35 L= 85)   250     17        0 (1380/1424-1423/1463) -71.6 W33-W36 N20-N23
      6  532 1345 (N37E27 L= 23)   405     91        1 ( 493/1312- 567/1374) -39.8 E30-E24 N33-N39
      7 1649 1179 (N00W42 L= 92)   773    214        0 (1624/1073-1668/1280)  74.7 W39-W44 S05-N06
      8 1127 1192 (N14W10 L= 60)   664    296        2 (1061/1108-1186/1259)  26.6 W04-W14 N09-N17
      9  469 1078 (N21E34 L= 16)   519    149        1 ( 403/1041- 524/1117) -29.1 E39-E30 N19-N24
     10  752 1047 (N14E15 L= 35)   475    124        0 ( 699/1024- 815/1076) -39.8 E18-E11 N11-N16
     11 1873  983 (S18W59 L=109)  1261    155        1 (1834/ 909-1911/1048) -59.0 W55-W65 S23-S13
     12 1828  936 (S19W53 L=103)   909    174        1 (1791/ 866-1866/1020) -61.4 W49-W57 S24-S13
     13  677  932 (N09E22 L= 28)   302     41        1 ( 658/ 901- 697/ 962)  56.3 E23-E20 N07-N10
     14 1408  764 (S17W17 L= 67)  1376    432        3 (1256/ 685-1535/ 903)  14.0 W06-W27 S19-S11
     15  354  793 (N06E49 L=  1)   690    172        0 ( 328/ 722- 374/ 877)  83.7 E51-E46 N02-N11
     16 1592  519 (S39W28 L= 78)  4830    628        5 (1345/ 359-1831/ 675)   7.0 W05-W57 S45-S33
     17  696  589 (S11E30 L= 20)   320     67        1 ( 663/ 573- 736/ 613) -26.6 E31-E27 S12-S09
     18  318  505 (S11E65 L=345)  1425    207        1 ( 294/ 414- 356/ 604)  52.4 E72-E59 S17-S04
  

    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