Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1410 1628 (N31W47 L=347)   583    120        0 (1362/1593-1459/1674) -74.1 W45-W48 N27-N35
      2 1663 1583 (N20W66 L=  6)   657     69        1 (1625/1567-1703/1598) -20.6 W61-W70 N18-N22
      3  985 1507 (N37W11 L=311)   263     66        1 ( 950/1481-1019/1524) -11.3 W08-W13 N35-N37
      4 1458 1433 (N19W39 L=339)   616    120        0 (1403/1408-1516/1454) -20.6 W35-W43 N17-N21
      5 1265 1169 (N10W17 L=318)   272    104        0 (1228/1122-1298/1210) -71.6 W16-W18 N06-N12
      6  502 1109 (N23E32 L=269)   805    310        1 ( 370/1072- 646/1138) -20.9 E42-E21 N18-N27
      7 1703 1072 (S08W44 L=344)   497     74        0 (1680/1036-1727/1121)  26.6 W41-W47 S09-S05
      8  150 1007 (N22E70 L=230)   928    128        1 ( 122/ 961- 180/1068) -33.7 E76-E65 N19-N25
      9 1314 1028 (N01W17 L=317)   224     68        1 (1293/ 990-1338/1061)  45.0 W14-W18 S01-N02
     10  344 1010 (N20E46 L=254)   675     99        1 ( 302/ 971- 377/1059)  33.7 E50-E43 N17-N22
     11 1836  935 (S20W55 L=355)   668    117        0 (1792/ 897-1859/ 975)  12.5 W49-W58 S22-S17
     12 1245  911 (S04W10 L=310)   301     70        2 (1213/ 884-1279/ 935) -45.0 W08-W11 S06-S02
     13  390  883 (N12E44 L=256)   346    119        0 ( 331/ 862- 438/ 898) -18.4 E49-E40 N09-N12
     14 1308  620 (S24W07 L=308)  1204    372        2 (1213/ 496-1419/ 768)  33.7 E01-W17 S29-S17
     15 1526  520 (S38W22 L=322)   380     81        0 (1512/ 479-1537/ 563) -81.9 W20-W23 S41-S34
     16  724  468 (S20E33 L=267)   332     99        0 ( 700/ 434- 756/ 512) -90.0 E34-E32 S22-S16
     17 1135  300 (S43E14 L=286)  2604    306        4 (1036/ 214-1213/ 424)  31.8 E26-E04 S48-S34
  

    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