Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1414 1734 (N55W07 L=159)  1404    156        3 (1370/1681-1454/1782) -37.9 W01-W12 N50-N59
      2 1082 1703 (N39E18 L=134)   666    119        0 (1025/1676-1142/1729)   6.3 E22-E12 N37-N40
      3  763 1532 (N19E30 L=121)  1292    242        5 ( 712/1455- 827/1650) -35.5 E37-E23 N14-N25
      4 1540 1603 (N48W22 L=174)   527     77        1 (1507/1582-1577/1625)   0.0 W18-W26 N46-N49
      5  323 1531 (N09E66 L= 86)  1674    205        2 ( 266/1469- 356/1593) -83.7 E75-E62 N04-N13
      6  904 1373 (N13E16 L=136)   529    200        0 ( 851/1317- 973/1472)  85.2 E17-E14 N08-N20
      7  994 1317 (N12E09 L=142)   255    103        0 ( 964/1272-1017/1372) -39.8 E12-E06 N09-N14
      8 1528 1228 (N21W25 L=177)   311     97        0 (1483/1204-1573/1245)   0.0 W22-W29 N20-N22
      9 1800 1118 (N23W51 L=202)   772    126        1 (1752/1093-1846/1147)  -5.7 W45-W56 N21-N25
     10  570  874 (S24E24 L=127)   282    114        0 ( 529/ 857- 627/ 903)  15.9 E27-E20 S24-S21
     11  666  831 (S24E17 L=135)   214     90        0 ( 626/ 807- 707/ 854)  38.7 E19-E14 S26-S21
     12  920  777 (S20W00 L=152)   202     60        0 ( 886/ 767- 951/ 791)  33.7 E01-W02 S21-S19
     13 1311  681 (S16W26 L=178)   611    157        2 (1276/ 608-1346/ 759) -78.7 W24-W28 S20-S10
     14 1528  622 (S13W43 L=194)   395     64        1 (1500/ 592-1562/ 645)  63.4 W40-W44 S15-S11
     15  745  595 (S36E06 L=145)   311     52        1 ( 717/ 567- 773/ 617)  45.0 E07-E04 S38-S34
     16  867  334 (S51W13 L=165)   363     45        1 ( 843/ 309- 892/ 362)  78.7 W11-W15 S53-S48
  

    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