Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  416 1603 (N15E61 L= 18)   482     78        0 ( 395/1569- 445/1636) -16.7 E65-E55 N13-N16
      2  287 1524 (N07E71 L=  8)  1629    114        4 ( 257/1477- 320/1572) -71.6 E77-E67 N04-N10
      3 1327 1509 (N33W05 L= 84)   430     74        1 (1294/1476-1363/1539) -26.6 W02-W08 N31-N34
      4 1894 1017 (N21W65 L=144)  2625    291        2 (1828/ 908-1953/1127)  35.8 W57-W85 N12-N30
      5  646 1015 (S14E22 L= 57)   228     76        0 ( 600/1003- 682/1028)   9.5 E25-E19 S15-S12
      6  316  820 (S34E46 L= 33)   314    122        0 ( 277/ 787- 370/ 838)  20.0 E50-E39 S37-S32
      7  557  702 (S35E22 L= 57)  1953    358        5 ( 452/ 599- 682/ 833)  57.0 E31-E15 S44-S23
      8 1572  596 (S13W48 L=127)  1389    384        3 (1485/ 481-1659/ 691)  82.9 W41-W52 S20-S04
      9 1226  462 (S31W31 L=110)  1516    356        3 (1099/ 387-1374/ 538)  59.5 W23-W37 S39-S22
  

    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