Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1520 1635 (N43W44 L=130)  2044    320        4 (1455/1569-1565/1699)  29.4 W35-W51 N38-N48
      2 1077 1505 (N35W04 L= 90)   571    153        1 (1054/1436-1099/1567) -90.0 W02-W05 N29-N39
      3  735 1518 (N35E22 L= 64)   370     47        2 ( 714/1489- 756/1543)  76.0 E24-E20 N33-N37
      4 1684 1320 (N21W48 L=134)   409    100        0 (1643/1279-1720/1349) -26.6 W44-W50 N18-N22
      5  723 1310 (N21E20 L= 66)   310     61        0 ( 701/1268- 741/1346) -68.2 E21-E18 N18-N23
      6 1060 1193 (N14W02 L= 88)   296    170        1 (1044/1110-1082/1274)  79.7 W01-W03 N08-N19
      7 1678 1234 (N15W46 L=132)   786    182        0 (1618/1194-1743/1263)   9.5 W40-W52 N13-N17
      8 1274 1191 (N14W16 L=102)   680    122        3 (1231/1141-1314/1259) -58.0 W13-W18 N10-N18
      9 1706 1111 (N08W47 L=134)   995    273        1 (1618/1078-1797/1176) -20.6 W40-W56 N05-N12
     10  773 1129 (N10E16 L= 70)   304    114        1 ( 719/1112- 827/1148)  -8.1 E19-E12 N08-N11
     11  656  981 (N01E23 L= 63)   435    108        1 ( 602/ 954- 713/1013)  15.9 E26-E19 S01-N02
     12 1700  659 (S21W50 L=136)  1741    160        4 (1596/ 543-1765/ 822)  62.1 W42-W54 S28-S10
     13  213  789 (S13E62 L= 24)   589    120        0 ( 178/ 767- 249/ 816) -14.0 E66-E57 S14-S11
     14  628  585 (S25E28 L= 59)   461    113        0 ( 574/ 565- 680/ 606)  14.0 E31-E23 S26-S23
     15 1267  429 (S36W19 L=105)   383    111        1 (1244/ 392-1300/ 471)  71.6 W17-W20 S39-S32
     16 1164  346 (S43W12 L= 98)  1130    204        4 (1123/ 303-1202/ 414)  84.3 W08-W15 S47-S37
     17  858  363 (S42E14 L= 72)   498     71        1 ( 825/ 340- 890/ 399)  11.3 E16-E11 S43-S38
  

    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