Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1271 1663 (N39W05 L=238)  1371    401        0 (1109/1575-1438/1754) -14.0 E09-W19 N35-N43
      2  361 1626 (N22E65 L=168)  1631    114        2 ( 318/1583- 399/1676)  59.0 E71-E58 N19-N26
      3  502 1121 (S08E34 L=199)   335    101        0 ( 478/1068- 533/1154)  71.6 E35-E32 S11-S05
      4  498  888 (S22E32 L=201)   802    255        2 ( 440/ 810- 543/ 974)  74.7 E36-E28 S27-S16
      5  141  833 (S28E72 L=161)   658     51        1 ( 110/ 819- 173/ 853)   3.8 E80-E65 S28-S26
      6 1288  796 (S15W20 L=253)   453    137        1 (1222/ 756-1339/ 838)  49.4 W16-W22 S18-S11
      7  745  778 (S25E14 L=219)   244    111        0 ( 699/ 748- 784/ 814)  45.0 E16-E11 S27-S22
      8 1271  660 (S23W22 L=255)  2838    817        4 ( 950/ 547-1642/ 811)  35.0 W03-W43 S35-S07
      9 1704  637 (S15W55 L=288)   723     98        0 (1662/ 596-1737/ 678)  63.4 W53-W57 S18-S12
     10  833  583 (S36E05 L=228)   486    128        0 ( 775/ 553- 887/ 626)  -7.1 E09-E01 S37-S34
  

    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