Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  476 1670 (N35E55 L=360)  1967    425        5 ( 435/1590- 537/1746) -24.2 E65-E45 N30-N40
      2  906 1591 (N36E17 L= 38)   555    159        1 ( 844/1556- 964/1648) -20.0 E23-E12 N33-N39
      3 1414 1496 (N36W23 L= 77)  1213    283        5 (1353/1392-1474/1602) -60.3 W17-W26 N28-N43
      4  721 1380 (N19E24 L= 30)  1912    425        4 ( 593/1239- 813/1515) -33.7 E37-E15 N11-N26
      5 1106 1432 (N27W00 L= 55)  1037    218        1 (1035/1330-1192/1513) -39.8 E05-W07 N21-N31
      6 1703 1192 (N19W45 L=100)   855    200        0 (1663/1097-1727/1302) -71.6 W42-W47 N13-N25
      7  252 1228 (N04E59 L=356)   381     85        1 ( 241/1186- 268/1276) -78.7 E60-E57 N01-N06
      8 1526 1197 (N17W31 L= 85)   754    250        0 (1455/1136-1565/1261)  49.4 W25-W33 N13-N21
      9 1332 1134 (N11W18 L= 72)   206     78        0 (1292/1105-1371/1160) -21.8 W15-W20 N09-N12
     10 1889 1071 (N13W66 L=121)  2199    280        4 (1857/ 986-1918/1151)  70.0 W61-W70 N07-N18
     11 1676 1044 (N09W43 L= 98)   806    161        3 (1646/ 971-1709/1137) -73.3 W40-W46 N04-N14
     12 1101 1062 (N04W04 L= 59)   602    165        1 (1020/1039-1185/1089)  26.6 W00-W09 N01-N06
     13  466  941 (S11E35 L= 19)  1327    431        3 ( 310/ 903- 649/ 997)   0.0 E48-E22 S13-S07
     14 1188  831 (S09W12 L= 67)   239    118        0 (1160/ 768-1219/ 882)  76.0 W11-W13 S13-S05
     15  317  814 (S21E48 L=  7)   346     59        0 ( 300/ 786- 344/ 844) -45.0 E49-E45 S22-S19
     16  263  778 (S24E53 L=  1)   377     91        0 ( 244/ 752- 281/ 813) -90.0 E55-E51 S25-S21
     17 1736  706 (S10W55 L=110)   867    179        1 (1706/ 666-1786/ 766) -56.3 W51-W61 S12-S06
     18 1522  378 (S35W50 L=105)  1419    325        4 (1436/ 356-1647/ 405)  16.4 W42-W60 S37-S31
  

    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