Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1624 1593 (N22W59 L=182)   529    128        1 (1600/1544-1640/1627)  33.7 W55-W61 N19-N23
      2 1759 1221 (S02W53 L=176)   379    114        2 (1720/1200-1798/1243) -26.6 W50-W56 S03-N00
      3  160 1207 (N35E62 L= 61)   626    204        0 ( 126/1179- 199/1227) -14.0 E65-E57 N32-N35
      4 1826 1088 (S11W56 L=179)   641    144        3 (1811/1011-1846/1161)  71.6 W53-W57 S15-S06
      5 1699  704 (S30W37 L=161)  2008    366        3 (1553/ 572-1825/ 842)  16.5 W23-W50 S34-S25
      6 1148  627 (S19E02 L=121)   226    108        0 (1096/ 614-1202/ 638) -21.8 E05-W00 S20-S17
      7  826  538 (S17E24 L= 99)   306     71        4 ( 798/ 500- 851/ 580)  80.5 E25-E21 S19-S13
      8  727  547 (S14E30 L= 93)   254     51        1 ( 699/ 531- 752/ 563) -18.4 E31-E28 S15-S13
      9  950  481 (S23E19 L=105)   512    135        0 ( 894/ 452-1021/ 528) -45.0 E20-E14 S26-S19
     10  705  432 (S21E36 L= 87)   817    216        1 ( 628/ 350- 747/ 520) -90.0 E37-E33 S27-S14
     11  416  460 (S14E57 L= 66)   847    137        2 ( 365/ 412- 473/ 498) -45.0 E61-E54 S18-S11
     12  534  399 (S20E50 L= 73)  1408    355        7 ( 487/ 334- 592/ 457)  31.6 E58-E44 S24-S16
  

    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