Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  935 1725 (N33E26 L=310)   980    218        1 ( 841/1660-1008/1777) -13.2 E35-E18 N29-N35
      2 1140 1532 (N25E06 L=330)   212     72        1 (1126/1495-1153/1571) -59.0 E07-E03 N22-N27
      3  985 1485 (N18E13 L=323)   238     63        0 ( 951/1474-1031/1500)  11.3 E15-E10 N17-N18
      4  179 1403 (S00E74 L=262)   684     39        0 ( 162/1385- 200/1427) -15.9 E79-E70 S01-N01
      5 1722 1314 (N28W39 L= 15)   293     61        0 (1690/1289-1753/1330)  39.8 W36-W42 N25-N30
      6 1335 1235 (N12W13 L=349)   225     46        2 (1308/1213-1353/1260)  90.0 W10-W13 N10-N13
      7  443 1161 (S11E39 L=298)   422    139        0 ( 424/1092- 463/1217) -81.9 E41-E37 S15-S08
      8 1142 1160 (N03W03 L=339)   558    261        0 (1010/1137-1250/1193)  32.5 E03-W08 N00-N07
      9 1724 1135 (N17W41 L= 17)   277     56        0 (1703/1114-1745/1165)  53.1 W39-W42 N15-N19
     10 1923  939 (N13W67 L= 43)  1045    106        1 (1908/ 879-1938/ 998) -74.1 W64-W70 N09-N16
     11  357  954 (S25E42 L=294)   421     71        3 ( 320/ 922- 393/ 978)  -9.5 E45-E39 S26-S23
     12 1665  846 (S00W43 L= 19)   512    117        1 (1616/ 825-1742/ 867)  11.3 W39-W49 S01-N01
     13  239  667 (S45E60 L=276)  1848    216        3 ( 154/ 626- 328/ 717)   0.0 E82-E47 S46-S41
     14 1591  659 (S12W45 L= 21)   259     49        0 (1564/ 634-1616/ 679)  45.0 W43-W46 S14-S10
     15 1360  575 (S23W33 L=  9)   873    298        2 (1256/ 553-1527/ 607)  17.5 W25-W44 S26-S18
     16 1517  427 (S27W53 L= 29)  2299    436        3 (1363/ 364-1745/ 510)  41.0 W44-W68 S34-S14
  

    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