Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1009 1781 (N41E21 L= 98)   298     52        1 ( 990/1756-1031/1799)  90.0 E22-E19 N39-N43
      2  867 1420 (N13E17 L=101)   451    138        1 ( 836/1363- 914/1485) -66.8 E20-E13 N09-N16
      3  501 1335 (N02E38 L= 80)   271     58        1 ( 486/1305- 520/1365) -45.0 E40-E36 N00-N03
      4  399 1272 (S03E45 L= 74)  1129    217        2 ( 318/1180- 476/1363)  56.3 E49-E40 S09-N03
      5 1509 1182 (N12W26 L=144)  1161    303        3 (1461/1031-1538/1337) -58.0 W20-W30 N04-N20
      6  174 1167 (S12E64 L= 55)  2337    335        5 ( 130/1069- 224/1248) -70.0 E69-E59 S17-S06
      7 1734 1037 (N10W45 L=164)   490    105        0 (1717/ 980-1749/1092) -81.9 W43-W46 N06-N13
      8 1889  961 (N10W62 L=181)  2367    459        6 (1839/ 875-1943/1053)  90.0 W57-W69 N04-N15
      9  669  852 (S24E18 L=100)   226     41        1 ( 642/ 837- 696/ 867)  26.6 E20-E16 S24-S22
     10 1240  480 (S34W28 L=147)   365    112        1 (1192/ 466-1300/ 506)   7.1 W24-W32 S35-S32
  

    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