Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  595 1803 (N45E52 L=264)  1257    128        5 ( 544/1776- 648/1829)   4.8 E63-E43 N42-N47
      2  847 1733 (N40E22 L=295)   375    108        1 ( 802/1708- 898/1751)  23.2 E25-E17 N37-N41
      3  648 1593 (N28E32 L=284)   288     64        1 ( 633/1564- 666/1626) -78.7 E34-E30 N25-N30
      4 1691 1557 (N37W50 L=  7)  1438    182        1 (1615/1491-1749/1612) -32.5 W43-W54 N33-N41
      5  429 1541 (N23E49 L=268)   363     46        2 ( 410/1518- 447/1566)  71.6 E50-E46 N21-N24
      6  787 1472 (N20E19 L=297)   210     67        0 ( 761/1441- 813/1493)  63.4 E20-E17 N18-N22
      7  948 1336 (N13E07 L=309)   691    292        1 ( 898/1239-1022/1470) -60.3 E12-E02 N07-N21
      8 1301 1040 (S02W16 L=333)   905    240        1 (1243/ 941-1375/1122)  61.4 W13-W20 S07-N04
      9 1732 1005 (N01W47 L=  4)   277     11        0 (1711/ 989-1751/1027)  33.7 W45-W49 N00-N02
     10 1623  968 (S02W39 L=356)   342     89        0 (1596/ 924-1646/1006)  90.0 W37-W40 S05-N00
     11 1554  679 (S20W39 L=356)   261    101        0 (1499/ 667-1602/ 698)  23.2 W35-W43 S21-S18
     12  918  551 (S35E03 L=314)  1515    371        3 ( 761/ 477-1046/ 685)  41.6 E14-W05 S41-S25
     13  340  630 (S34E51 L=265)   534    104        1 ( 324/ 582- 362/ 682)  49.4 E54-E47 S37-S30
     14  520  608 (S35E34 L=283)   392     44        2 ( 486/ 589- 556/ 625)   0.0 E36-E30 S35-S33
     15  965  390 (S47W03 L=320)   322     83        0 ( 919/ 372- 996/ 411)  18.4 E01-W05 S48-S44
     16 1501  315 (S44W57 L= 14)   487     49        0 (1466/ 302-1531/ 330)  15.9 W52-W61 S45-S43
  

    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