Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1336 1234 (N17W15 L=283)  3464    768        4 (1046/ 935-1565/1523) -31.2 E08-W35 N03-N30
      2  212 1048 (S14E59 L=209)   552     61        1 ( 187/1020- 235/1072) -50.2 E61-E56 S15-S12
      3  137  998 (S18E68 L=200)  1189    164        2 ( 114/ 948- 161/1052)  90.0 E71-E64 S21-S14
      4  254  978 (S18E53 L=214)   476    106        1 ( 232/ 927- 273/1026)  74.1 E54-E51 S21-S14
      5 1589  855 (S01W39 L=306)   279     50        1 (1560/ 832-1619/ 875)   0.0 W36-W41 S01-N00
      6 1368  775 (S09W26 L=293)   458     90        0 (1318/ 757-1420/ 795)  15.9 W22-W29 S10-S08
      7  576  767 (S25E24 L=243)   467    127        1 ( 515/ 733- 631/ 790)   0.0 E28-E20 S26-S23
      8 1073  694 (S20W09 L=277)   769    175        1 ( 997/ 672-1170/ 731)   5.2 W04-W15 S21-S18
      9 1335  304 (S40W43 L=311)  3110    504        3 (1112/ 276-1502/ 339)  12.0 W25-W58 S44-S35
  

    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