Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  935 1557 (N38W11 L=162)   714    165        1 ( 887/1481- 992/1617)  24.4 W06-W17 N34-N41
      2 1222 1539 (N28W28 L=179)   364     63        3 (1201/1507-1244/1570)  76.0 W26-W29 N26-N30
      3 1357 1512 (N23W35 L=186)  1021    183        4 (1273/1475-1438/1542) -41.2 W29-W38 N18-N25
      4  296 1211 (N34E45 L=106)   349     91        2 ( 268/1189- 327/1234) -38.7 E47-E42 N32-N36
      5  470 1160 (N27E31 L=120)   589    182        0 ( 399/1136- 558/1195) -13.0 E36-E23 N25-N29
      6  695 1093 (N18E17 L=134)   750    152        1 ( 619/1068- 772/1139) -29.1 E21-E12 N15-N20
      7 1347 1063 (S01W19 L=170)  1049    229        4 (1238/1036-1451/1101) -18.4 W12-W24 S04-N01
      8 1909 1048 (S20W62 L=213)   512     97        2 (1889/1025-1930/1070)  18.4 W59-W64 S21-S18
      9 1309  983 (S04W15 L=165)   369    126        1 (1255/ 946-1373/1000)  -7.1 W10-W18 S06-S03
     10 1525  898 (S15W26 L=176)   533    109        1 (1477/ 866-1575/ 939)   0.0 W22-W29 S17-S13
     11 1239  739 (S16W04 L=155)   713    188        0 (1208/ 668-1279/ 837)  79.7 W02-W06 S21-S10
     12  919  597 (S16E18 L=133)   691    150        0 ( 823/ 565- 980/ 627) -41.2 E22-E14 S18-S11
     13  596  435 (S18E44 L=106)   670    135        3 ( 540/ 403- 652/ 466) -12.5 E48-E39 S19-S15
     14  503  416 (S17E53 L= 98)   472    110        0 ( 462/ 378- 552/ 459) -66.8 E54-E50 S20-S13
     15  642  349 (S24E46 L=105)   635    103        2 ( 604/ 305- 676/ 393)  71.6 E47-E43 S27-S21
     16  905  171 (S45E43 L=108)   944    134        2 ( 883/ 129- 929/ 209)  38.7 E48-E38 S49-S41
  

    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