Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1570 1650 (N30W53 L= 92)  1286    153        1 (1539/1577-1601/1725)  54.0 W49-W59 N24-N35
      2 1322 1599 (N31W30 L= 69)   672    152        0 (1250/1580-1402/1613) -16.7 W25-W35 N28-N32
      3 1702 1535 (N21W59 L= 98)   460     44        1 (1675/1520-1736/1550)  -9.5 W55-W62 N19-N21
      4 1799 1456 (N15W65 L=104)  1520     49        1 (1774/1403-1827/1510)  33.7 W61-W72 N11-N17
      5 1376 1343 (N14W26 L= 65)   357    108        1 (1325/1319-1431/1369) -33.7 W23-W29 N11-N15
      6  763 1225 (N15E13 L= 26)   245     92        0 ( 733/1192- 796/1267) -59.0 E14-E11 N13-N18
      7 1775 1106 (S05W51 L= 90)  1282    139        2 (1708/1068-1833/1159)   8.7 W44-W57 S07-S02
      8 1313 1045 (S03W17 L= 56)   538    157        1 (1286/ 971-1357/1109)  45.0 W14-W20 S06-N00
      9 1488  944 (S11W27 L= 66)   561    202        0 (1471/ 861-1507/1048)  79.7 W25-W28 S15-S04
     10 1641 1002 (S10W39 L= 77)   458    119        0 (1595/ 980-1705/1037) -26.6 W35-W43 S11-S07
     11  929  713 (S17E10 L= 29)  3783    920       13 ( 739/ 361-1034/1036)  76.6 E23-E05 S39-N03
     12  184  704 (S07E68 L=331)  1064     77        3 ( 150/ 666- 209/ 746)  39.8 E73-E64 S09-S04
     13 1386  667 (S26W18 L= 57)  1201    281        2 (1254/ 618-1500/ 734)   6.3 W08-W26 S28-S23
     14  360  684 (S11E49 L=349)   615    108        1 ( 322/ 644- 391/ 717) -59.0 E52-E45 S13-S08
     15  465  663 (S13E42 L=357)   379     79        2 ( 441/ 613- 487/ 691)  68.2 E43-E39 S16-S11
     16  400  474 (S24E56 L=343)  2739    282        4 ( 296/ 406- 510/ 565)  11.6 E75-E41 S27-S19
     17  717  298 (S41E40 L=359)   510     57        1 ( 688/ 277- 749/ 323)  -9.5 E42-E36 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