Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1272 1271 (N14W21 L=278)  1410    388        1 (1144/1200-1445/1299) -14.0 W11-W31 N09-N16
      2 1817 1004 (S18W51 L=308)  1129    220        1 (1757/ 970-1881/1037) -12.1 W44-W58 S20-S15
      3  600  983 (N15E26 L=232)  2061    436        6 ( 425/ 930- 776/1036) -21.4 E38-E14 N08-N19
      4 1692  811 (S25W36 L=294)   776    165        0 (1646/ 734-1722/ 893)  37.9 W31-W40 S28-S20
      5  302  735 (N07E56 L=202)  1777    278        4 ( 261/ 624- 358/ 880)  60.3 E60-E49 N00-N14
      6 1240  748 (S14W05 L=262)   418    195        0 (1215/ 664-1283/ 839)  52.1 W02-W09 S19-S10
      7 1149  778 (S10W01 L=258)   587    190        2 (1104/ 699-1209/ 838)  32.0 E03-W05 S14-S07
      8 1674  634 (S36W33 L=290)   599    133        1 (1621/ 609-1716/ 654) -32.0 W28-W36 S39-S33
      9 1244  497 (S30E02 L=256)   694    154        0 (1217/ 426-1283/ 558)  66.0 E04-W02 S35-S26
     10  678  397 (S21E41 L=217)  2536    410        3 ( 577/ 245- 756/ 527)  77.7 E47-E37 S33-S10
     11  495  478 (S12E50 L=208)   372     63        0 ( 477/ 440- 511/ 505)  53.1 E51-E47 S14-S10
     12  560  329 (S23E54 L=204)   629    205        1 ( 537/ 274- 596/ 389)  66.0 E56-E48 S28-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