Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1129 1684 (N46E04 L=119)  1481    248        2 (1086/1595-1173/1774) -68.2 E08-W00 N38-N54
      2  607 1725 (N39E48 L= 76)   443    105        0 ( 582/1703- 639/1748)  -7.1 E51-E43 N37-N40
      3  516 1547 (N25E45 L= 78)  2174    325        4 ( 393/1473- 617/1667)  47.3 E54-E40 N19-N34
      4 1470 1581 (N44W26 L=150)   658    137        1 (1415/1559-1525/1605)  24.4 W21-W32 N41-N46
      5  546 1422 (N17E38 L= 85)   604    175        1 ( 484/1375- 589/1493)  83.7 E41-E35 N13-N22
      6 1458 1425 (N32W24 L=148)   439    111        0 (1431/1377-1476/1468) -90.0 W22-W25 N28-N34
      7 1168 1242 (N15W06 L=129)   835    234        2 (1101/1181-1203/1351)  85.2 W02-W08 N10-N22
      8  911 1305 (N16E11 L=112)   292    112        1 ( 855/1291- 968/1322)   8.1 E14-E07 N14-N17
      9  671 1236 (N08E25 L= 98)   292     94        2 ( 635/1207- 703/1273)  78.7 E26-E22 N05-N10
     10 1253 1050 (N05W13 L=137)   503    141        1 (1212/ 974-1281/1109) -54.5 W10-W16 N00-N08
     11  377 1011 (S10E43 L= 81)   306     86        0 ( 340/ 990- 409/1024)  -9.5 E45-E39 S10-S09
     12 1621  921 (N02W40 L=164)   610     85        2 (1593/ 878-1642/ 969) -68.2 W37-W42 N00-N05
     13  809  895 (S11E11 L=112)   318    139        0 ( 741/ 879- 866/ 921)  23.2 E15-E08 S11-S08
     14 1165  800 (S11W12 L=135)   739    139        2 (1125/ 735-1201/ 877) -71.6 W08-W14 S15-S06
     15 1607  754 (S08W43 L=166)   267     67        0 (1575/ 719-1633/ 786)  68.2 W40-W44 S10-S05
     16 1710  694 (S10W53 L=177)   447     85        2 (1677/ 660-1734/ 736)  68.2 W51-W55 S12-S07
     17  666  566 (S34E18 L=105)   295     61        0 ( 637/ 547- 691/ 587)  36.9 E20-E16 S35-S32
  

    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