Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  463 1559 (N18E49 L= 32)  1461    257        3 ( 410/1492- 541/1660) -30.1 E58-E39 N13-N25
      2 1116 1551 (N27E04 L= 77)   942    272        3 (1009/1461-1219/1615) -20.6 E12-W04 N23-N30
      3  495 1291 (N02E37 L= 43)   316     75        0 ( 459/1261- 525/1318) -18.4 E40-E34 N00-N03
      4  486 1128 (S08E35 L= 46)   906    292        3 ( 422/1051- 559/1243)  65.2 E39-E30 S13-N00
      5 1807 1140 (N16W52 L=133)   320     53        0 (1791/1114-1824/1164) -45.0 W49-W53 N13-N17
      6  253  960 (S20E54 L= 27)   384     66        1 ( 239/ 922- 268/ 997) -63.4 E55-E52 S22-S18
      7  660  940 (S17E21 L= 60)   271     52        0 ( 634/ 923- 689/ 966) -33.7 E22-E18 S17-S15
      8  533  797 (S27E29 L= 52)  1005    299        2 ( 443/ 703- 597/ 927)  55.5 E36-E24 S34-S18
      9  705  677 (S32E15 L= 66)  2123    554        6 ( 601/ 511- 833/ 924)  65.9 E23-E09 S45-S15
     10 1614  657 (S16W46 L=127)   284     50        0 (1585/ 640-1638/ 677)  45.0 W44-W48 S17-S14
     11  875  622 (S33E03 L= 78)   510    161        0 ( 797/ 607- 948/ 638)  21.8 E08-W02 S34-S30
  

    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