Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  888 1844 (N41E38 L=176)  1230    299        1 ( 852/1773- 927/1916) -35.2 E47-E30 N36-N48
      2  795 1532 (N16E27 L=188)   533    105        0 ( 754/1481- 842/1583)  82.9 E28-E25 N12-N20
      3  743 1335 (N04E24 L=191)   556    116        0 ( 690/1287- 774/1401)  74.1 E25-E22 N00-N07
      4 1692 1294 (N27W35 L=250)   866    129        1 (1653/1224-1735/1363) -40.6 W31-W39 N23-N30
      5  883 1261 (N02E14 L=201)   921    264        1 ( 738/1225-1012/1281)  14.0 E22-E06 N00-N05
      6 1046 1075 (S04E00 L=215)   911    263        2 (1008/ 952-1077/1222) -53.1 E05-W04 S10-N03
      7  908  981 (S12E05 L=209)   963    239        2 ( 800/ 939-1031/1020)  28.6 E11-W01 S17-S08
      8  623  641 (S40E16 L=198)  1768    421        2 ( 444/ 621- 823/ 675)  16.2 E32-E01 S44-S34
      9 1211  611 (S25W23 L=238)  1108    318        1 (1056/ 584-1370/ 644)  29.1 W13-W31 S29-S19
  

    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