Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  469 1448 (N46E33 L= 55)   625    173        0 ( 399/1419- 550/1470)   0.0 E41-E24 N44-N46
      2 1010 1041 (N09E00 L= 87)   272     97        1 ( 985/1006-1029/1090)  45.0 E02-W01 N06-N10
      3  285  922 (N15E53 L= 35)  1462    233        1 ( 237/ 836- 336/1030) -65.2 E57-E49 N09-N22
      4  193  894 (N15E65 L= 23)  2044    226        5 ( 134/ 831- 241/ 993) -78.7 E74-E59 N10-N20
      5 1501  826 (S16W24 L=112)   827    188        1 (1454/ 743-1562/ 904)  33.7 W20-W29 S20-S12
      6 1116  639 (S17E04 L= 84)  1147    423        1 (1058/ 486-1201/ 795)  50.5 E10-W04 S26-S09
      7 1708  647 (S35W40 L=128)   404     26        0 (1689/ 625-1724/ 673) -63.4 W37-W41 S36-S32
      8 1312  611 (S24W08 L= 95)   801    214        0 (1231/ 552-1405/ 672)  13.0 W01-W14 S26-S22
      9  280  620 (S03E64 L= 24)   636     90        1 ( 260/ 584- 308/ 661)  78.7 E65-E61 S05-N00
     10 1438  508 (S35W15 L=102)   472     99        0 (1401/ 475-1476/ 542)  18.4 W11-W17 S36-S33
     11  740  278 (S34E44 L= 44)  2234    289        1 ( 685/ 165- 800/ 417)  43.7 E60-E37 S44-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