Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  952 1641 (N48W07 L=326)   438     67        0 ( 925/1617- 989/1660)  15.9 W03-W10 N46-N49
      2 1042 1543 (N39W11 L=331)   672    135        1 (1013/1493-1088/1599) -77.5 W09-W13 N34-N43
      3 1455 1476 (N26W39 L=358)   266     71        0 (1416/1448-1491/1503) -51.3 W36-W40 N23-N28
      4  947 1439 (N33W02 L=321)   269     69        1 ( 929/1404- 965/1469) -90.0 W00-W02 N30-N35
      5  352 1093 (N18E46 L=274)   878    203        1 ( 271/1033- 436/1151) -26.6 E54-E38 N14-N22
      6 1573  801 (S16W32 L=352)   594    211        0 (1488/ 745-1634/ 876)  22.6 W25-W37 S18-S12
      7  290  842 (N03E54 L=265)   433     53        1 ( 263/ 815- 320/ 868) -21.8 E56-E51 N01-N05
      8  800  763 (S06E18 L=302)   823    184        2 ( 759/ 671- 840/ 835)  48.4 E21-E13 S11-S02
      9 1294  783 (S12W13 L=333)   377     79        1 (1270/ 735-1318/ 818)  59.0 W11-W14 S15-S10
     10 1680  736 (S23W41 L=  1)   865    190        1 (1656/ 680-1726/ 812)  90.0 W39-W46 S27-S18
     11 1304  546 (S28W11 L=331)   549    198        0 (1272/ 457-1344/ 629)  61.4 W08-W14 S34-S23
     12 1018  379 (S36E12 L=308)   302     30        1 ( 999/ 354-1036/ 402)  45.0 E13-E10 S37-S33
  

    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