Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1345 1229 (N17W16 L=283)  3555    781        4 (1053/ 928-1575/1520) -30.6 E08-W36 N03-N29
      2  215 1035 (S15E58 L=209)   724     71        1 ( 193/ 986- 237/1072)   9.3 E60-E55 S17-S12
      3  139 1001 (S18E67 L=200)  1699    221        4 ( 115/ 943- 172/1070)  63.4 E70-E63 S21-S13
      4  258  978 (S18E53 L=214)   520    102        1 ( 235/ 927- 276/1025)  74.1 E54-E50 S21-S14
      5  553  974 (S13E28 L=238)   214     49        0 ( 542/ 946- 564/1010) -76.0 E29-E27 S14-S10
      6  335  923 (S20E45 L=222)   279     65        0 ( 324/ 882- 345/ 960) -78.7 E45-E43 S22-S17
      7 1344  760 (S11W25 L=291)   696    197        0 (1245/ 720-1439/ 793)  32.5 W19-W30 S14-S07
      8  587  763 (S25E24 L=243)   462    123        1 ( 526/ 730- 640/ 786)   0.0 E28-E19 S26-S23
      9 1081  692 (S20W10 L=277)   780    182        1 (1000/ 671-1182/ 727)   4.8 W04-W16 S21-S18
     10 1335  302 (S41W44 L=310)  2872    479        3 (1117/ 274-1503/ 339)  14.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