Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1690 1481 (N15W59 L=273)  1233    183        3 (1647/1391-1730/1557)  70.0 W54-W64 N08-N19
      2  323 1421 (N46E51 L=163)   423     83        2 ( 291/1402- 362/1443) -10.3 E55-E44 N44-N47
      3 1183 1402 (N25W19 L=233)   211     60        0 (1164/1374-1201/1432) -90.0 W18-W20 N23-N26
      4 1456 1262 (N10W31 L=245)   325     88        0 (1433/1225-1487/1307) -78.7 W30-W32 N07-N12
      5  435 1266 (N34E36 L=178)   290    114        0 ( 388/1226- 481/1287) -26.6 E40-E31 N30-N35
      6 1407 1030 (S01W22 L=236)  1311    409        1 (1289/ 860-1517/1136)  21.8 W11-W31 S08-N03
      7 1726  973 (S14W43 L=257)   506    124        0 (1705/ 910-1767/1024)  35.5 W40-W47 S16-S11
      8 1766  954 (S16W47 L=261)   860    182        1 (1726/ 865-1811/1020)  35.0 W41-W51 S20-S12
      9  906  753 (S06E13 L=201)   630    137        0 ( 889/ 678- 931/ 853)  73.3 E14-E10 S10-N00
     10  967  752 (S07E09 L=205)   267     97        1 ( 933/ 715-1007/ 796) -63.4 E10-E07 S09-S03
     11 1131  729 (S12E01 L=213)   572    230        2 (1076/ 645-1204/ 786)  29.1 E05-W04 S16-S09
     12  398  751 (N04E46 L=168)   307     67        0 ( 363/ 738- 431/ 770)   8.1 E49-E42 N02-N04
     13  523  451 (S17E48 L=166) 10071      0       26 ( 393/ 235- 674/ 693)  41.6 E77-E33 S33-S01
     14  670  594 (S10E31 L=182)   355     53        1 ( 633/ 569- 697/ 618) -33.7 E33-E30 S12-S08
     15 1444  551 (S32W16 L=229)   267     55        0 (1423/ 521-1462/ 584)  63.4 W13-W17 S33-S29
     16  950  529 (S20E16 L=198)   329     92        0 ( 918/ 494- 981/ 576)  33.7 E19-E13 S22-S17
     17  735  337 (S29E39 L=175)  1303    275        0 ( 647/ 252- 843/ 441)  21.0 E52-E26 S33-S23
  

    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