Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1126 1837 (N53E25 L=154)  1639    165        1 (1041/1802-1191/1861)  19.7 E31-E17 N48-N56
      2  830 1784 (N38E41 L=138)   408     48        1 ( 805/1757- 853/1808) -15.9 E44-E37 N35-N39
      3 1324 1701 (N48E01 L=178)   426     54        1 (1292/1677-1356/1722)  53.1 E02-W01 N45-N50
      4  463 1616 (N17E57 L=122)   886    133        1 ( 438/1578- 509/1666) -50.2 E62-E51 N14-N20
      5  506 1500 (N11E46 L=133)   445     69        0 ( 479/1469- 538/1532)  78.7 E48-E45 N08-N13
      6 1755 1407 (N41W46 L=225)   490     89        1 (1705/1395-1793/1421)  18.4 W39-W51 N38-N42
      7 1154 1394 (N20E02 L=177)   409    155        1 (1073/1365-1215/1415)  12.5 E07-W02 N18-N21
      8 1754 1060 (N18W47 L=226)   573     63        1 (1729/1018-1769/1099) -78.7 W44-W48 N15-N20
      9  400  998 (S21E39 L=140)   479    165        0 ( 374/ 932- 415/1079) -71.6 E40-E37 S24-S15
     10  297 1018 (S22E49 L=130)   292     48        1 ( 276/1002- 324/1034)  36.9 E51-E46 S23-S20
     11  921  928 (S12E03 L=176)   567    200        1 ( 878/ 835- 946/1012) -84.8 E04-E02 S18-S07
     12  524  914 (S23E28 L=151)   305     58        0 ( 491/ 893- 551/ 937)  26.6 E30-E26 S23-S21
     13 1188  811 (S12W15 L=194)   270     53        0 (1159/ 788-1218/ 830)   0.0 W12-W16 S12-S10
     14  451  739 (S35E32 L=147)   575    137        1 ( 408/ 698- 498/ 799)  52.1 E35-E28 S39-S30
     15  992  658 (S26W08 L=187)   246     85        0 ( 956/ 639-1031/ 685)   0.0 W05-W11 S26-S24
     16 1665  488 (S16W63 L=242)   440     58        1 (1635/ 464-1692/ 513)  53.1 W61-W64 S17-S13
     17 1496  423 (S25W52 L=231)   666    124        0 (1438/ 372-1547/ 477)  77.5 W51-W54 S29-S20
     18  699  404 (S52E04 L=175)   309     21        1 ( 676/ 385- 721/ 426) -45.0 E06-E01 S53-S50
  

    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