Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1395 1739 (N53W42 L=207)   541     62        0 (1360/1715-1426/1764)  13.0 W35-W48 N50-N55
      2 1039 1735 (N54W02 L=167)   357     46        0 (1017/1710-1064/1757)  45.0 W00-W04 N51-N55
      3  935 1731 (N53E09 L=156)   449     42        1 ( 901/1714- 966/1746)   8.1 E12-E05 N51-N54
      4  811 1551 (N38E17 L=148)   233     49        0 ( 785/1535- 841/1565) -11.3 E19-E14 N36-N39
      5  969 1491 (N34E04 L=161)   258     31        0 ( 948/1471- 996/1506)  14.0 E05-E01 N32-N34
      6  242 1462 (N29E75 L= 90)   819     30        1 ( 223/1442- 261/1484) -33.7 E82-E69 N28-N30
      7 1196 1456 (N31W13 L=177)   218     30        0 (1169/1440-1219/1469) -14.0 W10-W14 N30-N32
      8  366 1348 (N23E50 L=115)   397     75        0 ( 350/1307- 382/1386) -51.3 E52-E47 N20-N25
      9  311 1277 (N18E54 L=111)   329     57        0 ( 295/1248- 325/1309)  76.0 E54-E52 N16-N20
     10  488 1268 (N18E37 L=128)   479    146        0 ( 430/1241- 552/1290)   0.0 E41-E32 N16-N19
     11  215 1065 (N05E61 L=104)   517     55        1 ( 190/1037- 245/1095)  33.7 E64-E57 N02-N06
     12  178  970 (S02E66 L= 99)   912     82        2 ( 140/ 946- 213/ 991)   9.5 E72-E60 S02-N00
     13  521  568 (S26E37 L=128)  1850    429        2 ( 401/ 415- 657/ 714)  33.3 E59-E24 S38-S15
     14 1396  683 (S18W25 L=190)   266    129        0 (1325/ 667-1448/ 699)   6.3 W19-W28 S19-S17
     15  801  487 (S32E17 L=148)  3448    596        9 ( 602/ 347-1002/ 677)  34.8 E37-E01 S43-S18
     16 1609  511 (S31W47 L=212)   579     64        0 (1576/ 485-1638/ 544)  53.1 W44-W49 S33-S28
     17 1528  320 (S47W53 L=218)  1517    108        2 (1496/ 285-1569/ 365) -19.8 W44-W69 S51-S42
  

    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